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A Further Study of the y-Radiation from Polonium 


By B. ZAJAC, E. BRODA anp N. FEATHER 
Department of Natural Philosophy, University of Edinburgh 


MS. received 7 February 1948 


ABSTRACT. Absorption experiments using lead, gold and tungsten have shown that, in 
addition to the known y-radiation of 0°77 Mev. energy, polonium emits other (soft) radiations 
of which the most intense has a quantum energy of 84+4 kev. The intensity of this 
radiation is of the same order of magnitude as that of the hard radiation (roughly one 
quantum per 10° disintegrations). Experiments by the recoil method indicate that the 
emission of the polonium y-radiation is not delayed by more than 107! sec. The results of 
Chang (1946) concerning the fine-structure of the a-particles of polonium remain uncor- 
related with all other experimental evidence. Whilst present information regarding y-ray 
energies and intensities may be reasonably explained, satisfactory explanation of the 
a-particle fine structure “‘ data’’ appears as remote as ever. 


Sila LINPINIRKOMDMOKC AMICON 

NTEREST in the y-radiation of polonium has recently been revived by the 
[ pivtcesn by Chang (1946) of quite unexpected evidence for a many- 

lined fine structure in the «-particle spectrum of this body. ‘Two experi- 
mental studies of the polonium y-radiation have since been described (de 
Benedetti and Kerner 1947, Siegbahn and Slatis 1947) and a critical discussion 
of Chang’s results in relation to earlier information on the y-radiation (Bothe 
and Becker 1930, Webster 1932, Bothe 1935, 1936) has been given by one of us 
(Feather 1946). The experiments now to be described were, like the others 
already mentioned, prompted by a consideration of Chang’s results; they deal 
more particularly with the question whether a true nuclear radiation of quantum 
energy less than that of the now (Siegbahn and Slatis 1947) well-established 
“line” of 0:77 Mev. energy is in fact emitted, and secondly whether or not the 
emission of the polonium y-radiation is measurably delayed in relation to the 
primary «-disintegration process. It should be remarked that the absorption 
measurements of Siegbahn and Slatis (1947) and de Benedetti and Kerner (1947) 
would appear to rule out the possibility of any significant emission of quantum 
radiation in the energy range 0-1 to 0-7 Mev. (or even of somewhat lower energy) ; 
on the other hand the only possible remaining way of explaining Chang’s results 
on the basis of nuclear excitation followed by y-ray emission is to postulate the 
relatively frequent emission of quanta of low energy (Feather 1946). It may be 
said in anticipation that no adequate low-energy quantum emission could be 
established in our experiments, although we have shown, in apparent contra- 
diction of a statement of de Benedetti and Kerner, that weak nuclear radiations 
of about 84 kev. energy are present. No evidence that any appreciable fraction 
of the polonium y-radiation is delayed by as much as 10~!sec. after «-emission 
was obtained. It will be recalled that the possibility of such delayed emission 
was discussed in an attempt to interpret (Feather 1946) the reported observation 


> 
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of y-radiation of quantum energy 0:2 to 0-43 Mev. (Bothe 1935, 1936) in terms 
of excitation of nuclear states of corresponding energy only to the extent of a 
few times per million disintegrations. Now that the emission of such y-radiation 
has been effectively disproved (de Benedetti and Kerner 1947, Siegbahn and 
Slatis 1947) the reason behind this suggestion obviously disappears. 


§2. ABSORPTION EXPERIMENTS 

Absorption of the y-rays in lead was investigated over a range of thickness 
from zero to 3:4gm/cm2, and in a separate experiment a careful comparison 
of the absorption in lead, gold and tungsten was made up to thicknesses of 
450mg/cm?. The first experiment extended over a period of 15 days, and the 
second (comparison) experiment over 48 days. No significant discrepancy 
appeared in the process of correlating all the observations in these experiments 
by simple correction for decay, taking A=5-74x 10-8sec"?, The strength of 
the source at the beginning of the absorption measurements was about 6mc., 
and a rough estimate based on f-particle counting showed that at this stage 
the radium (D +E) contamination was responsible for about 10-° disintegrations 
per disintegration of polonium. ‘The source was deposited by evaporation 
of an HCI solution, drop by drop, on a polythene disc 50 mg/cm? thick over an 
area 0-5 cm. in diameter. 

A thin (end) windowed counter was employed for detection of the y-radiations, 
the polythene disc carrying the source being supported on wood at a distance 
of approximately 2cm. from the counter window, which was 1:-4cm. in diameter 
and 2:6mg/cm? thick. ‘The window material was mica. A standard argon— 
alcohol filling was employed and the counter pulses were registered by a “scale 
of 32” counting set. The average background counting rate of the (unshielded) 
counter was about 26min-!. ‘Throughout the experiments aluminium absorbers 
of 480 mg/cm? thickness were placed directly over the source to cut out the 
primary f-particles from the radium E “impurity”’, and a further absorber of 
aluminium of 30mg/cm? was invariably inserted immediately in front of the 
counter window so that the secondary electrons entering the counter should always 
originate in the same material. All absorbers were square foils of 3-2 cm. side. 

The lead absorption curve obtained with this experimental arrangement 
is given in figure 1, Analysis of this curve by ‘‘least square’? and graphical 
methods indicates component radiations as follows : 


hard component (./p) = 0-086 cm?/gm., intensity 75%, 
soft component (./p) = 2-2 cm?/gm., intensity 16%, 
more absorbable radiation, intensity 9°. 


The value of the mass absorption coefficient of the hard component is to be 
compared with 0-0815 cm?/gm. obtained by de Benedetti and Kerner in measure- 
ments up to 62gm/cm?, and 0-081 cm?/gm. deduced by Siegbahn and Slatis 
from observations over 33 gm/cm? absorber thickness *. It should be entirely | 
unnecessary to point out that our value cannot claim anything approaching the © 
independent weight of either of these other two determinations, made with much 
stronger sources under better geometrical conditions, but the numerical agreement 
is at least satistactory. On the other hand, our observation of softer components 
of quantum radiation, transmitted through 0-5 gm/cm? of aluminium and res- | 


* As long ago as 1931 rough absorption experiments extending to 51 gm/cm? of lead were | 
made by Blau and Karz-Michailova (1931). 
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ponsible for 25% of the total y-ray effect in our counter, is a matter of more interest. 
It is explicitly stated by de Benedetti and Kerner that they could detect no 
radiation softer than the main component of about 0-8 Mev. energy until their 
aluminium absorber thickness was reduced below 0-5 gm/cm?. Likewise the 
lead absorption curve of Siegbahn and Slitis shows little or no trace of any softer 
component. However, these authors used a counter with a gold cathode, and 
if this counter was the same as that later employed by Slatis (1947) in another 
investigation then, from information provided in the later paper, 1t would appear 
that absorption of our soft components in the wall of the counter would have 
amounted to at least 75°% in the arrangement used. It cannot be held, therefore, 


Soft components 


Total y radiation 


55 


0 0-5 10 15 2-0 

Absorber thickness (g/cm?) 
Figure 1. 

x Observations with 50 mg/cm? polythene above the source. 


that serious conflict with the results of Siegbahn and Slatis is established, and 
without further information as to the degree of attention given to the initial 
portion of their absorption curve by de Benedetti and Kerner (their published 
curve shows no plotted point below 1 gm/cm? thickness of lead) it is impossible 
to say whether there is any outstanding disagreement with their observations, 
either. It can, however, be stated quite definitely that the total intensity of the 
soft components observed in our experiments is at least 7 times that of any effect 
which might be attributed to quantum radiations from the radium (D +E) 
impurity in the source—and in any case the abso :bability of our soft components 
is quite different from that of the general quantum radiations from a source 
of radium (D+E+F) in equilibrium. With such a source, under conditions 
similar to those obtaining in the present experiments, the effective mass absorption 
coefficient of the radiation transmitted by 0:5 gm/cm? of aluminium is roughly 
1-1 cm2/gm. in aluminium, or about 25 cm?/gm. in lead. This is entirely different 
from (j/p)py =2:2cm?/gm. which characterizes the main soft component now 
under discussion. 

It should also be stated that there is no possibility that the soft radiations 
arise from disintegration effects, or as a result of inelastic scattering of the 


N 
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w-particles, in the aluminium absorber placed over the source. Observations | 
made with 50mg/cm? of polythene directly above the source (see figure 1) 
showed that these radiations still persisted with the same intensity as_ 
before. 
The comparison measurements with absorbers of lead, gold and tungsten 
were carried out in an attempt to reach more definite conclusions concerning 
the quantum energies of the soft components. On the basis of the lead absorption 
measurements alone it would be possible to assign an energy of about 145 key. 
or 84kev. to the main soft component, depending on whether or not it were 
assumed that K-shell absorption was involved. Figure 2 shows the results. 
of the comparison measurements, corrected for decay of the polonium as already | 
described. It is at once evident from the figure that in general the soft radiation | 
is more strongly absorbed in gold and tungsten than in lead—and slightly more 
so in gold than in tungsten. This result immediately rules out the possibility 
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Figure 2. 


that the main component of soft radiation has quantum energy greater than) 

87-9 kev., the K-absorption energy for lead or less than 80-7 kev., the K-absorption 

energy fot gold. A more detailed (graphical) analysis of the curves of figure 2 
was carried out on the following assumptions : 

(i) the hard component, deduced from figure 1 as responsible for 75°, of 

the total y-ray effect in our first experiment, was responsible for the} 

same fraction of the total effect in the comparison experiment (figure 2)34 

(ii) the mass absorption coefficients of the hard component in gold and) 

tungsten could be calculated from the value obtained in the first experiment: 

for absorption in lead, using standard theoretical expressions. Valuesif| 

of 0-084 cm?/gm. (gold) and 0-079 cm?/gm. (tungsten) were deduced on} 

this assumption. 

The results of this further analysis are summarized in table 1 which include 


the corresponding results from the previous analysis of the lead absorptio 
curve of figure 1 for additional comparison. 
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Table 1 
Experiment Se 56 ss ce = 3 I II II II 
Absorber : a es oe a ee Pb Au WwW 
Main soft component ule (cm?/gm.) A an Se ae oy oe 
effective intensity (°%) ke 16 15 19 19 
More absorbable radiation—effective intensity (°%) .. 9 10 6 6 


It is at once evident from the figures given that it is impossible to provide a 
complete explanation of all these results in terms of two monochromatic soft 
components only; the term ‘“‘main soft component”? which has been used 
hitherto has obviously no unique reference to a single radiation. his latter 
conclusion derives from two considerations, first that the apparent intensity 
of this hypothetical component is greater on the basis of gold (or tungsten) 
absorption than on the basis of lead absorption measurements, and secondly 
that corresponding to (/p)p, =2:2cm?/gm. we should expect (j./p),,=8-4 and 
(/p)w=6-8cm?/gm., instead of 6:1 and 5:7cm?/gm. respectively. Clearly 
there must be present radiation more absorbable than the most prominent soft 
component radiation according to lead absorption measurements, but less 
absorbable than that component in gold and tungsten. Such a radiation must 
obviously have a quantum energy greater than 87:-9kev., the K-absorption 
energy of lead. On the other hand, the presence of a radiation more absorbable 
in all three elements than the most prominent soft component must also be | 
postulated. A quantum energy of less than 35kev. would thus appear to be 
indicated. We may therefore accept a very rough three-component inter- 
pretation of the soft radiation as follows (table 2) 


Table 2 

Quantum energy (kev.) .. Sho oe 8444 105+15 

Effective intensity (°% a an 6 15 + 
Further discussion of this interpretation is reserved for a later stage, here it 
should merely be re-emphasized that the figures given under heading “‘ effective 
intensity”’ in table 2 represent percentages of the total y-ray counting rate in 
our experimental arrangement, not absolute, or even relative, values of true 
quantum intensities per disintegration of polonium. 


§3. DELAY EXPERIMENTS 

These experiments were designed in the first instance to investigate the 
possibility that the residual nucleus 2°6Pb, which results from the «-disintegration 
of polonium, is occasionally produced in a metastable state of lifetime of the 
order of minutes; they were later extended to search for possible y-ray lifetimes 
as short as a few milliseconds. The first experiment originated in the suspicion 
(Feather 1946) that the unidentified activity of 1:6m. half-value period, which 
Waldman and Collins (1940) have reported as produced by the action of high- 
energy x-rays on lead, might result naturally in a rare mode of the disintegration 
210Po->206Ph, As already stated, this experiment gave a negative result. It 
was carried out by exposing a brass disc collector of 1-4cm. diameter at a distance 
of 1cm. and at a potential of —700 v. with respect to a “‘clean”’ source of about 
12mc. of polonium deposited on a platinum foil of 4cm? area. Exposures 
lasted for a standard time of 12 minutes. Any activity collected by the disc 
was then observed over a period of 10 minutes immediately following the end 
of the exposure, using a thin-windowed counter and an absorber of 15 mg/cm* 
paper to cut out the effect of the #-radiation of the polonium which was collected 
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by aggregate recoil. As a result of 14 exposures it was concluded that the — 


collector carried no activity giving an initial counting rate of more than 1 count 


per minute under our experimental conditions. The a-activity of the polonium | 


collected by aggregate recoil was estimated, by removing the paper absorber, 
as 1-4min~! per minute exposure in the same geometrical arrangement, this 


figure corresponding roughly to one polonium atom collected per 3000 «-disin- 


tegrations of the source. 


The next experiment, on possible short-period delays, was carried out in — 


two different arrangements exemplifying the same basic principle. This prin- 


ciple was to examine the activity of a collector permanently exposed to a polonium | 


source, the source and collector being so disposed in the neighbourhood of the 


end-windowed counter that the counter was considerably more sensitive to | 
y-radiation emitted by atoms disintegrating on the collector than to y-radiation — 


from the primary source.* Transference of recoiling atoms from the source 
to the collector in the applied field was calculated to take place in an average 
time of 2-4 x 10-sec. in the first arrangement and 5 x 10-*sec. in the second. 
In spite of the longer transference-time the second arrangement had the great 
advantage over the first that the collector subtended a much greater solid angle 
at the counter window than before: in the second arrangement, in fact, the 


collector was separated from the window only by a0-5 cm. air gap and analuminium jf} 


absorber 30 mg/cm? thick. In each arrangement the counting rate with source 
and collector in position was determined both when the source and collector 
were at the same potential and also when the collecting field was applied. A 
“genuine” increase in the counting rate on the application of the field would 
in this experiment provide clear evidence for the production of an activity in 
the recoil nucleus 2°*Pb of any lifetime between a few milliseconds and 1 minute. 
In fact, the mean increase in three sets of observations using the first arrangement, 
was 0-80 + 0:53 min“ and, in five sets of observations using the second arrange- 
ment, 0:-40+0-56min+. The y-ray counting rate determined with the primary 
source in the collector position in arrangement (i) was 20-5 +2:-4min4, and it 
was estimated that in arrangement (1i) it would have been at least 60 min-?. (As. 
already implied a direct determination of the y-ray effect of the polonium source 
through only 30mg/cm? aluminium could not easily be made on account of the 
radium (D+E) impurity). The efficiency of collection of 2°*Pb recoil atoms. 


was estimated—in the only way available, in terms of the efficiency of aggregate | 


recoil collection—as about 3 of the efficiency in the standard exposure arrangement 


of the long-period delay experiment (see above) in arrangement (i), and as about. 


4 of this standard efficiency in arrangement (ii). The latter result was all the 


more gratifying since the average distance travelled by a recoil atom in the ] 
field—and that not in a straight line—was about 10cm. in this arrangement, } 


as against 1 cm. under “standard” conditions. 


§4. DISCUSSION 


In discussing the experimental results described above we may deal first, | 


quite briefly, with the results of the delay experiments. Always assuming 


that the recoil efficiency of the primary polonium source was high—which the ]} 
appearance of the activated platinum foil and our observations on aggregate | 


recoil might be held to justify—it is obvious that these results show that not more 


* Essentially the same principle was used by Ellis (1932) and also by Wright (1932) in work 
on ThC and RaC’, respectively. 


———————— 
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than 10% of the y-ray effect in our counter arises from activities of lifetime 
greater than 5 x 10-%sec. and not more than about 3° of the total effect from 
activities of lifetime greater than 10-'sec. It is significant to remark that this 
second statement includes y-ray activities which because of a high degree of internal 
conversion might show up only through the emission of L x-rays (Curie and 
Joliot 1931), since the aluminium absorber in the second arrangement was 
sufficiently thin to transmit a considerable fraction of such fluorescent radiation. 

Turning to the results of the absorption experiments one important conclusion 
can be stated at the outset. Because of the close similarity of the absorption 
curves obtained with gold and tungsten absorbers it must be concluded that no 
significant component of the “soft” radiation from polonium has quantum 
energy intermediate between the K-absorption energies of tungsten and gold, 
i.e. in the energy range between 69-4 and 80:7kev. Now the Kz, and Ka, 
radiations of elements 81 to 84 inclusive lie within this range, thus we have the 
important result that the soft radiations which we have detected do not include 
the K-radiations of lead. That these radiations must be emitted following the 
internal conversion of the 0-77 Mev. polonium y-ray is undeniable, but they are 
evidently below the limit of detection with our arrangement : present considera- 
tions, in fact, support previous estimates of the absolute intensity of the photo- 
electron lines (Feather 1946) and show that the internal conversion coefficient 
for the 0-77 Mev. radiation cannot be abnormally large (not greater than about 
10~*). Since they are not K x-rays obviously the soft radiations of table 2 must 
be predominantly of nuclear origin. Before we can make much progress towards 
an understanding of their mode of emission, however, we must make some 
attempt to estimate their intensity. ‘This can only be done in the present instance 
by a rough comparison with the harder radiation. 

Let us consider first the intensity of the main hard component. ‘The “ best”’ 
estimate in this case would appear to be (1-0 + 0-4) x 10~° quanta per disintegration 
(Feather 1946)—and no difficulty would ensue if this were the only y-radiation 
emitted. For, if the same change of spin were involved in an «-disintegration 
leading to the 0-77 Mev. excited state as in that leading to the ground state of the 
nucleus 2°$Pb, we should expect an excitation of the 0-77Mev. state of 
(1:8+0-5)x10-* per disintegration (Feather 1948). The agreement between 
this figure and the observed y-ray intensity is remarkably close, and in the 
absence of the soft radiation we might well be satisfied that the rare nuclear 
processes were completely understood. We have now to see whether it is possible 
to explain at the same time the emission of the softer y-radiations in terms of the 
excitation of other states having spin values consistent with our observation 
that the soft y-radiation is not appreciably ‘“‘delayed”. In this connection 
we should note the estimate of Curie and Joliot (1931) that the intensities of the 
L- and M-radiations of polonium (or lead) which they observed were about 
4x 10-4 and 1-5 x 10-3 quanta per disintegration, respectively. 

Now it is certain from our results that the soft radiations are not considerably 
more intense (in quanta emitted per disintegration) than the hard component. 
Using the results of Bradt and collaborators (1946) in making our comparison 
we should in fact conclude that our 84+4kev. soft component has roughly the 
same intensity as the 0-77 Mev. (hard) radiation.* Our other soft components 


* Bradt et al. found that an aluminium counter is roughly five times as sensitive to y-rays of 
(-8 Mey. energy as to 0-1 Mev. radiation. 
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would then appear as less intense than this. As concerns the main soft component 
two possibilities must now be considered: either the major portion of the 
L-radiations of Curie and Joliot follows the internal conversion of this soft 
y-radiation (of 84+4kev. energy) or it does not. On the first assumption the 
84 kev. nuclear transition has an intensity of roughly 4 x 10~ per disintegration, 
the radiation has an internal conversion coefficient of about 40, and because of 
intensity considerations it cannot be coupled in any way with the emission of 
the 0:77 Mev. radiation. Separate excitation of a nuclear state of 84kev. energy 
has thus to be assumed, and a difference of spin between excited and ground 
states of some 5 units is necessary to explain the intensity. The assumption 
has therefore to be ruled out if the radiation is not appreciably delayed. 

The second assumption avoids this particular difficulty by allowing coupling 
of the hard and soft radiations—and so permitting a single high energy state 
excitation. We might, in fact, on this view, imagine that the 0-77 Mev. state 
was fed chiefly by y-ray transitions from a state of 84kev. higher energy (and 
possibly to a smaller extent from other states in the same energy range) excited 
directly by «-disintegration. The difference of 84kev. in the energy available 
for x-disintegration would be expected to introduce an adverse factor of no 
more than 3-5 in the partial disintegration probability (which would not seriously 
disturb the numerical agreement already noted), and if it were assumed, for 
example, that the ground state of 2°*Pb and the excited state of 0-85 kev. each 
had zero spin, and the 0-77 Mev. state had spin of 2 units, no major interpretative 
difficulty need remain—except the explanation of the intensity of the L- and 
M-radiations observed by Curie and Joliot and the whole problem of the fine- 
structure spectrum of Chang (1946). Concerning the first difficulty it might 
be possible to accept the suggestion, frequently. put forward in the past, that 
the fluorescent radiations arise in the (infrequent) internal ionization of the 
radioactive atom by the escaping «-particle, having no connection, therefore, 


either with the appearance of discrete lines in the «-particle energy spectrum, ~ |] 


or with nuclear y-ray transitions, but the problem posed by Chang’s results 
would seem as far from solution as ever. Amongst the rare effects associated 
with the polonium disintegration only the fluorescent x-ray emission itself 
appears to have an intensity of the same order of magnitude as the fine-structure 
lines in Chang’s spectrum. 
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ABSTRACT. ‘The experimental data on the daily variation in cosmic-ray intensity are 
first surveyed. It appears to be established by Rau’s experiments that the energetic 
meson component at sea level shows a marked semi-diurnal variation which is in phase with 
the semi-diurnal variation shown by the barometric pressure. Several authors have sug- 
gested that this phenomenon is explicable in terms of the Pekeris theory of atmospheric 
oscillations. The implications of this explanation on the process of meson formation are 
here examined quantitatively. It is concluded that the explanation is only possible -if 
mesons arise mainly at about 60 or 70 km. above sea level, which is highly unlikely since 
the corresponding cross-section for meson production would be much larger than appears 
plausible. 

In addition the experimental data on the diurnal variation are briefly discussed. This 
variation seems to be more complicated in character than is usually assumed and not 
attributable to factors such as a heliomagnetic field and fluctuations in the geomagnetic field, 
which are often held responsible. 


ol IN LRODUC LION 

F the experimentally observed solar daily variation in cosmic-ray intensity is 

analysed into Fourier components, it is usually found that the only important 

components are the diurnal (24-hour period) and the semi-diurnal (12-hour 
period). It is with the semi-diurnal component that we are chiefly concerned 
here. 

Certain observers (Regener and Rau 1939, Ehmert 1939, Rau 1940, Sarabhai 
1945) have found a daily variation in meson intensity which is mainly semi-diurnal 
and shows a strong positive correlation with that in atmospheric pressure, and have 
pointed out that the theory of atmospheric oscillations outlined by Pekeris (1937) 
and described in §3 below provides at least a qualitative explanation of this 
phenomenon, which is difficult to understand on any other basis. If this explan- 
ation is valid, then the experimental data on the amplitude of the cosmic- 
ray effect, together with the quantitative predictions of the Pekeris theory, 
provide a means of determining the approximate height at which mesons are mainly 
produced and hence the approximate value of the average cross-section for meson 
production. At the present stage of theory and experimental data relating to 
meson creation, even a very rough estimate of the value of this cross-section would 
be useful. 

It was clear from the general features of the oscillation predicted by Pekeris 
that this method of explaining the semi-diurnal variation in meson intensity would 
certainly lead to rather a large value for the cross-section, but whether this value 
would be too large to be plausible could not be decided without quantitative 
investigation. ‘The main purpose of this paper is to describe such investigations 
and their result (§3). One of the subsidiary purposes is to make clear the impli- 
cations of the Pekeris theory, for the results of this theory have been incorrectly 
interpreted by certain writers on the semi-diurnal variation in cosmic-ray intensity. 
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Since there is much apparently conflicting evidence on the daily variation of | 
cosmic radiation, it is essential to survey all the experimental data before proceeding 
to consider theories to explain any part of these data. We shall now make such a | 
survey and then consider what facts may be looked upon as experimentally 
established. 


§2. EXPERIMENTAL DATA 

It is important to classify the various experimental results according to the | 
nature of the intensity to which they refer, for the daily variation may be expected _ 
to depend upon the constitution of the radiation observed. For example the 
mesonic and electronic components may not be affected in the same way by 
short-period fluctuations in meteorological conditions, since these fluctuations may 
be entirely different in character in the upper atmosphere where the mesons mainly | 
arise and in the lower atmosphere where most of the electrons recorded at sea-level 
arise. 

In addition it is important to make clear whether the observer in presenting his © | 
results has tried to allow for changes in meteorological factors such as temperature | 
and pressure, for it often happens that the so-called ‘‘ corrections”’ applied are so | 
large that they completely change the character of the daily variation. For 
example, the ‘external temperature” correction, applied by many workers in an 
attempt to eliminate the effect of fluctuations in the temperature of the air above the 
recording apparatus, may be sufficient to shift the phase of the diurnal variation by 
eight hours; it is therefore important to consider the justification for this correc- 
tion. ‘The correction is usually made before harmonic analysis is carried out, by 
adjusting the hourly recordings by means of a temperature coefficient determined 
from the correlation of long period (usually average monthly) variations of intensity. 
with temperature. Now the variations in monthly mean intensity at sea level are 
probably caused mainly by changes in the average height of origin of mesons; 
this height depends upon the temperature distribution throughout the atmosphere 
below it and also on the pressure at sea level. Consequently, even in these long 
period variations, there is not always a simple dependence of the intensity on the 
temperature at the ground; even if there were it would have little bearing on the 
relation between hour-to-hour changes in intensity and ground temperature. 
Test made by Forbush at Cheltenham (U.S.A.) suggest that at this location at 
least, no external temperature correction ought to be made. - Consequently the |} 
value of results to which this correction has been applied seems very doubtful. 

The practice, often adopted, of correcting hourly readings for fluctuations in 
barometric pressure by using a barometric coefficient based on comparison of the 
mean daily intensity and barometric pressure, may not be entirely justified either; |] 
for although the barometric effect is certainly due partly to mass absorption, it may | 
also depend upon factors such as movement of meson-producing layers, which are | 
differently related to short and long period pressure changes. A third correction | 
is sometimes applied for fluctuations in the internal temperature of the apparatus, 


but in the best designed experiments this is made unnecessary by thermostatic ]} 
control. 


2.1 Collection of results 


Results are collected in table 1. ‘The notation used therein is as follows :— 
A= geomagnetic latitude of the station. 
h=altitude of station in metres above sea level. 


i 
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é = angle of cone formed by observed radiation; since information on 
this is not very precise it will be classified simply as wide or 
narrow. 

@=inclination of axis of observed cone to the vertical. 

D(4,c); S(4,c) = diurnal (semi-diurnal) variation of amplitude 5°/ and maximum 
atc hours. In cases where the harmonic dial method of analysis 
has been used the probable error in b is also given. 

+S; —S imply that the maximum of the variation is displaced by less 
than one hour from the maximum (minimum) of the semi- 
diurnal pressure variation. 
I.C.; G.C. indicate use of ionization chamber; Geiger counter apparatus. 
B.; E.T.; 1.1. indicate correction for barometer effect ; external temperature 
effect; internal temperature effect. 


2.2 Discussion of the experimental evidence on the S variation 


In view of the smallness of the daily variation, observations over a very long 
period are required before one can be certain that the variation obtained is of 
physical significance. For this reason the reliability of many of the results quoted 
above is somewhat doubtful. In the case of the S variation the only results of 
convincing statistical significance are those of Duperier (1945), based on three 
years’ observations and those of Rau (1939), although the results found by Hoerlin 
under 10 cm. Fe do carry some conviction on account of the large amplitude of the 
S effect recorded and the closeness with which it follows the pressure variation. 

Rau used two independent ionization chambers of different size suspended 
under 40 metres of water in a narrow vertical fissure in the rock forming the floor 
of Lake Constance. The dimensions of this fissure were such that the intensity 
observed was mainly vertical. ‘lhe data from the two chambers showed a similar 
variation. The results given above were obtained with the larger chamber. The 
first result S (0-18 + 0-04, 11" 30) was found by a harmonic dial analysis of the 
data recorded on 192 days, while the second result 5 (0-14, 10") was obtained by 
taking the average of the hourly readings over a period of 15 months, and then 
harmonically analysing the resulting daily variation curve. ‘The two results when 
considered together, suggest that the amplitude lies between 0-14 and 0-18%, and 
we shall in numerical work take 0-16°% as the amplitude established by Rau. 

The interpretation of Rau’s results is much less ambiguous than that of 
Duperier’s for the following reasons : (i) Only mesons are measured instead of the 
total intensity made up of mesons and soft component. (ii) The angle ¢ is small 
and consequently the radiation can be regarded as entirely vertical to a good 
approximation. (iii) Only particles of initial energy >10!° ev. can penetrate the 
absorber above Rau’s apparatus (50m. water including the atmosphere), which 
means that the effect of magnetic variations has been eliminated. ‘This 1s con- 
firmed by the absence of any change in intensity during magnetically disturbed 
days in the case of Rau’s experiments. 

In view of the difference in nature of the intensity measured in the two cases, it 
is not surprising that the variations found by Rau and Duperier are of very different 
character. 

At the equator there are no completely convincing results. ‘Thompson’s are 
uncertain on account of their short duration (only about 6 day’s readings for a 
zone 23° in latitude) and the fact that a relatively large correction is made necessary 


‘Table 1 (i) Variation of radiation 


Observer r h Apparatus fo) 6 
Barnothy and Forro 47° N. 124 CHO, wide 0) 
(1939) ; 
Duperier (1945) 54° N. 0 Git wide 0 
Alfvén and Malmfors 58° N. 0 G.C; narrow S0NE 
(1943) 30 Se 
60° N 
60S: 
Hoerlin 2NE 5500 to sie wide 0) 
reported by 6100 
Rau (1940) 2500 to Ce wide 0 
5500 
Kolhorster (1941) 52:4°N. 0 GG: 45°N. 
ATE Sy 
45° E. 
45° W. 
Sarabhai (1945) SeaNt 900 EGG wide | 0) 
(ii) Variation of total 
Steinmaurer (1935) 47°N. | 590 LC. wide 0 
Hess and Graziadei 47°N. 2300 1G: wide 0 
(1936) 
Doan (1936) (2) S2aeNe 150 1eKor 0 
Forbush (1937) (3) 50° N. 72 1Ge wide 0 
Schonland, Delatizky 8221S 122 mes wide 0 
and Gaskell (1937) (4) 
Thompson (1938) (5) 40-55° N. 0 1G. wide 0 
25-40° N. 
10-25° N. 
10° S.-10° N. 
10-25° S. 
25-40° S. 
Sarabhai (1945) SeeNE 900 ice narrow 0) 


(111) Variation of 
Barnothy and Forro 47° N. 124 ar) 0 
(1936) (6) # 50° 


Rau (1939) (7) i 50° N. 360 LG. 


narrow 0) 


(iv) Variation of | 
Barnothy and Forro ; 
(1939) iS | | 0 


SS SS a 
(1) aoe eet Appreciable variation with 0; time of maximum varies 
(2) Very small barometric variation on 10 days selected. Results from 7 

meters running simultaneously. 
(3) D variation well established. 
(4) Results made uncertain by large temperature fluctuations of apparatus itself. 
(5) Recording made during 12 voyages. On average 6 days readings for each zone 23° in | 
latitude. Uncertainty introduced by necessity of correcting for latitude effect. 


v7 aeNG 124 


identical recording 


under little or no OSE So eee eae al ee 


Pot sete ler Sa | 


Absorber ee Corrections Results 
Begefie: applied 
S 
0 2 months 
0 3 years None (O25 2ai7e sey —(0-18, 0312) 
B. (O: 2. 17h 20) | 
0) 24 years B. (0-:12+0-03, 10h) 
(0-08 +0: 03, 13h) 
(0-16 +0- 02, 13h) 
(0-17-£0- 02, 19h) 
0 10 days Bayelels —(0-06, 04h 30) 
10cm. Fe}] 8 weeks Beaalers +(0-78, 10 47) 
0) 8 months None 
0 15 days None (0:3, 14h 30) 


meson component 


10cm. Pb] 9 months B. (0-25, 8h) 
By, 1oyIe (O;25,012)0) 
12 cm. Pb 3 years Bypbers (0:2, 104 40) (0-5, 00h) : 
12tcmea eb 10 days B. (0:19-+0-04, 09h) 
12cm. Pb | 273 days B. (0:17+0-02, 11h) (0-04+0-01, 004 30) 
12 cm. Pb 3 years B. (0:2, 16h) 
Be, 1eTe (0:2, 6h) 
Ion With, TAI (0:2, 14h) 
12 cm. Pb | 11 months B. (0:33, 14h 20) 
in all (0-22, 14h 30) 
(0:26, 148 20) 
(0:23, 13h 50) 
(0-21, 14h 30) 
(0-17, 13h 10) 
Scm. Pb* 12 days None 0) +(3-0, 10h 30) 
energetic meson intensity 
36 cm. Pb | 4 months B. 
36 cm. Pb | 3 months B. 
40 m. H,O} 15 months None 0 (a) +(0:18+0-04, 11h 30) 


(b) +(0-14, 102) 


electronic shower intensity 


? 12 months B. (0-41, 15h) 


(6) Curves suggest reversal of phase of S variation as 6 changes from 0 to 50°, but physical 
significance of these results is not convincing. 
(7) (a) Harmonic dial method used on 192 days observations. 
(b) By harmonic analysis of mean of 15 months observations. A second smaller, indepen- 


dent I.c. was used simultaneously and showed similar variation. 
* In this case the method of shower anticoincidences was actually used for excluding the soft 
component, but the arrangement is essentially equivalent to a counter telescope with about 


8 cm. Pb akssorber. 
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by the latitude change in a 24-hour period. Sarabhai’s experiments are unfortu- 
nately also too brief to give convincing results; they suggest the presence of a 
large +S variation. Hoerlin’s results, obtained during the 1932 Cordillera 
Blanca Expedition and re-examined by Rau (1940), carry the greatest conviction. 
They suggest very strongly that the variation under an absorber of 10cm. Fe is 
semi-diurnal, with large amplitude and phase strikingly close to that of the semi- 
diurnal pressure variation (meson maximum at 10"47, pressure maximum at 
10°50). It seems likely, therefore, that at the equator the meson component 
shows a large + S variation but further experiments are necessary to establish this. 

However this may be, Rau’s conclusion that at latitude 50° the energetic 
mesons are characterized by a +S variation, appears to have good foundation 
and it seems worth while to examine quantitatively whether such a variation can 
be explained in terms of the Pekeris theory of atmospheric oscillations, particularly 
as this explanation of the effect has been given in several independent papers. 


§3. EXPLANATION OF THE -35 VARIATION 

A predominant semi-diurnal variation is not likely to result either from 
temperature or from magnetic variations, as both would produce mainly diurnal 
rather than semi-diurnal effects. Apart from this consideration, Rau’s results 
show clearly that magnetic fluctuations are not responsible, as the particles he 
observes have energies too high to be influenced by them *. 

It has already been pointed out independently by Rau and Ehmert (1940) and 
by Sarabhai (1945) that a + S variation would result from atmospheric oscillations 
of the kind predicted by Pekeris, provided that the bulk of the mesons arise above 
30km., at which height the Pekeris oscillation has a nodal surface. Before 
working out the implications of this interpretation we will briefly indicate the 
nature of these atmospheric oscillations. 


3.1 The theory of atmospheric oscillations 


The diurnal component of the daily variation in atmospheric temperature is 
considerably larger than the semi-diurnal component, as one would expect. In 
spite of this, however, in the variation of the barometric pressure at sea-level, the 
semi-diurnal component is dominant, being about 14 times as large as the 
diurnal component. An explanation of the large semi-diurnal pressure variation 
has long been sought. Kelvin suggested that the atmosphere has a free oscillation 
of period about 12 hours, so that the semi-diurnal tide excited by the sun is magni- 
fied by resonance. Since then many attempts have been made to determine the 
period of free oscillation of the atmosphere. 

G. I. Taylor (1929, 1932, 1936) showed that this period can be calculated from 
the speed of propagation of waves of explosion such as that caused by the Krakatau 
eruption in 1883, and that for this eruption the speed with which the wave was 
propagated corresponds to a free period of about 10} hours; a similar speed was 
found by Whipple (1934) for the wave of the Great Siberian Meteor. A periodof | 
105 hours is, however, fatal to the resonance theory of atmospheric oscillation. 
In 1937 Pekeris extended the theory which had been outlined by Taylor in 1936, 


* In the case of Rau’s experiments, tides in Lake Constance provide a possible cause for an 
S effect, since any change in surface level of the lake alters the thickness of the water absorber above 
the apparatus. ‘The change in intensity produced by such tides appears to be negligible however , 
as their magnitude is very small indeed—the change in level from high to low water being of the 
order of 0-1 inch. 
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and was able to prove that, when the temperature distribution with height sug- 
gested by experimental data is assumed, the atmosphere can have a free oscillation 
of 12-hour period in addition to one of 104-hour period. 

As early as 1924 Chapman discussed the semi-diurnal oscillation excited partly 
tidally and partly thermally by the sun, and concluded that if the atmosphere has a 
natural period of 12 hours, then this oscillation will be magnified about a hundred- 
fold by resonance. Recently a more detailed investigation of the magnitude of the 
semi-diurnal pressure oscillation at different heights has been carried out by 
K. Weekes and M. V. Wilkes, who have extended the basic theory outlined by 
Pekeris and have used the Cambridge Differential Analyser to obtain numerical 
solutions of the equations involved. An account of this work together with 
detailed numerical results is to be published shortly. Meanwhile, Weekes and 
Wilkes have very kindly made their results available to the present writers. 

The mode of oscillation of 12-hour period is characterized by a horizontal nodal 
surface at a height of about 30km., the atmosphere above and below this level 
swinging in opposite directions. ‘The changes in pressure and the horizontal 
velocity of winds caused by the oscillations are reversed in phase in crossing this 
level; the vertical velocity of winds has the same phase at all levels. * The 
amplitude of the pressure oscillation is small and almost uniform up to nearly 
30km.; at high altitudes it increases rapidly with increasing rarefaction of the 
atmosphere. ‘The feature of particular significance in considering the effect of 
such pressure oscillations upon meson intensity at sea level is the variation AA in 
the height of an isobaric surface situated at mean height h. ‘Table 3 shows the 
values of AA as a function of h calculated by Weekes and Wilkes for an atmosphere 
having the temperature distribution given in table 2, a distribution based as far as 
possible upon experimental data. ‘The values given for AA correspond to an 
oscillation of amplitude 1 mm. Hg at sea level at the equator, which is very close 
to that experimentally observed. ‘The amplitude of the oscillation decreases with 
increasing latitude by a factor approximately cos*A; Af is given for A=0 and 50°. 


Table 2. The variation of temperature with height adopted in 


the calculation of table 3 
‘Temperature (° K.) 


Height (km.) Evidence 


O- 10 288-220 ; linear fall. Up to about 25 km. well established 
10— 30 220 by radio-sonde data. 
30— 55 220-350 ; linear rise. Inferred from the abnormal propa- 
gation of sound (Whipple 1935). 
55— 60 350 | : : 
Pre Suggested by radio exploration 
Oe USI 9 eee (Budden, Rarclificand Wilke 1039) 
77-100 190 J 
Table 3. The semi-diurnal oscillation of isobaric levels at various altitudes 
h (km.) 0 10 20 30 40 50 
A=0 +0-010 +0-012 +0-012 —0-002 —0-086 —0-272 
See O 2250? 0-003 +0-003.2 40-003 . =0-001- —0-023" . —0-072 
h (km.) 60 70 80 90 100 
RO —0:-49 —0:51 —0-60 —1-05 —2:22 
Ah (km.) A=50° —0-129 —0-135 —0-160 —0-28 —0:59 


* It is interesting to note that this type of oscillation is compatible with that suggested by the 


“ dynamo ”’ theory developed by Chapman (1919) to account for the semi-diurnal variation of the 
earth’s magnetic field. Chapman finds that the pressure oscillations in the upper conducting layer 
where the dynamo effect is produced are nearly 180° out of phase with those observed at the ground. 
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These values do not agree with those calculated by Mailvaganam (1946). 
From the following considerations it seems that Mailvaganam’s results cannot be | 
valid. It can be seen from the results of Pekeris’ own analysis that near resonance | 
the horizontal velocity is proportional to the pressure oscillation at any height, | 
and consequently that the curve given by Pekeris (1937, figure 1) for the horizontal _ 
velocity also represents the pressure oscillation. Now the pressure oscillation | 
calculated by Mailvaganam does not agree even in character with that given by 
Pekeris. ‘The most obvious discrepancy is that Mailvaganam’s results do not 
show the node at 30km. which characterizes all Pekeris’ results. One of the 
causes of the discrepancy seems to be that Mailvaganam has used incorrect numer- 
ical values for constants (C and D in his notation). We are indebted to Weekes | 
and Wilkes for investigations on this point. 

It is apparent from the nature of the oscillation described above that mesons 
created above 30 km. will show a +S variation at sea level, for if the isobaric levels 
in the lower atmosphere move down, corresponding to a fall in sea level barometric 
pressure, isobaric levels above 30km. move up; mesons are therefore formed at 
higher levels and the meson intensity at sea level is decreased as the probability 
of decay is enhanced. Also the +S variation arising in this way should be much 
more pronounced at the equator than at latitude 50°, as the pressure oscillation 
increases by about a factor 4 in going from A=50° toA=0°. The experimental 
data at present available at the equator are compatible with the occurrence of such _ 
an increase in amplitude of the S variation, although not sufficient to establish it. 

We shall now calculate the average height at which mesons would have to arise 
to produce a +S variation of the magnitude observed by Rau. 


3.2 Calculation of the average height of creation of energetic mesons 


3.2.1 Treatment assuming all mesons arise at one height 

In this calculation we shall assume that all mesons arise at height X above 
sea-level, and that at this height they have a differential energy spectrum propor- 
tionalto E-”. It seems probable that mesons at a high level will have a spectrum 
of this form; the assumption is unlikely to introduce a large error, as the results do 
not appear to depend critically upon the energy spectrum. Usually n will be taken 
equal to 3, but one case with m = 2 has been worked out to test the dependence of the 
result upon n. 

The probability P(E, X) that a meson of energy E will travel distance X without | 


decaying is 
7X 
pe dx 
exp( Nites ms) 


where p, 7 are the mass and the lifetime at rest of the meson. The only mesons 
we shall be concerned with in interpreting Rau’s results have E101 ev. and for |} 
these the energy loss in penetrating the atmosphere may be neglected so that | 
P(E, X) = exp (—«X/E) where « is constant (=e/7). | 

The total number W(X) of mesons recorded at sea level 


Fie | - E-"exp(—aX/E)dE, 
Ey 


es . | 
where Ey is the lowest energy recorded by the apparatus and K is a constant. | 
] 


) 
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‘The change 6M in M produced by a change 5X in height of creation of the mesons is 


| Kx Ee exp —2X/B) dE | Se ee (1) 
Ey 
so that 

a “dX | pe exP ( —2X/E) dE 

ue | E exp(—aX/E) dE 

Ey 
=a0X . F(Ey,X) say. 

or 


SX=8M/(aM F(EyX)). see (2) 


For 7 an integer the integrands involved in F(E£,, X) can be integrated by parts 
and the function F can be readily tabulated against X for a given E). 

If 5M/M is put equal to the amplitude of the semi-diurnal variation experi- 
mentally observed, then 6X gives the movement of the meson-producing layer 
required to account for the observed variation at sea level. On the average, 
mesons will be produced when the primary particles have penetrated a certain mass 
of air, 1.e. at a certain isobaric level in the atmosphere. The mean height of this 
level may be found by comparing the 5X determined from equation (2) with Ah 
given in table 3. 

As Rau’s results represent the most accurate determination, this process was 
carried out using his value for 6M/M and Ey. For H,=10ev. and n=3, 
F(E), X) is a very slowly varying function of X and 5X is about 160 metres for all 
values of X between 30 and 100km.* Using table 3 this gives about 80 km. as the 
mean height of meson production. 

Ehmert (1940), in a very rough calculation using a mean decay coefficient, 
estimated that a movement of 270 m. in the meson-producing layer was necessary 
to give a variation of 0-18%, in meson intensity at sea level. The discrepancy 
between this value and 6X found above seems to be largely due to the fact that 
Ehmert, following Euler and Heisenberg (1938), has taken different values for 
and 7 (160 instead of 180 x mass of electron for , and 2-7 x 10~® instead of 
2:1 x10-®sec. for 7). If Ehmert’s value for 6X is used however, the height of 
meson formation is only altered by about 8km. . 

In the above calculation it has been assumed that mesons all arise at the same 
height, but in actual fact meson production will be spread over a considerable range 
inheight. The effect of this spread will be to lower the height at which maximum 
production occurs, for the oscillation of the isobaric levels given in table 3 increases 
rapidly with height above 65km., and mesons produced higher than the average 
have a much greater effect on the variation at sea level than those produced lower. 
Although the result of the simple calculation described above is not seriously 
modified by this consideration, a more accurate calculation can be made as shown 


in the next section. 


* All mesons are assumed to have the same mass (200 X mass of electron) and the same lifetime 
at rest (2°15 x 10 sec.). 
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3.2.2 Treatment taking account of the distribution of meson production with height 


Suppose that 
P(h) =number of primary particles per unit area at height h above sea level. 
a = cross section (assumed constant) for production of mesons by these primary | 
particles. 
p(h) =number of nuclei per unit volume at height h. 
f(h)=number of mesons produced per unit length. 


Then 
f(h) =dP/dh=c P(h) p(h). wae Sa 
For an isothermal atmosphere, 
p(h) = poe et, ee (4) 
where H is a constant often called the scale height; in this case 
éPidh=pycP De ae eae (5) 
By integrating (5), 
InP=—Ae-Y where A=Hpo; y=h/H.  _...... (6) 
Hence 
dPjdy=Ae * &xp (— Ae) =2e4, |) eee (7) 
where . 
z=Ae¥Y and Ingzg=lInd—y=-=—y-+constant. ..<..- (8) 


It follows that a plot of ze~* against In z will represent the rate of meson 
production, f(y), as a function of y; this is shown in figure 1. The position of 
sea level (y=0) in this diagram is determined by the value of In A which depends 
only upon the cross-section—a change in cross-section displaces the curve hori- 
zontally but does not alter its shape. The fact that the same curve applies for all |]} 
values of o greatly simplifies the numerical work involved. 

The curve shown 
in figure 1 strictly 0-4 
applies only for an iso- 
thermal atmosphere ee 
(H constant). The f) 
calculation of this y 
curve is considerably O2 
more complicated 
when H is a varying 0-1 
function of height as 
in the true atmo- 
sphere; it has been -8 EF “4 -2 0 
carried out for the In(Aby 
true atmosphere for Figure 1. Rate of meson production per unit length as a function | 
one particular value of height. 
of co. As the curve 
obtained in this case was very similar in shape to that for H constant and only 
displaced slightly in height, it has been assumed that the variation in H may b 
neglected for present purposes. — 


2 4 6 
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Figure 1 gives the rate of meson production for a production process with a 
particular cross-section «, In actual fact o will depend upon the energy of the 
primary particle and also upon the energy of the meson it produces, and the meson 
production actually occurring in the atmosphere will be the result of superposing 
curves corresponding to a whole range of values of c. However, as there is no 
reliable theoretical or experimental knowledge of the dependence of ¢ on energy, 
it seems reasonable to assume for the entire produc tion an average cross-section 
independent of energy. 

Suppose the mesons have a differential energy spectrum E~” at their point of 
origin (7 will be taken equal to 3 in numerical work). Then M, the total number of 
mesons reaching sea level 


nD 


= | a) i i exp aa hgh ae. ae ay (9) 


We require to find the change tn M resulting from a movement Ah of the isobar, 
whose mean position is at height # above sea level. We shall assume that the rate 
of production of mesons remains the same at a particular isobaric level throughout 
the semi-diurnal atmosphere oscillation. This is justified provided the density at 
a given isobaric level does not vary significantly, which is the case below 90 km. at 
least. With this assumption 


Sie: | fla)Ah dh i E-texp(—ah/E)dE, «10... (10) 
h=0 Ey 
‘so that 
elas | _ fh)Ah F,(h) ah 
= nad RG Mir eck Vt (11). 


OO 


where Fy=| E-*exp(—ah/E)dE, F,=| E-*exp(—ah/E)dE. 
Ey Y By 
ae (12) 


F, and F, can be found by integration by parts and are readily tabulated as 
functions of h. The integrals occurring in (11) are then evaluated by numerical 
integration. This evaluation has been performed for three different values of a, 
using the values of Ah provided by Weekes and Wilkes at latitude 50° and taking 
E, = 10" ev., in accordance with Rau’s measurements. ‘The results are given in 


table 4. 


Table 4. The semi-diurnal variation in meson intensity at sea 
level for three different values of o 


In A ae ae aes ef 6 8 9 
APs a ee a 5S e&~400 e8~3000 e®~8100 
Ht. of max. meson production (km.) 48 64 72 
6M/M ie ae ts Ne 0-0008 0-0014 0:0017 


Rau’s measurements give 5M/M=0-0016, and so to explain these measure- 
ments by means of the Pekeris theory the mesons have to be created mainly at a 
height of about 70 km.* 


* Before the calculation which resulted in table 4 was carried out it was thought that the height 
.of meson formation might be reduced to as low as 50 or 60 km. by taking account of the spread in 
~ height of meson formation, and ina brief statement on this work at a Physical Society Conference in 
July 1946, about 50 km. was suggested ‘as the mean height of meson formation required to explain 
Rau’s results (Sarabhai and Nicolson 1947). 
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3.3 The high cross-section 


It is concluded above that the bulk of the mesons must arise at about 70 km. | 
above sea level if Rau’s results are to be accounted for. At 70km. the pressure is. 
about 0:01cm. Hg, which is roughly 1/8000 of atmospheric pressure. This : 
means that primary particles have to create mesons in traversing an amount of air 
equivalent to only 1/8000 of the atmosphere. The corresponding cross-section | 
for the production process is considerably larger than seems likely if protons are’ 
the meson primaries. ‘The maximum conceivable cross-section in this case, in 
air, is of the order of the geometrical cross-section of the oxygen nucleus and | 
corresponds to a mean free path for the proton of about 1/20 of the atmosphere, 
so that the cross-section derived above is over 100 times the largest reasonable 
value for proton primaries. 

it has been suggested (Hoyle 1947) that the primary hard radiation is not) 
composed of protons but of very heavy nuclei of atomic number Z about 80. 
Since the energy a particle needs to penetrate the earth’s magnetic field is pro-. 
portional to its charge, such a nucleus would need energy about 101” ev. to reach the: 
equator. Its fate upon impact with oxygen and nitrogen atoms at the top of the 
atmosphere is unknown. It is probable that collisions with air nuclei willlead to.} 
rapid disintegration of the heavy nucleus but it is uncertain whether mesons will} 
arise directly from these disintegrations or only indirectly through the agency of 
protons and neutrons which will be liberated. For direct collisions involving|} 
meson interaction, the maximum cross-section for meson production by nuclei with |} 
Z about 80 cannot be much larger than about twelve times that for production by] | 
protons. (This estimate corresponds to the appropriate geometrical cross 
section, on the assumption that the nuclear radius is proportional to (atomic4 
weight)'). This is smaller than the cross-section required to explain the +S} 
variation by a factor of more than 10. Larger cross-sections for meson production} 
by a heavy nucleus might arise as a result of the Coulomb interaction, but it seems# 
very unlikely that for Z=80 the cross-section could be as large as that needed tol 
account for the +S effect. : | 

It is hard to see how this serious difficulty in the explanation of the + S effect 
may be resolved. The Pekeris theory of atmospheric oscillations appears to bell 
well-established ; Weekes and Wilkes believe that the results quoted in table 3 arelf 
substantially correct up to about 60km. On the other hand, for a semi-diurnallf} 
variation which is positively correlated with pressure, it is difficult to think of anif 
alternative explanation to that in terms of the Pekeris theory. The first need isi 
for more experimental evidence of a reliable kind on the existence and magnitude 
of this +S variation. Since the pressure oscillation decreases rapidly with in 
crease in latitude, reliable evidence is most readily obtained at the equator. Onef} 
of us (Sarabhai) hopes to obtain such evidence in India shortly. 


§4. THE —S VARIATION 

The experimental data suggest that, although the intensity measured under 4 
thick absorber shows a + S variation, that measured without an absorber shows af 
—S variation. ‘This implies that the easily absorbed radiation, consisting of the 
soft component and low energy mesons, is characterized by a —S variation which ig} 
large enough to swamp the +S variation shown by the energetic mesons. | 
If the soft radiation arises in the lower atmosphere as is certainly the case for 


the electronic component observed at sea level, it will tend to show a —S variation 
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on account of the ordinary mass absorption effect. However, even for the elec- 
tronic component it is difficult to say what would be the magnitude of the —S 
variation produced purely by pressure changes. ‘The effect on electrons arising 
from a primary electronic component could probably be estimated fairly well 
from the cascade theory, but the effect on the remaining electrons which arise from 
mesons and which probably form about 2 of the electrons at sea level, is extremely 
difficult to assess. It is likely that about equal numbers of these electrons arise by 
knock-on processes and by decay, and as their production is spread over various 
strata of the lower atmosphere, individual electrons may be influenced by different 
pressure changes in reaching sea level. In addition the point of origin of the 
mesons which produce the electrons is affected by pressure changes, again in a way 
depending upon altitude. In fact the position is so complicated that it seems 
almost impossible to calculate the resultant effect (which is known to be small) with 
any certainty. 

An alternative approach is from the experimental measurements of the day-to- 
day variation of intensity and the corresponding pressure variation, i.e. deter- 
mination of what is usually called the barometric coefficient. Here again no 
definite conclusion may be reached, for as already mentioned, it is uncertain how 
much of the observed barometric effect is a true mass-absorption effect and how 
much is produced by other factors, such as movement of meson-producing layers, 
which may be differently related to semi-diurnal and to day-to-day pressure 
changes. Duperier finds a barometric coefficient of 0-35°% per mm. Hg and 
deduces that about 60°% of this is a mass absorption effect and that the remainder is 
caused by movement of meson-producing layers; this result is dependent upon 
mesons being created at about 16km. which may not be the case. It does appear 
fairly certain, however, that the —S effect observed by Duperier with no absorber 
is too large to arise purely as a mass-absorption effect from pressure changes. 

Although the explanation of the —S effect is not clear, the suggestion made by 
Duperier (1945) that the main factor responsible is the atmospheric oscillation 
predicted by Pekeris is untenable when the magnitude of this oscillation is con- 
sidered. A negative S effect can only be produced by the Pekeris oscillation if 
mesons arise below 30 km. but the maximum displacement of isobars below 30 km. 
is only about 3 metres at latitude 50° and this movement is insufficient to produce 
an observable cosmic-ray effect. Duperier appears to be considering the vertical 
movement of individual particles, which indeed may be as large as 20 metres, as he 
states, but which is irrelevant in a consideration of the height of meson production, 
where the important quantity is the vertical movement of isobaric level. 

Further confusion on this point has arisen from recent publications by two 
workers in Ceylon, which apparently confirm Duperier’s conclusion that the —S 
variation can be accounted for by the Pekeris theory (Mailvaganam 1946, Kid- 
napillai and Mailvaganam 1946). It has already been mentioned (§ 3.1) that the 
results obtained by these workers for the movement of isobars are incorrect. 
Apart from this, their final result of 0-14°% for the amplitude of the —S effect seems 
to apply at \=0, whereas Duperier’s value of 0-18°%, with which they compare it, 
applies at A=50° and if produced by pressure oscillations would be four times as 
large at A=0. 

Soa Eee DP LUIRINA EV AR TAT LON 

From §2.1 it will be seen that although the D variation at high latitudes, 

recorded by different observers using similar apparatus, usually has an amplitude 
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about 0-2°/, the phase varies appreciably even at stations having the same latitude. 
For example, the results of Hess and Graziadei (1936).and of Forbush (1937),) | 


which both appear to be well established show a discrepancy of several hours in 
phase when allowance is made for the fact that with no external temperature 


correction the maximum observed by Hess and Graziadei would probably occur — 


about four hours earlier. Duperier (1941) observes a maximum as late as 174, but 


here the discrepancy may arise partly from the difference in thickness of absorber | 


used. Doan (1936) whose results are fairly convincing, places the maximum at 9". 


In fact all that can be concluded in general about the time of maximum intensity 


seems to be that it occurs during the daytime rather than during the night. 

It is possible that the recordings made with ionization chambers, even when 
shielded with 10cm. Pb, are affected by fluctuations in the gamma-ray content of 
the atmosphere. The latter, in some localities at least, appears to show a daily 
variation with maximum about noon. Forbush (1937) estimates that this factor 
might cause as much as one third of the D variation recorded by his apparatus and 
suggests that it might produce quite different distortion of the true cosmic-ray 


D variation in experiments performed in different localities. However, it appears — 
very unlikely that the wide phase variation recorded in the experiments mentioned | 


above can be entirely due to this effect. The results of Alfvén and Malmfors (1943) 
recorded in Stockholm with Geiger counter apparatus and supported by those of 
Kolhorster (1941) in Berlin, indicate that the D effect is appreciably different in 
different directions with the same zenith angle, and in view of this it does not seem 
improbable that the true D variation at a particular zenith angle will be different in 
different localities. 


The variation in character of the experimentally observed D effect has been 


ignored by most workers seeking an explanation of the effect; following the |] 
early results of Hess they usually assume that the maximum occurs at midday and |] 
try to account forthis. For instance, Vallarta and Godart (1939) suggest that such |} 
a variation can be accounted for at high latitudes by a heliomagnetic field and at low |} 


latitudes by fluctuations in the geomagnetic field. It is clear, however, that the 
D effect indicated by a wider survey of experimental data could not be mainly 
produced by either a heliomagnetic field or by geomagnetic fluctuations since it 


appears to show an appreciable variation in phase at different stations at the same |] 


latitude. Malmfors (1945) shows that such explanations are also incompatible 


with the results found by Alfvén and himself (1943) in Stockholm. These results |}} 
are particularly interesting as they also eliminate meteorological changes as the |]} 
main cause of the D variation. Malmfors suggests that occasional small distur- |} 


bances of the isotropy of the primary radiation in remote space may be responsible 
for the Stockholm results. | 


Before the cause of the D effect can be considered further, more experiments at | 
various latitudes and longitudes are needed to establish the variation of the phase | 


of the effect under a given absorber with geographical situation and also the vari- | | 
ation at a given station of phase and amplitude with thickness of absorber and J} 


with angle @. 
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The Neutrons Emitted in the Disintegration of Nitrogen 
by Deuterons 
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ABSTRACT. The photographic plate method has been used in a study of the neutrons 
produced in the disintegration of nitrogen by deuterons. 

A single group of neutrons, corresponding to an energy release of (5:15+0-10) mev., 
was detected, but the existence of the low-energy neutron group reported by Stephens, . 
Djanab and Bonner (1937) was not confirmed. The new results show that, if this group 
does exist, its intensity is not more than 15% of that of the principal neutron group, and that 
no excited levels of the 1°O nucleus less than 5 Mev. above the ground state are effective in 
this reaction. 

The neutron flux from the N--D reaction was compared with that of D+D neutrons 
produced by disintegration of deuterium present in the surface layers of the target employed. 
The data obtained were used in determinations of the angular distributions of the D+D 
and N--D neutron fluxes, and in tests of the reliability of the photographic plate method for 


relative flux measurements. 
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§1. INTRODUCTION 
"|e neutrons produced in the disintegration of nitrogen by deuterons 


according to the scheme : ; 
WN 42D->40 + in 

have been investigated by Stephens, Djanab and Bonner (1937), who used an 
expansion chamber and detected three groups of neutrons emitted from nitrogen 
targets bombarded by deuterons. They attributed two of the groups to the 
above reaction, and the effective Q values were found to be (1:1 + 0:2) MeV. and 
(5-1 +0-2) Mev. respectively. The third group was ascribed to the disintegration 
by deuterons of deuterium present in the surface layers of the targets employed, 
and the O value of this group (3-3 Mev.) was in agreement with the known energy 
release of the D + D reaction. 

Further investigation of the N+D reaction is desirable in order to confirm 
the existence of an excited level in the 420 nucleus, which is indicated by the 
presence of the low-energy group in the N+D neutron spectrum. In addition, 
the 5 Mey. neutron group may be useful in experiments where a neutron source 
of this energy is required, and for this purpose the characteristics of the neutron 
spectrum must be determined at various angles to the direction of the deuteron 
beam. 5 ; 

The use of photographic plates in the determination of neutron energies 
has been described by Powell (1940, 1942) and Richards (1941 a), who have studied 
the distributions in energy of the neutrons emitted in the bombardment of a 
number of light elements by deuterons. Their results were in general agreement 
with those obtained by the expansion chamber method. ‘The photographic 
plate technique has the advantage that plates may be exposed simultaneously 
at different angles to the direttion of the primary beam, so that information 
may be obtained concerning the angular distribution of the neutron yield 
(Richards 1941 b, Livesey and Wilkinson 1948). The photographic technique 
has been improved considerably as a result of the introduction by Ilford Ltd. 
of concentrated emulsions (Powell, Occhialini, Livesey and Chilton 1946), 
with which greatly increased accuracy of track measurements is attainable, 
and at present it is the most convenient and precise method for the determination 
of neutron spectra. This method was therefore adopted for a more detailed 
investigation of the neutrons emitted in the disintegration of nitrogen by deuterons. 


§2. EXPERIMENTAL METHOD 


A thick target of aluminium nitride was bombarded with a beam of deuterons 
accelerated through 930 kilovolts by the Cavendish Laboratory high tension 
equipment. 

Four light-tight boxes containing photographic plates of the Ilford C.? type 
with emulsions 100 microns thick were set up around the target so that the plates 
recorded neutrons emitted at 0°, 30°, 90° and 150° respectively. 

In order to obtain reasonably good geometrical conditions with the relatively 
large neutron source, the plates were placed 15cm. from it and were irradiated 
for 3 hours with an average beam current of 80ua. This arrangement gave a 
maximum spread of + 3° in the directions of the neutrons incident upon the plates. 
A stop could not usefully be employed to reduce the size of the N+D neutron 
source, as the stop itself acted as a source of D+ D neutrons. 
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A second series of plates was irradiated with neutrons emitted from a pure 
aluminium target bombarded with deuterons, in order to estimate the contribu- 
tion of the aluminium to the spectrum of neutrons from the aluminium nitride 
target. ‘The conditions of exposure were the same as for the nitride target, 
except that the target surface was cooled less efficiently and was therefore main- 
tained at a higher temperature by the beam. 

After the plates had been processed, detailed microscope measurements 
were carried out on the proton tracks occurring in the plates by two observers, 
using a Leitz microscope fitted with a 1/12 inch oil-immersion objective which, 
in conjunction with twin eyepieces of power x 10, gave an overall magnification 
of x1200. ‘Traverses were made across each plate in a direction parallel to that 
of the incident neutrons, and all proton tracks were recorded which originated 
in the area covered by the traverse, which were completely contained within 
the emulsion, and which made angles of less than 20° with the neutron direction. 
Between 300 and 500 tracks were recorded in each plate. The area examined 
was calculated from the total length of the traverse; it' was approximately 20 mm? 
in al] cases. 

For each selected track the total length @ projected on the focal plane of the 
microscope was measured, and the initial direction of the proton was derived 
from measurements of the angle of scattering projected on the focal plane « 
and the initial angle of dip . In order to obtain % it was necessary to correct 
the distances measured in the vertical plane for the effects of shrinkage of the 
emulsion caused by fixation. The correction consisted in multiplying all vertical 
distances by a factor s, which was assumed to be equal to the ratio of the initial 
to the final thickness of the emulsion. A series of experimental determinations 
of the shrinkage factor indicated s=2-0+0-2, and this value was used in all 
calculations of the angle ys. The actual angle of scattering y for each proton 
track was calculated, to the nearest degree, from the relation secy =secz secu. 

The proton range R was calculated from the equation R=asecy. This 
formula applies strictly only to tracks which are straight; in fact, some of the 
tracks were deflected from straight lines by nuclear collisions, but these deflections 
did not produce appreciable errors in the estimates of R and y. The proton 
energy E£,, was derived from the range—energy relation obtained for the B.1 type 
of emulsion by Lattes, Fowler and Cuer (1947). A comparison of the properties 
of the B.1 and C.2 emulsions has shown that the stopping-powers of the two 
types for «-particles are very nearly identical over a wide range of energies. 
It is therefore reasonable to assume that the stopping-powers for protons are the 
same: The neutron energy F, is related to the proton energy EF, by the formula 
E,=E, secty. The quantity £, was calculated for each individual track 


recorded. 


§3. NEUTRON SPECTRA 

The energy distributions of the neutrons emitted from the aluminium nitride 
target at the four angles of observation ($=0°, 30°, 90° and 150°), and from 
the aluminium target at 6 =0°, were derived from the numbers of tracks recorded 
in successive intervals of 0-1 Mev. inenergy E,. The data were divided, according 
to the magnitude of the scattering angle y, into three groups defined by the con- 
ditions : 

Ge Oy 97) 1G) 10 14 Oi)? ISP. 
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A typical track distribution is shown in figure 1, which refers to tracks measured. 
in the 150° nitride plate. 

It was found that, in the 
energy distributions obtained 
for 15°<y<19°, the D+D 
and N+D neutron groups 
were broader than those 
occurring in the distributions 
obtained for low values of y. 
This effect was due to small 
errors in definition of the 
directions of the incident 
neutrons, caused by the use 
of a comparatively large 
neutron source. ‘This “ang- 
ular straggling” is always 
more pronounced in the case 
of protons scattered at large 


angles, and measurements on 5 
these tracks are less reliable for 0 
the determination of neutron 0 | 2 3 4 S 6 


i ergy in Mev 
energies. Despite the strag- Energy e 


gling effects the principal Figure 1. | Number of tracks plotted against neutron 
energy for 150° plate: 


neutron groups were clearl ; 
coe 2 <y<9°; (ii) 10° <y < 14°; (it) 15° Sy <1. 


resolved in the distributions  ° 
obtained for y<20°, and no x 
difficulty was experienced in estimating the number of proton tracks which fell 
within each group. The data for 0°<y<9° were used in all energy determina- 
tions, since the effects of angular straggling were very small in these cases. 
Since the total number of tracks measured in each plate was comparatively 
small, statistical variations in the number of tracks at each energy were considerable. 
The observed distributions were therefore smoothed by the addition to each 
ordinate on the energy scale of the mean of the two neighbouring ordinates. 
In order.to obtain the actual neutron energy distributions, the smoothed experi- 
mental distributions were corrected, firstly for the variation of the total neutron- 
proton scattering cross-section with energy, the formulae of Kittel and Breit 
(1939) being used, and secondly for the escape of protons from the emulsions. 
The correction for the latter effect was calculated geometrically on the assumption 
that the proton range varies with energy according to the Geiger rule, which, 
for a fixed neutron energy E,, yields the following relation between the proton 
range R and the scattering angle y: R=R,cos*y where R, is the range of a 
proton with energy E,. The fraction f of the total number of protons, scattered 
in the range y; <y <j, which is recorded in an emulsion of thickness 2d, is given by 


f=l- 
Ro (2 =) + 4(sin 27, —sin 271) —}(sin 4yy —sin 41) — (sin 6y2 —sin 6y,)} 
8nd (sin? yp —sin? y,) ; 
provided that Ry cos* y, siny,<2d. Although the Geiger rule does not accurately 
represent the range-energy relation for protons, the correction factor is not 


Neutrons emitted in the disintegration of nitrogen by deuterons car 


appreciably affected by the 
departures from it which 
_ actually occur. The cal- 
culated factor is plotted as 
a function of the neutron 
energy in figure 2 for the 
three ranges of y used in 
the experiment. f 
The corrected neutron 
spectra for the aluminium 
nitride target at four angles 
of observation and for the 
aluminium target at 6=0°, 
are shown in figure 3. 
The four nitride spectra 
show two principal neutron 
groups, the mean energy 
of the one varying between 1-9 Mev. at 150° and 3-5 Mev. at 0°, and that of the other 
lying between 5 Mev. and 6 Mev. at all angles. In addition, a third group below 
1 Mey. was detected in the 0° spectrum, the investigation of which was extended 
to. lower energies than in the other cases. ‘This group is ascribed to the disinte- 
gration of carbon present on the target surface, according to the reaction : 


2C+7D—jn + iN --0-27 Mev. 


“0 l 2 3 4 5 6 7 8 
Neutron energy in MeV. 


Figure 2. Escape fraction f as a function of neutron energy : 
G07 <9; (1) 107 yol4" 5 (i) 15° 7 < 197: 


A similar group was found in the 0° aluminium spectrum. ‘The occurrence 
of the #C +D neutron group indicates that high-energy neutrons were emitted 
also as a result of the disintegration of the isotope #C : 


1C +2D—jn + 17N +5:5 Mev. 


It has been shown by Bonner and Brubaker (1936) that the abundance of the 
13C +D neutron group is only 1% of that of the C+D group. It is clear, 
therefore, that the #3C + D reaction cannot account for more than a small propor- 
tion of the 5-6 Mey. neutrons occurring in the nitride spectra, and these may be 
identified with the N+ D neutron group investigated by Stephens, Djanab and 
Bonner (1937). Nevertheless, a small neutron component was found near 
6Mey. in the 0° aluminium spectrum, and this may be ascribed to the ®"C+D 
reaction. 

The remaining principal group in the nitride spectra was produced by the 
disintegration of deuterium present in the surface layers of the target, as may be 
seen by comparison with the D+D neutron spectra obtained by Livesey and 
Wilkinson (1948), who deduced a value of 3-23 + 0-02 Mev. for the energy release 
of the reaction. The shape of the D+D neutron group in the 0° aluminium 
spectrum is different from that of the D+D group in the 0° nitride spectrum, 
and this is probably due to the high temperature at which the aluminium target 
was maintained: under these conditions the material closest to the surface 
of the target retained only a very small amount of deuterium, and so the yield 
of neutrons of maximum energy was low. 

Between 0:8 Mev. and 2-0 Mev. no clearly defined group was detected in the 
nitride spectra; this result disagrees with those of Stephens, Djanab and Bonner. 
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A comparison of the 
nitride and aluminium 
spectra at 0° shows 
that the distributions 
of neutron energy be- 
tween ~ 0:8mMev. and 
2 UME Vie> “aver avery 
nearly the same in the 
two.cases. ‘The ratio 
of the neutron fluxes 
in the two spectra (Al 
flux/AIN flux) was 
determined directly, 
and the values ob- 
tained for the energy 
ranges 1-0-1-5 Mev. 
and 1-5-2-0 Mev. were 
1:-6+0-2 and 1:-4+0-2 
respectively. 

The fact that the 
aluminium target pro- 
duced more neutrons 
in the energy range 
1-0-2-0 Mev. than the 
aluminium nitride 
target shows that ni- 
trogen did not contri- 
bute appreciably to 
this group of neutrons. 
It is known that the 
carbon isotope #C 
emits neutrons of this 
energy when bom- 
barded with deute- 
rons, but the abund- 
ance of the C3) 
group, as calculated 
from the results of 


Bennett, Bonner, _. : 2 " : , 
Hudepctim eRicha de Figure 3. Corrected neutron spectra, with ordinates in arbitrary units : 

ee ; (a) AIN target at 150°, (b) AIN target at 90°, (c) AIN target at 30°, 
and Watt (1941), is too (d) AIN target at 0°, (e) Al target at 0°. 


0 2 3 4 5 6 
Energy in Mev 


small to account forthe 


observed effect. It is possible that the aluminium present in both targets produced 
most of the neutrons in this energy range. If this hypothesis is correct, the ratio 
of neutron fluxes, as calculated from the chemical compositions of the targets, 
should be 1-5, which is in agreement with the experimental value. An upper 
limit may therefore be set to the abundance of the low-energy N+D group 
reported by Stephens, Djanab and Bonner (1937). On the assumption that 
the maximum variation in the neutron flux ratio is three times the probable 


Neutrons emitted in the disintegration of nitrogen by deuterons 529 


error, this upper limit is found to be 15°, of the abundance of the high-energy 
N+D group at the same angle, namely 0°. This same conclusion is valid for 
all four angles of observation, and it disagrees entirely with the results of Stephens, 
Djanab and Bonner, who estimated the abundance of the low-energy group: 
to be 150% compared with that of the high-energy group, at an angle of 90° 
to the direction of the deuteron beam. 

The high-energy N+ D neutron group is well-defined and is clearly resolved 
from other components in all the neutron spectra investigated. The mean 
thick-target energy varies between 5-1 Mev. at 150° and 5-7Mmev. at 0° when 
the bombarding energy is 0-93 Mev. 


§4. THE ENERGY RELEASE OF THE N+D REACTION 


In order to obtain the energy release of the N+D reaction, it was necessary 
to extract, from the observed thick-target spectrum for each of the four angles 
of observation, the corresponding thin-target energy. A careful examination 
of the spectra of figure 3 showed that the effects of both angular straggling and 
range straggling of the proton tracks could be represented analytically by the 
superposition of a generalized Gaussian error function on the actual thick-target 
neutron spectrum. A method described by Livesey and Wilkinson (1948) 
was therefore employed to derive the thin-target energy for each spectrum. 
The N+D spectrum at each angle was re-plotted, with the same corrections 
and smoothing procedure as those previously described, but with the neutron 
energy specified at 0-05 Mev. intervals, and the integrated spectrum was con- 
structed in each case. ‘The correction required for estimating the thin-target 
energy (£,) from the value given by extrapolation of the straight part of the 
integrated spectrum was obtained by comparing the shape of the experimental 
spectrum with those of a set of theoretical curves calculated for different values 
of the Gaussian straggling parameter. ‘The extrapolated energies, the derived 
thin-target energies, and the corresponding values for the energy release QO are 
included in table 1. 


Table 1. Extrapolated energies, thin-target energies and 
calculated Q values in Mev. 


Angle ¢° .. ee ea at aE 0 30 90 150 

Extrapolated energ i: ae a 6-00 5-94 5-82 5:39 
Correction im ye A ..  —0:08 —0-10 —0-13 —0-10 
Thin-target energy, E, .. te re Sy 5-84 5-69 5-29 
O value SS nt aE ts 5:07 DAMS el 5-20 


These figures were checked by independent determinations of FE, from the 
spectra corresponding to angles of scattering y between 10° and 14°, and no 
discrepancies were found. If the four estimates of the Q value are of equal 
weight, the mean value is 5-15 Mev., and the probable error, derived from the 
above results only, is +0-06Mev. The uncertainties in the range-energy relation 
of Lattes, Fowler and Cuer (1947) must also be taken into account; the probable 
error in their data for energies between 5 Mev. and 6Mev. is +0-07Mev. In 
addition, the error of +10°% in tne shrinkage factor of the emulsions leads to a 
probable error of + 0-03 Mey. in the final Q value. ‘These errors in combination 
give a total probable error of +0-10Mev., and the energy release is therefore 


quoted as Q=5:15 + 0:10 Mev. 
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This result agrees with the value of 5-1+0-2Mev. derived by Stephens, 
Djanab and Bonner (1937), and it may also be checked by calculation of the 
maximum energy of the positrons emitted by the oxygen isotope, ’O. The 
(d, p) reaction in nitrogen has been studied by Cockcroft and Lewis (1936) and by 
Holloway and Moore (1940), and the most accurate results give a value of 
8-55+0-08mev. for the energy release. Comparison of the two reactions : 

MAN +?2D->1H +13N + 8-55 + 0-08 Mev. 

AN +2D-> in +180 +5-15 + 0-10 Mev. 
yields the relation 30—%N= —(jn—{H)+3-40+0-13mev. The neutron— 
hydrogen mass difference is equivalent to 0-76Mev., according to Mattauch 
(1942) and Stephens (1947), hence the difference in the masses of 430 and 47N 
is equivalent to 2-64+0-13 Mev. ‘This exceeds the maximum positron energy 
from 1°O by an amount equivalent to the mass of two electrons, that is, 1-02 Mev., 
provided that the mass of the neutrino is negligible. ‘The maximum positron 
energy is therefore 1-62 + 0-13 Mev., and this result agrees with the experimental 
value of 1:72 Mev. deduced from the data of Konopinski (1943). 

The results shown in figure 3 did not confirm the existence of the excited 
state of 150 which was indicated by the results of Stephens, Djanab and Bonner 
(1937). A detailed examination of the nitride spectra, and especially of the 
150° spectrum, leads to the conclusion that the energy release of the N+D 
reaction is single-valued and that no excited levels of the *O nucleus less than 
5 mev. above the ground state are operative in this reaction. 


§5. ANGULAR DISTRIBUTIONS OF NEUTRON YIELD 


In making relative measurements of neutron flux by the photographic plate 
method, it is assumed that the number of proton tracks occurring per unit area 
within a defined solid angle, is, after correction for the escape effect, proportional 
to the product of the neutron flux and the total neutron—proton scattering cross- 
section. ‘This assumption is valid if the angular distribution of neutron—proton 
scattering is independent of the neutron energy in the range investigated, and 
if the concentration of hydrogen in the emulsion is constant. That the latter 
condition was satisfied in this experiment was shown by the fact that no abnormally 
large variations in the abundance of proton tracks were detected in examining 
different parts of the plates. The angular distribution of neutron-proton 
scattering has been determined experimentally by Dee and Gilbert (1937), by 
Bonner (1937) and by Champion and Powell (1944) for neutron energies below 
10 mev., and it is generally accepted that the distribution in the centre-of-mass 
coordinate system is effectively isotropic within this energy range. 

The reliability of the method of observation and of the escape correction used 
in this experiment may be checked by comparison of the observed angular distri- 
bution of proton tracks with that corresponding to isotropic scattering in the 
centre-of-mass coordinates. “ Provided that the effects of angular straggling 
may be neglected, the number of proton tracks recorded in the range of angles 
defined by y;<y<y-2 should be proportional to the integral 

Y2 

‘ 2 cosy sin ydy =(sin? y, —sin? y,). 
The experimental data were tested by dividing the total numbers of proton 
tracks in the D+ D and N+D groups, corrected for escape, by the function 
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(sin? y, —sin? y,), to see whether the same quantity was obtained for the different 
angular ranges. ‘The results are shown in table 2. 


Table 2. Ratio of experimental and theoretical proton abundances, 
in arbitrary units 
Spectrum OF S7<=9° “10° Sye14° 15°<y)<19° 


0° 3543 41+3 31+3 
30? 3343 2442 22-2 
90° 18+2 2142 17+2 
150° 24-2 28+2 Do eee 


The 0°, 90° and 150° results are internally consistent, but the 30° figures show 
possible deviations from isotropic scattering. ‘These deviations are probably 
not significant, but they may be regarded as an indication that the 30° results 
are less reliable than the others. 

In the four neutron spectra obtained with the nitride target, the N-+D and 
D+D neutron groups were everywhere resolved, and the number of neutrons 
detected in each group was readily estimated. In order to test more stringently 
the accuracy of the escape correction at the higher neutron energies, and also 
to detect whether there was any tendency on the part of the observers to count 
too many long tracks at high scattering angles, the ratio of the D+D flux to the 
N+D flux was determined for the three ranges of the scattering angle y. In 
the calculation of this ratio, the mean scattering cross-sections for the two neutron 
groups weie derived from the data of Kittel and Breit (1939). The results 
are shown in table 3. 


Table 3. Ratio of D+ D and N+D neutron fluxes 
Spectrum OFS 9° NOB SiS 15 ry 98 Mean Ratio 


0° 4140-9 4840-8 4-4+0°8 4-4 40-5 
30° 3-4+40-6 4341-0 59413 4-4 4055 
90° 1:7-+0:3 1-8-+40°3 1340-2 1-5040-13 

150° 1-8+40-2 2-5-40:3 2-5-40:3 2-120-18 


The internal consistency of these results shows that the method of observation 
and of correction for escape did not give rise to serious errors. 

The consistency tests already described indicate that any errors occurring 
when the photographic plate method is used for comparing neutron fluxes can 
be reduced to low values by the methods used in this experiment. ‘The data 
were used to determine the angular distribution of the D+D neutron flux in 
order to compare the results obtained with those of Allen, Livesey and Wilkinson 
(1948), who determined the angular distribution of neutrons emitted from a 
thick deuterium target bombarded by deuterons of energy 0-93Mev. The 
angular distribution was obtained by estimating the neutron flux at each angle 
relative to that at 0°, and the results are shown in table 4. 


Table 4. Angular distribution of D + D neutron flux 
Angle ¢° 0° <¥<9° 10°<y<i4° 15°Sy<19° 0°<y<19° 


0 1:00-0:08 1:00-0:07 1:00-40-07 1-:00-+0-05 
30 1-15=-0-12 0-76-+0:07 0:90 -+0-09 0:90-+0:05 
90 0:29-0:04 0:28+0-03 0:26-+0-03 027220:02 

150 0:62+0-06 0:45 -+0-06 0:54-40:06 0:52+0-04 


The results for the angle range 0°<y<19° may be compared with the experi- 
mental curve of Allen, Livesey and Wilkinson in figure 4. 
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In general the agreement be- 
tween the two sets of results is 10 
close. 

The above results show that 
it is possible to make reliable 
determinations of angular dis- 
tributions of neutron flux by the 
photographic plate method. The 
angular distribution of the N + D 
neutron flux was calculated from 
the data referring to the entire 
range of scattering angles be- 
tween 0° and 19°; for ¢ equal 


fe} 

° 

fe} 
—— 


Neutron flux relative to flux at O 
Oo 
wn 


to 0°, 30°, 90°, 1507 the: values vo a 2OmnEAO m6 go 100 «20 «4 0 180 
obtained were 1:00+0-10, Angle in degrees 
0-91 + 0-10, 0-90 + 0-08, Figure 4. Angular distribution of D+D neutron flux 


. hick t ts 
1-37 + 0-15 respectively. : from tie se an 
The 150° fi ee ehat Experimental curve of Allen, Livesey and Wilkinson, 

; 2 AUER Vee SES with 3 photographic plate results for the angle 
higher than the others, but the range 0<y < 19°. 
statistical accuracy is lower than 
that in the D+D distribution, and the deviation from an isotropic distribution 
of the neutron yield is hardly significant. 
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Dissociation Energy of the NO Molecule * 
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ABSTRACT. During a new investigation of the NO band spectrum many bands of the 
B and y system were photographed and analysed, and the rotational constants evaluated. 
Intensity drops and perturbations were found and two predissociation limits at 53800 and 
50 500 cm were established. In this way the dissociation energy of NO is fixed at 4:29 ev. 
The x*II ground state of NO dissociates in the N(?D)+O(P) atomic term combination. 
The new dissociation scheme is strongly supported by the results of photochemical experi- 
ments. 


many bands of the f, y and « system were photographed with great intensity 

and resolving power. ‘The rotational analysis of the « system was given 
in a preceding paper (Gero, Schmid and Szily 1944); after remeasuring the whole 
of the NO spectrum between 1950 and 2700 a. some new data for the 8 and y 
bands can also be given. 

Light source, spectroscope and experimental details were described in 
the paper mentioned above. ‘The rotational analysis of the following y bands 
(a?x —x°I1) was carried out up to quite large rotational quantum numbers: 
(0,0), (0,1), (0,2), (0,3), (0,4), (0,5), (1,0), (1,1), (1,3), (144), (155), (1,6), (2,0), (2,2), 
(2,3), (2,7), (3,1), (3,4), (3,5). Wave numbers of the (1,0), (2,0) and (3,1) bands 
are given in the tables 1-3. Some bands of the B-system (B?II — x?I) were also 
present between the much stronger y-bands; the (1,5), (1,6), (2,3), (2,4) and 
(4,2) B-bands were analysed up to K ~35-40, ie. rather high rotational quantum 
numbers. When photographed in active nitrogen, which is especially favourable 
for exciting the B-bands, they are observable only up to K=25-30 (Jenkins, 
Barton and Mulliken 1927). 

Averages of the upper-state combination differences of all bands with common 
initial states were evaluated. The rotational constants have been determined 
by the usual graphical method, which yields the B,, constants given in table 4. 
The data are not accurate enough to determine the variation of the D,, constants 
with v’; D’=6-0. 10-* cm~! was therefore used throughout. 

The vibrational terms included in table 4 are the heights of the first rotational 
levels with.J = 4 of the different rotational term series above the x*I1, (v=0, J =) 
ground level, determined graphically using all bands and branches. The 
1otational constant and the height of the first rotational level of the pD?&, v=0 
state, evaluated from the v’ =0, «-bands, are also given in table 4. As can be 
seen, the B,’ of the e-bands is greater than each of the B’-values of the y-bands ; 


Ip the course of a new investigation of the band spectrum of the NO molecule, 


* This paper was written by L. Ger) in 1943. In early 1945, during the period after the siege of 
Budapest, the first copy of the manuscript went astray, and figures 1 and 2 had to be reconstructed. 
In figure 2 we were able to mark only schematically the places where the accidental predissociation 


effect was observed.—J. G. VALATIN. 
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Table 1 
A= Or band AX —xX*7T1, ; 1,0 band 

J SRoy R, OF P; J R; QO, Py OP i. 

$ =46552:28 46544-61 46540-51 

14 rgeR.80 547-66 539-28 46535:57 14 46439-04 46427-87 46419:77 46415-95 
24 566-48 551-09 538-73 531-18 24 446-70 431-12 419-08 411-21 
34 574-70 554-97 538-73 527-39 34 454-34 434-81 419-08 407-16 
44 583-36 559-61 539-28 524-27 43 462-26 439-04 419-08 403-31 
53 592-20 564°85 541-10 521°58 54 470-97 443-83 419-77 400°31 
64 602-34 570-91 543-13 519-40 64 480-04 448-97 421-06 397-60 
7 612-57 577-11 545-50 517-84 74 490-00 454-34 422°91 395-42 
84 623-42 584-21 548-57 517-08 84 500-19 460:°71 425-13 393-87 
94 634-88 591-52 552-28 517-08 94 510-46 467-19 427-87 392-45 
104 646°85 599-84 556-44 517-08 103 521-58 474-36 431-12 391-87 
ie 659-67 608°51 561:-45 517-84 113 533-10 482-01 434-81 391-87 
124 673-24 617-92 566°48 519-40 123 545-50 490-00 439-04 391-87 
134 686-40 627°87 572:56 521-58 134 557-81 498-91 443-83 392-45 
144 700-99 638-34 579-19 524-27 144 570:91 508-00 448 -97 393-87 
154 716°11 649-35 586-44 527-39 153 584-21 517-84 454-88 395-81 
163 7S lesa 660-87 594-04 531-18 164 598-62 528-02 461-23 398-25 
174 747-65 673 :24 602-34 SOD 7 174 613°55 538-73 467-92 401-22 
183 764-47 686-07 611-05 540°51 184 628-35 549-83 475-22 404-58 
194 781-62 699-27 620-53 546-02 1924 643-97 561-45 483 -06 408-54 
203 799-31 TAB 28 630-66 552-28 204 659°67 573-99 491-26 413-00 
214 817-79 TOUT 641-19 558:89 214 676:68 586-44 500-19 418-00 
224 836-76 743-06 652-35 566-48 224 694-01 599-84 509-55 422-91 
234 856-39 758-50 664-16 573-99 234 711-66 613-55 519-40 429-21 
243 876:47 774-71 676-68 582:34 244 729°87 627:87 529-68 435-50 
254 897-11 791-50 689-26 591-52 254 748-40 642-82 540°51 442-90 
264 918-24 808-76 702-67 600:°75 264 767-64 658-01 552-28 450-11 
273 939-74 826-74 716:55 611-05 274 787°47 674-13 564:-14 458-40 

1 


283 962:12 845-01 731-31 621-64 284 807-64 ~ 690°30 576°63 467-19 
293 985-18  864:09 746-32 633-07 293 828:36 707-05 589-46 

303 47008:58  883:59 761-97 644-67 304 849-60 724-54 602:°92 485-53 
313 032-46 903-87 778-21 657-32 313 871-39” 742-36" 617-169 49558 
323 057-01 BEES WM Srilil 670:21 323 893-72) 100.84 mS 136 506-40 
333 082:05 945-58  812:44~ 683-68 334 -916°57 779-80 646-85 517-84 
343 107:87 967-26 830-40 697-63 344 939/45 e799 3H 662-17 a5 29-63 
353 133-96 989-66 848:90 712-07 353 963-51 819:35 678-60 


363 CMOS XG Pichi 363 Seleia — wexsSSii SOSA Ns 554-97 
373 035-93 887-52 743-06 37% 47012°55 860-87 712:07 567-91 
383 059°84 907-67 = 759-44 384 038-08 882-31 129372) DS1c65 
393 084:53 928-37 776-39 393 063-94 904-48" 748-405 595-9038) 
403 109-46 949-50 793-88 403 090°34 = 927-25 767-01 611-05 
413 1351 Seo lees 811-65 413 117-33 950:12) =" 786:3 5 026.59) 
423 161-27 993-57 830-40 | 424 144-74 973-83 806-15 642-82 |} 
433 187-93 47016-49 849-60 433 PZ. SkSEO} $26:74 "1659 -67—) 
444 215-34 040-09 868-99 44% 201-27 47022-48 847-56 676-68 |} 
453 243:07 063-94 = 889-35 453 230-20 047:93 868-99 694-01 
464 271:35 088-47. 910-14 | 464 259] Se eOUS dS 890-82 
474 300-20 =113-59 931-59 | 474 290-00 2099-995 913-23 

483 329:47 139-39. 953-10 |, 484 320°37 = 126-83 936-39 

493 359-64 =165:-49 975-51 493 351-48 154-21 959-80 

504 390105 192723> 9998-08 503 3832125 182-01 983-81 

514 421:32 219-58 47022-48 513 415-32 210-39 47008-:58 

524 452°83 247-36 046-07 523 AF 799 2395 ie 03879 

533 485-00 275-69 070-78 533 481-14 268-80 = 059-84 

543 517-63 304-65 095-97 543 514-78 298-66 085-79 

553 550°85 334-12 «= 121-77 553 549-29 329-47 ISD 

564 584-48  364:17 148-08 563 583-78 360-36 139-39 

57% 618-93 394-78 174-82 573 618-93 391-86 167-76 

585 653-45 425-89 . 202-01 583 654:°90 424-11 196-04 

594 689:02 457-43 230-20 593 C912 OR 45 Os 71 OOH, 

604 724:90 489-62 258-46 603 728:20 489-62 254-37 

613 760-94 522-46 287-60 613 765-44  523-59° 284-65 <= 

623 JQ8"2) 595-02 eesiyel 9 623 803:37 557-82 314-93 

634 835-75 589-44 347-49 633 841:30 592-66 346-05 

644 873-78 623-69 378-74 643 880:47. 627-70 377-66 

654 912-34 658-95 409-50 654 919:90 663-65 409-50 

663 951-20 694:07 440-88 665 669:62 442-39 

67% 729°98 473-14 67% 4D) 

683 505-60 
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Table 2 

A? —x*I1,; 2,0 band Ao Xe 2.0) bang 

J SRoy R, O1 P, J Ry OQ» P, Py. 
2 48861-73 48854-14 48850-11 

1d 868-28 856-62 849-08 48844-83 14 48737-33 48729-52 
23 875-42 860°16 848-29 839-85 23 739-90 728-59 
34 883-30 863-88 848-29 836°10 33 48763-05 743-46 727:90 48716-22 
4} 891-89 868-28 849-08 832°78 44 771-40 747-43 727-90 M235 
54 900-71 873-82 850-11 830-33 53 779-30 752-00 728°59 709-08 
65 910-23 879-70 851-92 828-06 63 78839 756:89 729-52 706:28 
7s 920-38 885-05 854-14 826-39 7% 797-35 762:27 731-08 703-67 
83 930-65 891-89 856-62 825-08 83 807-06 768-12 733-01 701-53 
9s 941-81 898-97 860-16 824-03 94 817-50 TE oy 735-31 700:75 
103 953-43 907-00 863-88 824-03 103 828-05 781-23 738-33 699-16 
114 965-54 Bil So2)| 868 °28 825-08 113 839:23 788-39 741-68 698-87 
123 978-40 923-83 873-08 826:39 123 850-51 796-20 745-43 698-87 
134 991-93 933-09 878-77 828-06 133 863-09 804-42 749-76 699-16 
144 49005-38 943-06 884-67 _ 830-33 | 144 875-42 812-97 754-42 700-75 
153 019-58 953-43 891-10 832-78 Sp 888-09 822-06 759-70 701-53 
164 034-56 964-58 898-08 836-10 | 164 901-60 831-54 765-33 702-93 
lo 049-88 975-97 905°85 839-85 173 OIA | 841-61 771-40 705-07 
183 065-75 988-14 914-09 844-14 183 929°-911 851-92 778-21 708 06 
194 082-25 49000:79 922°85 849-08 193 944-66 863-09 785°23 711-33 


20 099-36 013-84. 932-04 854-14 | 20$ 960-07 874-54 792-70 714-94 
214 117-04 027-49 «941-81 860-16 | 214 975-97 886-34 = 800-75 719-19 


224 041-99 952-24 866-77 | 22% 991-93 898-97 809-21 723-63 
234 056-48 963-10 873-82 | 234 4900867 911-76 818-34 728-59 
244 071-86 974-25 881-32 | 244 025-73 924-98 828-06 734-69 
254 087°55 986-20 889-28 | 254 043-63 938-87. 837-65 
264 103:55 998-85 898-08 | 264 061-63 953-43. 848-29 
273 120-71 49011-58 907-00 | 274 080-39 967-89 859-25 
284 137-99 025-28  916:86 | 283 099-36 983-15 870-50 
294 155-56 039-30 926-87 | 294 119-09 998-85 882-32 
303 174-31 053-80 937-67 | 303 138-98 49015-02 894-78 
314 193-09 068-87 949-05 | 314 159-35 031-79 907-58 
324 212-61 084-56 960-75 | 324 180-12 048-94 920-93 
334 232-31 100-39 972-86 | 33 20176 066-56 934-82 
344 252-63 117-04 986-20 | 344 223-75 084:56 949-05 
354 273-65 134-15 998-85 | 354 246-30 103-55 963-91 
364 295-18 152-12 49012-71 | 363 269-72 122-54 979-21 
374 317-13. 170-20 027-49 | 374 142-24 995-01 
384 339-13 188-78 041-99 | 384 316-40 49011-58 
394 362-66 207-57 057-66 | 393 340-45 182-71 

404 385-88 227-80 073-26 | 404 365-27 203-96 045-36 
414 410-36 248-00. 090-13 | 414 390:54 225-37 063-07 
424 434-71 268-96 106-98 | 424 41651 247-51 081-21 
434 459-95 290-00 124:52 | 434 442-58 269-72 099-36 
444 485-14 311-69 142-24 | 444 469-74. 292-77. 119-09 
454 511-26 334-12 161-14 | 454 497-27 316-30: 138-98 
464 537-63 356-74 180-12 | 464 525-91 340-45 159-35 
474 564-89 380-17. 199-82 | 473 55352 365-27 180-12 
483 592-94 404-01 219-68 | 48% 582-74 390-54 201-22 
- 493 621-05 428-19 240-62 | 493 61231 416°51 222-93 
503 649-14 452-98 261-03 | 504 642-36 442-58 244-94 
514 678-99 478-38 282-42 | 514 673-32 469-74 267-63 
521 504-47. 303-66 | 524 497-27 

533 325-65 


this is a strong argument against the hypothesis that the y and « bands belong 
to the same band system. It must also be determined whether there are homo- 
geneous perturbations in one or in another of these bands which deform the 
rotational constants. From figure 1 it is seen that the B’—B” curves (Gero 1935 a) 
for the bands y (0,2), (1,0), (2,0) and (3,1) and ¢ (0,3) are nearly horizontal straight 
lines which means that no considerable perturbations are present in thes¢ bands. 

Predissociation phenomena were found in the upper state of the y-bands. 
‘The rotational term series show sudden intensity falls at v=0, K=74; v=1, 
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Table 3 
Ax —x?I1,; 3,1 band a?d —x?I1,; 3,1 band 
J R, OQ, P, J Rs Q2 Pz 
$ 49254-40 
13 257-02 49249-43 13 49137-99  49130-29 
De 260-52 249-43 23 140-48 129-01 
34 264-21 249-43 $ 144-22 129-01 
43 268-96 249-43 4h 147-99 129-01 
54 273-65 DEOS 49231-53 $ 152-12 129-01 
$ 278-90 2A 228-91 63 157-49 130-29 
74 284-87 254-40 227-80 4 162-95 131-98. 
83 291-43 257-02 226-05 84 168-61 133-09 
$ 298 -46 260°52 225-37] $ 174-82 135-99 
103 306 :34 264:21 225-37 104 181-60 138-98 
113 314°75 268-96 226-05 114 188-78 142-24 
123 323-60 273-65 227-80 124 196-36 146-23 
134 332-80 278-90 228-91 134 203-96 150-41 
143 342-66 284-87 Doles 144 212-61 155-56 
153 353-04 291-43 PPT) 154 49287-60 222-05 160-44 
163 363-74 298-46 236-91 164 300-74 PENISS) 166-09 
173 375 :34 305-95 240-62 174 314-75 241-45 172-23 
184 PSI 314-14 244-94 183 328-89 252-12 179-01 
193 399-76 SL27 3 249-43 194 343-66 262:61 185-90: 
203 412-99 331-93 Dayal: 204 358-67 274-16 193-09 
214 426-45 341-49 261-03 214 373-98 285-92 201-22 
225. 440-77 352-14 267°63 224 390:54 298-46 209-70: 
234 454-86 362-66 274-16 235 407-11 310-98 218-81 
243 469-74 373-98 281-98 244 424-50 324-29 228-91 
254 485-14 385-88 290-00 254 441-68 338-09 238-15. 
264 501-75 397-91 298 -46 264 459-95 352-14 248-00 
274 518-74 410-99 307-58 274 478-38 366°88 259-38 
283 535-99 424-50 Silos! 284 497-27 382-15 270°62. 
293 SES 438-28 327-50 293 516°56 397-91 282-42 
304 O72 42 452-98 338-09 303 413-77 295-18 
314 590-71 467-22 348-97 314 430-31 307°58. 
324 609-77 483-03 360:67 324 447-21 320-73 
334 629-99 498 -84 373-03) 3350 464-62 334-12 
344 650-01 515-90 385-88 343 483-03 348-97 
354 670-63 532-78 398 -96 354 501-75 363-74 
364 692-03 550-25 412-99 363 380-17 
3743 568°31 427-50 
383 441-68 
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Table 4 
Vibrational terms 
v=0 44199-2 
2341-4 
1 46540-6 
2309°9 
A2D 2 48850°:5 
2276:7 
3 51127:2 
D> v=0 53291:-9 


e040 U2 Ki 2 Uo, = 38: 
joining these points in figure 2, a 
limiting curve of predissociation can 
be drawn, which intersects the ordinate 
axis at a value of about 53800 cm“. 
An N +O atomic term combination must 
lie at this energy, which gives rise to the 
predissociation phenomena observed. 
The upper-state rotational term series 
of the 6-bands intersect in many places 
the much steeper upper term series of 
the y-bands. In the neighbourhood of 
the intersections perturbation phenomena 
should be discovered. ‘Tables 5 and 6 
contain parts of the corresponding 
branches of the (4,2) 8-(2,2)y and 
(2,4) B-(1,4) y bands respectively, which 
belong to common lower states. At 
the intersection points of the upper states 
the differences of corresponding wave 
numbers change their signs; these places 


1o’cm 


SEi0f 
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0 1000 


2000 3000 4000 5000 
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in the tables are marked with arrows. Figure 2. 
Table 5 

J ae Joe ee ioe ee Ate, 
214 45519-76 45320:10 45428-05 45153-05 45440°58 45269-54 
224 495-79 335-97 399-50 161-19 416-22 287-42 
234 469-47 352-28 369-32 169:72 390:21 305-87 
244 442-75 369-32 338-31 178-49 364-38 324-56 
254 414-45 386-66 30587 188-60 336°77 aa 344-22 
264 384-74 ie 404-90 272-49 199-12 307-58 364-38 
274 354-48 423-63 237-61 210-06 277-62 384-74 
284 S237, 442-75 202-43 ae 221-63 ' 246:49 405-79 
294 290-44 462-84 165-15 233-76 
304 483-19 127-00 246°49 
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marked with squares in figure 2, full lines show where both intersecting term 
series were observed, and broken lines, where only the upper series of the y-bands 
were observed. ‘There are, however, three points where intensity anomalies 
were observed in the branches: of the y-bands; these points are marked with 
circles. ‘Two circles lying close together in a term series mean that as a result of 
intersection with II-states, the F, and F, perturbations of the a2 term series 
do not take place at the same J-values. These points also lie near the y-8 upper- 
term intersections, but at higher energies than those intersections at which no 
effect could be detected. Between the circles and squares a limiting curve (similar 
to the limiting curve of predissociation at 53800 cm~') can be drawn. Above 
this limiting curve, intensity diminutions are present, but no shifts of rotational level ; 
below the curve no such effects are observed. As the (4,2) band of the B system comes 
to an end witha sudden intensity drop just at this limiting curve and no B-bands with 
v’ >4 were found, the limit at 50500 cm™! is a real predissociation limit. This 
limit corresponds to a continuum, which is in strong interaction with the B?II 
state. At the intersection points above the limiting curve, the eigenfunctions 
of the B?I] state are mixed with the eigenfunctions of the a2 state and this 
mediates the intensity-weakening for the y-band lines. Such a combination 
of predissociation and perturbation effects was called by Coster, Brons and 
van der Ziel (1933) ‘“‘accidental predissociation”’ and was treated theoretically 
by Ittmann (1934) in the case of the Second Positive bands of the N, molecule; 
but neither the perturbing states, nor the cause of predissociation were known 
in the N, spectrum, while in NO both are well established. (See also Geré 
1935 b, and Geré and Schmid 1940.) 

The energy difference of the two predissociation limits found is about 3300 cm7?. 
Now in the energy scheme of N + O atomic term combinations, only one difference 
comes near to this value, namely Be 
the difference N(?D)+ O(P)— 
N(4S)+ OCD) =3355 cm; all 
other differences are much larger. 
Relating these two combinations 
to the limiting curves, the term 
scheme of figure 3 is obtained. 
According to this scheme, the 
dissociation energy of the NO 
molecule is 34 600 cm~ | = 4-29ev. 
The x?II ground state dissociates 
poreinto” the ~ N(*S)-OGP) 4 
ground atomic term combina- 
tion, but into the N?D)+O@P)  ,, 
combination, with a dissociation 
energy of 53800 cm™'. (The  ,, 
extrapolated value, which in 
most cases is 10-25% greater 2 
than the real one, is w,?/(4w,%,) 
=62500 cm~', and, therefore, a 
dissociation energy of 34600 
cm~! for the ground state would 
be too small.) Figure 3. 
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The above dissociation scheme of NO is strongly supported by the photo- 
‘chemical experiments of Flory and Johnston (1935). They irradiated NO with 
the 1832 a. line of mercury and found considerable photochemical decomposition. 
The energy of this mercury liné, 54567 cm, is slightly greater than the dissocia- 
tion limit of the NO ground state at 53800 cm~'; hence the efficiency of this 
wavelength. 
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Collision Broadening of the Inversion Spectrum of Ammonia: 
II. The Collision Cross-sections for Self-broadening 
and for Mixtures with Non-polar Gases 


By B. BLEANEY anp R. P.. PENROSE 
Clarendon Laboratory, Oxford 


ABSTRACT. The widths of the lines of the inversion spectrum of pure ammonia near 
1 cm™ previously determined experimentally (Bleaney and Penrose 1947) are shown to vary 
from line to line in a manner consistent with the assumption of a s:mple dipole-dipole inter- 
action for the collision mechanism. Measurements of the line breadth constant for the line 
(J,K)=(3,3) in mixtures of ammonia with six different non-polar gases are described, and 
from them the collision diameters are calculated. They are found to be approximately the 


same as the diameters obtained from the kinetic theory, whereas in pure ammonia the 
diameters are two to four times greater. 


A survey of microwave spectra indicates that the collision diameter is only significantly 
greater than the kinetic theory diameter for encounters between permanent electric dipoles. 


§1. INTRODUCTION 


HE analysis of the ammonia spectrum has been described in an earlier 
| paper (Bleaney and Penrose 1947, hereafter called ““A”). Of twenty- 
nine lines which were identified, seventeen were so well resolved at a 
pressure of 0-5 mm. Hg that they could be examined in detail. By measuring 
the absorption coefficient at a number-of frequencies the shape of each line could 
be determined and the line breadth constant calculated. The results of these 
measurements were given in table 4 of paper A; for convenience, this table is 
repeated here (table 1), with the order rearranged so that the lines are given in 
order of increasing breadth. 

It will be seen that the values of the line breadth constant at the pressure 
of 0-5mm.Hg vary from 1:8x10-4cem=! to 4:7x 10-4cm—. From them the 
appropriate collision frequencies can be calculated since the collision frequency f 
is associated with a line breadth constant Ay=f/2mc. This collision frequency, 
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which is the number of times per second that the absorption of radiation by one 
molecule is interrupted by the approach of another molecule, varies from 
3:4 x 10" to 8-9 x 10? per second. The collision frequency at 0-5 mm. pressure 
deduced from viscosity measurements by means of the kinetic theory is only 
4-0 x 10® per second, which is considerably smaller. It follows that the inter- 
ruption of absorption of centimetre-wave radiation occurs when another molecule 
approaches to a distance which is considerably greater than that required for an 
appreciable transfer of momentum from one molecule to the other. The relative 
distances can be found by calculating the collision diameters () of the ammonia 
molecule from the formula f= /2vz02N, where N is the number of molecules 
per cm® and v is the average velocity of the molecules. The values of o range 
from 8-7. for the narrowest lines to 14-2 a. for the broadest, whereas the kinetic 
theory diameter is only 4-4. 


Table 1 
Av at 
Line (J,K) p=0-'5 mm. Hg ¥o (db/km.) Lexp db/mol Tineor db/mol 
SC O-= cane’ 
Syl 1:8 7 3:3 3°5 
3,1 2-4 29 14-5 14 
peel 2:6 51 25 Di 
S22 2°6 20 i183 15 
11,9 2-9 8:5 5:7 6-4 
6,3 3-1 44 3595) 34 
See 3:2 100 62:5 59 
Des 3:3 96 71 735 
6,4 36 42 35 34°5 
US 3-9 113) 100 95°5 
Wes 4-1 28 27 7 
8,7 4-1 44 34 32°5 
10,9 4-2 35 25°5 23 
358) 4-5 360 289 289 
44 4-5 190 149 151 
35 4-7 186 148 137] 
6,6 4-7 276 198 210 


Av is measured value of the line breadth constant (em). 
Yo is measured value of absorption at centre of line (db/km.). 


n 
T is integrated intensity = | dv reduced to a length of path containing one mol of gas per 
unit cross-section. °° ” 

Texp is calculated from the relation Ipxp=¥o?7ypAx(10-°Vo), where Vp>=molar volume at a 


pressure of 0:5 mm. Hg. 
oe 4m NyKVo| HrK |? 
Ineor is calculated from the relation Itheor= sea ee 
Nyx=number of ammonia molecules per cm* in rotational level (JK). 
lever | 2—-square of dipole moment associated with the line (J,K). 


These large values for the collision cross-section may be attributed to the effect 
of the large electric field of the ammonia dipole, which will produce a considerable 
perturbation of the energy levels of an approaching molecule at distance of the 
order of 10-7cm. ‘This assumption is supported by consideration of the large 
variation in the widths of the different lines. This is at first sight surprising, 
since each line arises from an exactly similar transition between the levels of the 
inversion doublet. The significant difference arises, however, from the fact 
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that each line is due to molecules in a particular rotational energy level, and the 
line widths vary with the state of rotation of the molecule in the same way as does 
the Stark splitting due to an external electric field. In the following section 
of this paper the variation of line width is discussed, and it is shown that it can 
be explained by the simple hypothesis that a collision occurs when the electrical 
forces between the molecules due to their strong dipoles cause the energy of 
interaction to reach a certain value. The magnitude of the interaction energy 
required to produce a collision is evaluated from the observed line widths by 
averaging over the various settings of the dipoles with respect to the line of centres. 
at a collision. 


§2. THE COLLISION CROSS-SECTION FOR SELF-BROADENING 


The magnitude of the splitting or perturbation of the energy levels of a 
molecule possessing a permanent electrical dipole moment when subjected 
to a constant electric field is determined by the steady component of the dipole 
moment. In the general state of rotation of a symmetrical top molecule, the 
dipole moment will be precessing about the axis of total angular momentum 
of the molecule, and the steady component is determined by the projection of 
the moment on that axis. The frequency of precession is of the same order as. 
the rotational frequency, and one may therefore expect that the perturbation 
due to a transient electric field is still determined by the steady component of 
the dipole moment if the duration of the transient is long compared to the period 
of rotation. Since the experimental collision diameter is about 10-7cm. for 
ammonia, the duration of a collision is of the order of 10-1 sec., while the rotational 
periods are less than 10-sec. ‘Thus it is reasonable to assume that the inter- 
action between the electric fiekd of one molecule and the dipole moment of an 
absorbing molecule is determined by the steady component of the latter. In 
ammonia this varies as K/(J?+J)?, where J,K are the rotational quantum 
numbers of a symmetrical top. The interaction should therefore be largest 
when K=, and smallest when K is small compared with J. Inspection of 
table 1 shows that the line widths increase steadily as K/J approaches unity. 

The variation of the line breadth constant with the rotational state of the 
molecule may be given a quantitative aspect by assuming that a collision occurs. 
when two molecules approach to a distance such that the energy of interaction 
between them reaches a certain value W. For the molecules in a particular 
rotational state defined by the quantum numbers (J, K), the average distance 


at which a collision takes place will then be determined by the relation Wa p/r3. 
The variation in the line breadth 


constants depends only on the 
variation in the collision cross- 
section. Since the latter is pro- 
portional to 7?, the line breadths 
should vary, on this hypothesis, 
as (u/W)*. Hence, assuming 
that W is independent of the 
rotational state of the molecule, 
Av should vary as (K2/(J?+.J))$ 
for the lines of this inversion 
spectrum. In figure 1 the widths 


Figure 1. Widths of lines in ammonia inversion. 
spectrum near 1 cm7!. 
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(Av) of the seventeen lines measured at 0-5 mm. Hg pressure are plotted against 
(K?/(J? +J))*. Within the experimental error (about 5° in Av) the points fall 
on a straight line, thus confirming that the interruption of the absorption is due 
to dipole-dipole interaction. 

The success of this simple theory encourages a more exact calculation of the 
magnitude of the interaction energy between two ammonia molecules required 
to produce a collision. ‘This requires an averaging process in order to take 
account of the variation of the interaction energy between two dipoles with the 
angles at which they lie with respect to the line through their centres. A simple 
method of obtaining the average cross-section for collision is given in the Appendix 
to this paper. It is found that the average cross-sectional area for collisions 
of this type is 

Sie ee eee (1) 
where pu,? = .?{K2/(J? +.J)} is the squared dipole moment of an ammonia molecule 
in the rotational state (J, K), 13 =the mean square of the dipole moment of all the 
ammonia molecules, and W=interaction energy required to produce a collision. 

Hence the line breadth constant Av is 


DS ae Ii rwe 
Nps oe JF) ee ee ee ravers terse (2) 


where N is number of molecules per cm’, M is molecular weight, T absolute: 
temperature, and c and R have their usual meanings. 

At room temperature, computation of the average value of {K?/(J?+.J/)} 
over all rotational states of ammonia gives 13=0-40yu. Inserting the expression 
for S,,, and taking w=1-44 x 10°18E.s.U., we obtain 

N= AS x Alt pee a( Kei) =) ))ee (cm). | Sane aee (3a) 
where Pym is the pressure of ammonia in mm. Hg. 

The experimental value of Av at a pressure of 0-5 mm. Hg, obtained from 
figure 1, is 

Ai = 0 50ocl0e (Kea (J2-es))> (cmt), 9 oienaee (3d) 
Hence the interaction energy W is 3:0 x 107! ergs, or W/he=1-5 cm“, which. 
is approximately equal to twice the separation of the two energy levels of the 
inversion doublet, between which transitions are taking place. It is thus hardly 
surprising that the absorption or emission of radiation is interrupted. 

The mechanism of this collision process is rather different from that assumed 
in the simple theory of collision broadening of spectral lines in the optical region. 
For such lines the periodic time (10-!°sec.) of the radiation is much shorter 
than the duration of a collision between two molecules (10- to 10-™ sec.). The 
‘classical picture”’ of a collision is that of a small change in the natural frequency 
of the molecule due to the influence of an approaching molecule, which causes. 
the oscillating molecule gradually to get out of phase with the exciting radiation. 
This effect persists over some hundreds of oscillations throughout the duration 
of a collision, and the integrated effect causes such a large phase shift that the 
absorption or emission of radiation is interrupted. The perturbation of the 
energy levels caused by the approaching molecule is, however, small compared 
with the separation of the two levels between which transitions are being induced. 
This is very different from the situation in the ammonia inversion spectrum. 
Here a collision is caused by a large perturbation of the energy levels acting 
for less than one period of the radiation. 
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§3. THE CROSS-SECTION FOR COLLISIONS WITH 
NON-POLAR MOLECULES 

Since the abnormally large collision diameters in pure ammonia are ascribed 
to the effects of dipole-dipole interaction, it would be expected that the cross- 
section for collisions between an absorbing ammonia molecule and a non-polar 
molecule would be considerably smaller. Measurements were therefore under- 
taken of the line breadth constant of the strongest line in the inversion spectrum 
of ammonia (J,K =3,3) for admixtures of six non-polar gases. Since none 
of the gases used had appreciable absorption at centimetre wavelengths, the 
effect of diluting the ammonia with one of them is to increase the collision fre- 
quency for the same partial pressure of ammonia and thus to broaden the line. 
Although it was possible to determine the line breadth constant for the mixture 
by the method already described (Bleaney and Penrose 1947) for the inversion 
lines of pure ammonia, the close agreement with the theoretical intensity obtained 
in the previous measurements on this line made it possible to use a different 
technique which was not only simpler but rather more accurate. Instead of 
measuring the absorption coefficient at a number of frequencies in order to 
delineate the line, the absorption coefficient was measured only at the centre 
of the line. For a mixture with a non-absorbing gas the absorption there is 
smaller than in pure ammonia, varying inversely as the line breadth constant 
at pressures where the absorption due to tails of neighbouring lines is negligible. 
The formula for the absorption coefficient at the centre of a line is 

vol Sirs | LK | vo N. JE 

i 3RT Av 
where | 7x |?, Nyx and Av are as already defined and v is the wave number of the 
line in cm-1. 

The line breadth constant Av=f/27e, where f is the collision frequency. 
For a given mixture both fand N,, vary directly as the pressure, and the absorp- 
tion coefficient at the centre of a line is therefore independent of the pressure. 
The addition of a non-absorbing gas to a certain concentration of ammonia 
increases f from the value for pure ammonia, f,;, to f,; +f,2, where fi, is the fre- 
quency of collisions between an ammonia molecule in the state (3,3) and molecules 
of the foreign gas. The absorption is therefore reduced in the ratio 


9/04 =fir/(fr +h) 
where «, is the absorption coefficient for pure ammonia and «, that for the mixture. 
Hence fyo/fir = (%/%2) —1. 
If Sy. and Sj, are respectively the cross-sections for ammonia-foreign gas | 
collisions and for ammonia~ammonia collisions, the expressions for f,, and for 
fra are 


per cm. Per G)) 


hie =2N.S}(2RT/7)*(My* + M3"), 
fir =4N,S1,(RT/7M,)', | 
where N,, N, are the numbers of ammonia and foreign gas molecules per cm? | 
respectively, and M,, M, the molecular weights of ammonia and of the foreign | | 


gas. 
Hence the ratio of the cross-sections is 


Sie N, 12 Lif Sx 
so Gd MM = (2-1) GC eM...) 
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§4. EXPERIMENTAL 


The six non-polar gases used in the mixtures were helium, hydrogen, nitrogen, 
oxygen, argon and carbon disulphide. The ratios in which the mixtures were 
made up were chosen so that the line breadth constant for the mixture would be 
approximately double that for the same partial pressure of pure ammonia. Two 
slightly different techniques were used to prepare the mixtures, the choice 
depending on whether the non-polar gas would condense in liquid air or not. 

The mixtures were prepared in a bulb of about one litre volume which had 
a small side tube which could be immersed in liquid air. In cases where the 
foreign gas would not condense at liquid air temperature, the bulb was first 
filled with ammonia at a known pressure. After the bulb has been isolated, 
the ammonia was frozen out in the side tube (the vapour pressure of ammonia 
at 90°K. is about 10-*’mm.Hg). Then a suitable pressure of the foreign gas 
was admitted to the bulb, after which the bulb was again isolated. Finally 
the liquid air bath was removed, and the ammonia in the side tube allowed to 
evaporate and mix with the other gas, several days usually being allowed to 
elapse to ensure proper mixing. Since the volume of the side tube was negligible 
in comparison with that of the bulb, the error due to its different temperature 
was negligible. 

When the foreign gas would condense at liquid air temperature, the following 
procedure was adopted. A second bulb was filled to a known pressure with the 
foreign gas, which was then transferred to the first bulb by freezing it into the 
side tube by means of a liquid air bath. ‘This operation was repeated with 
ammonia. ‘The first bulb was then isolated, and the gases allowed to evaporate 
and mix as before. 

In each of these methods a given volume was filled to a known pressure 
with each of the constituents of the mixture. ‘The ratio (NV,/N,) of the numbers 
of molecules per cm? of the two gases was therefore accurately known. 

The intensity of absorption of the mixture of gases was determined by 
observing the reduction in power transmitted through a cavity resonator of 
known Q, when the gas was admitted. A full account of this technique has 
already been given in paper A. In terms of the detector readings 5) and 6, 
before and after admitting the 
mixture, the absorption  co- 
efficient is given by the expression 
a= (27/AQo)((So/81)! 1). The 
absorption was measured for a 
number of pressures such that 
the partial pressure of ammonia 
lay in the range 0:2 mm. to3 mm. 
Figure 2 shows the absorption— 
pressure curve for an ammonia-— 

- oxygen mixture in the proportion 
N,/N,=0:107. The falling-oft 
in the absorption at low pressures 
is due to the disturbance of 
thermal equilibrium (Bleaney and 0 Pieeciet 
Penrose 1948), while the rise at Figure 2. Absorption in ammonia-oxygen mixture 
high pressures is attributable to at 0:7964 cm=2. 
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the contributions from other lines in the spectrum. The shape of the full curve 
includes an allowance for this rise, and provides a means for smoothing the 
points. The dotted line represents the corrected value (179 db/km.) for the 
line (3,3) in this mixture. The absorption coefficient for the centre of this line 
in pure ammonia is 360 db/km. (see A, or table 1 of this paper). From the ratio 
of these two absorptions, the ratio of the collision cross-sections for ammonia- 
oxygen and ammonia—ammonia can be calculated by means of equation (5). 
Since the latter is known accurately from previous measurements by the authors, 
the collision diameter (,,) for ammonia—oxygen can be evaluated from the 
formula 045 = ~/(S}_/7). 

From similar measurements and calculations for the other mixtures, the 
value of the collision diameter in each case was determined. 


§5. RESULTS AND DISCUSSION 


The results of these measurements are shown in table 2. 


Table 

(1) (2) (3) (4) (5) (6) (7) (8) - 
Helium 4-0 0-21 0:0295 17-6 2-35 5220 3-6 
Hydrogen 2-0 0:78 0-064 38-2 3-505 0325 Same 
Nitrogen 28-0 1:72 0-205 122 6-4 4:09 5:0 
Oxygen 32-0 iholsyil 0-123 WES 4°85 402 4-9 
Argon 39-9 1-74 O-tat 66:3 4-6 4:04 4-0 
Carbon disulphide 76:1 8-6 0-296 177 TS —- 6-9 
Ammonia 17-0 — — 597 13-8 4-4 — 


(1) Gas mixed with NH; ; (2)molecular weight ; -(3) polarizability (in cm*® x 10-*4) ; 
(4) S'42/S1, measured ; (5) S;, measured (in cm? x 1071) ; (6), (7), (8) collision diameters 
in A. ; (6) measured; (7) according to kinetic theory ; (8) calculated. 


In the first three columns are shown the non-polar gas, its molecular weight 
and its polarizability. ‘The experimental value of the ratio of the collision 
cross-section to that for the line (3,3) in pure ammonia is given in column (4), 
and the actual cross-section in column (5). The last three columns contain res- 
pectively the experimental value of the collision diameter, the kinetic theory 
molecular diameter and the diameter calculated by an extension of the theory 
of §2 of this paper. The values of the kinetic theory diameter are not, of course, 
measured directly for each mixture; the figure given is obtained by taking the 
sum of the radii of ammonia and of the admixed gas (Kennard 1938). 

The calculated diameter given in the last column is obtained as follows. 
A molecule of a non-polar gas becomes polarized in the electric field of an ammonia 
molecule, and its induced moment produces a field which acts on the permanent 
dipole of the ammonia molecule. It is easy to show that the energy of the latter 
interaction is r~a| wx |?(1 +3 cos?6), where a= polarizability of non-polar mole- 
cule; r= distance between centres of ammonia molecule and non-polar molecule; 
#=angle between dipole moment of ammonia molecule and radius vector. 

If the assumption is made, as before, that a collision occurs when the energy 
of interaction reaches a certain value W, then by an averaging process similar to 
that given in the Appendix, one obtains for the collision cross-section 


Sy =3-90(u2a/W)*, 


Collision broadening of the inversion spectrum of ammonia: III 547 


If the value of W necessary for a collision is assumed to be the same as that found 
for an ammonia—ammonia collision, i.e. 3-0 x 10-16 ergs, the value of S,,, and 
hence the collision diameter, may be calculated. 

Inspection of table 2 shows that the experimental values of the collision 
diameter are of the same order as the kinetic theory values, as, also, are the values 
calculated on the basis of the theory outlined in the preceding paragraph. The 
latter means that at the distance between the molecules required by the purely 
electrical interaction for a collision, other forces come into play, and that the 
collision mechanism suggested is superseded by more powerful interactions. 
It is satisfactory, however, that the theory does not give values of the collision 
diameter which are much bigger than those observed. 

The collision cross-section for ammonia-—nitrogen appears to be exceptionally 
large in comparison with those for ammonia—oxygen and ammonia~argon. 
This was unexpected, since the polarizability and the kinetic theory diameter 
of each of these molecules is very nearly the same. The first measurements 
were made with commercial nitrogen, and gave a collision cross-section of 
128 x 10-®cm?. ‘The experiment was then repeated with a mixture made from 
specially pure and dry nitrogen, for which the value shown in table 2 was obtained. 
The two values agree, within the experimental error, and the abnormally high 
value seems, therefore, to be correct. ‘There appears, however, to be no obvious 
explanation. 

In the rotational state (3,3) the dipole moment precesses about the axis of 
total angular momentum at a fairly small angle, and its projection on that axis 
islarge. Its effect in producing interactions that may cause collisions is therefore 
almost a maximum, corresponding to the fact that the line (3,3) in pure ammonia 
is one of the broadest in the inversion spectrum. ‘The collision diameter is 
14-2., whereas the line (5,1) gives a collision diameter of only 8-74. in pure 
ammonia. ‘This wide variation between the different lines seems unlikely to 
occur for collisions between ammonia and non-polar molecules, since the collision 
diameters are approximately the same as those given by the kinetic theory. 
For the same reason it is unlikely that there will be any significant difference 
between collisions which broaden the line and those which produce thermal 
relaxation (Bleaney and Penrose 1948), since a kinetic theory collision must 
surely be effective in producing thermal relaxation. 


§ 6. COMPARISON WITH OTHER CENTIMETRE WAVE SPECTRA 


There are as yet few data on the widths of lines in centimetre wave spectra 
of molecules other than ammonia. Accurate measurements have been made 
on the water-vapour line at 0-742 cm™! for mixtures of water vapour and air by 
Autler and Becker (1946), and for pure water vapour by ‘Townes and Merritt 
(1946). Atatemperature of 45°c. the former workers find Av = 0-087 + 0-001em™ 
for zero concentration of water vapour in the atmosphere, rising to 0-107 + 
0-001 cm-! at a concentration of 50gm/m°. From these figures the value of 
the collision diameter for water-air is 5:-4.4., which lies between the values for 
ammonia—oxygen and ammonia-nitrogen given in table 2. The diameter for 
water—water collisions is found to be 10-2.., while the value calculated from a 
direct measurement by Townes and Merritt on pure water vapour at a pressure 
of 0:103 mm.Hg is 10-0a. These values are comparable with those for the 
narrower lines in pure ammonia; this is satisfactory, since the water molecule 
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is an asymmetrical top, and the dipole moment will have a fairly small steady 
component in the two rotational levels between which the centimetre wave 
transition takes place. 

The measurements of Beringer (1946) on the oxygen spectrum at 2cm™ 
show that the line breadth constant at atmospheric pressure must lie between 
0-02 and 0:05cm-!. ‘This gives for the collision diameter 2:8 to 4-5 a., whereas 
the kinetic theory diameter is 3-64. The collision diameter for oxygen—nitrogen 
was the same. 

The main conclusion to be drawn from these results is that the collision 
diameter is only significantly greater than the kinetic'theory diameter for encounters 
between permanent dipoles. For collisions between permanent dipoles and 
non-polar molecules the diameter is approximately the same as the kinetic theory 
diameter, and the measurements on the oxygen spectrum suggest that this is 
also true (as would be expected) for collisions between two non-electrically 
polar molecules. 


APPENDIX 
Calculation of the collision cross-section 


Let p,=pu(K?/(J?+J))? be the mean dipole moment averaged over the 
rotation of the molecule absorbing radiation, and let p, similarly be the mean 
dipole moment of an approaching molecule which causes an interruption of the 
absorption. Then the latter produces at a point whose coordinates are (r, @) 
with respect to the direction of its mean dipole moment an electric field E whose 
magnitude is 

| Z| =per(1 +3 cos? 6). 
If the dipole », makes an angle ¢ with this field, the interaction energy is 


LE cos ¢ = py per? cos (1 +3 cos? 6). 


> 


The assumption is made that a “collision”? occurs when this energy reaches 
a certain value W. ‘Then the dipole x, can be imagined as surrounded by a 
surface at which the interaction with the dipole ., reaches the required value 
for a collision. Owing to the directional nature of the dipole this surface is not 
spherical, but has the form of a prolate spheroid whose major axis is twice as 
long as the two minor axes. 'To find the mean cross-section of the molecule 
it is sufficiently accurate to find the volume of this spheroid, and to use the cross- 
section of a sphere of equal volume. Since the interaction energy varies as 
1/r?, this leads to a very simple mathematical procedure. 

The equation for the surface of the spheroid is 


7? = [Ly fy cos f(1 +3 cos? 6)?/W. 


Hence its volume V is 
| 2 74/2 
aa Be ode - Hts SOS |" (+3 cos? 6)! sin a0 = e 1-38, Haba cose 
~ 0 


since the value of the integral is 1 +(In(2 + +/3))/2./3 =1:38. 

This expression for the effective volume of an approaching molecule requires 
to be averaged over the mean dipole moments of all approaching molecules, |} 
and over the angle ¢ which the absorbing dipole 1, makes with the field of the || 
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approaching dipole. ‘The former averaging requires simply that ft be replaced 
by j4z, the average of the mean dipole moments of all the molecules present ; 
the second averaging gives simply a factor of one-half. Hence 

a 2a 
rat 


Equating this volume to that of the equivalent sphere, one finds for the mean 
cross-section area of the approaching molecules 


Sy. =7(3V/47)8 = 7(0-69 14 129/ W)* = 2-45 (1,208) W2). 


V 1-38 puype9/W. 
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ABSTRACT. A review is made of the derivation of formulae for the efficiency of counting 
systems, comprising in the general case a counter and amplifier, a scaling circuit and a 
recorder. In the case of recorders of constant resolving time no complete formula has 
been found in the literature, and this has been derived. A number of existing formulae 
are shown to be valid only fora particular type of variable recorder. A discussion is given 
of the use of the formulae in practical counting experiments. 


§ 1 
N various branches of nuclear and atomic physics it is necessary to measure the 
intensity of some kind of incoming radiation, the measurement taking the form 
of a count of the number of particles or photons arriving within a given time. 
‘For this some form of counting system is employed. In nuclear physics a counting 
system normally consists of a counter and amplifier feeding a scaling circuit 
; (hereafter termed a “‘scaler’’) which in turn actuatesarecorder. ‘he counter may 
be an ionization chamber, Geiger-Miiller tube etc., depending upon the type of 


PROC. PHYS. SOC. LX, 6 36 


550 M. Blackman and 7. L. Michiels 


radiation to be detected. In all cases its function is to produce an electrical 
impulse upon the passage through it of a particle or photon of the incoming radi- | 
ation. ‘These impulses are then fed into the remainder of the system. In sucha | 
counting system the various components have finite resolving times. ‘That is, 

after they have been actuated, a certain recovery time must elapse before they can — 
again be actuated. Because of this, a certain fraction of the incoming radiation _ 
which arrives during the recovery time will not be detected. A large number of 
theoretical and semi-theoretical discussions occur in the literature in which the | 
efficiency of counting systems is treated. In the simplest of all cases, in which the | 
counting system consists of a counter, amplifier and a simple recorder, it is usual to | 
assume an infinitely small resolving time for all parts of the system except the | 
recorder which has a constant resolving time 7. The efficiency £ of this system is 

then defined as the ratio of the counts recorded in a given time to the number of © 
particles arriving in that time. It is rather astonishing to find, even in this simple | 
case, that there are at least three different expressions for this efficiency E to be | 
found in the literature, namely :— 


(a) E=1/(1+pr) (Ruark and Brammer 1937), —...... (1) 
(b) E=exp(—pr) (Volz 1935, Schiff 1936),  —..... (2) WW 
(c) E={1—exp(—pr)}/ur (Locher 1933). a. (3) 


In these formulae, p is the average number of particles arriving per second, the. 
distribution of particles being supposed to be a random one. ‘There are also 
generalizations of (a) and (6) due to Ruark and Brammer (1937) and Alaogiu and 
Smith (1938) for cases where the finite resolving time of the counter is taken into 
account and where a scaler is introduced. In view of the different versions which 
are current, it seems worth while to examine the problem afresh in order to find 
which formula is correct, and how it should be generalized. * 


§2 

In discussing the efficiency of counting systems, it is essential to adopt a clear 
terminology for the impulses in various parts of the systems. In this paper, the 
incoming radiation will be said to consist of a stream of ‘ particles’. The electrical 
impulses produced by the counter will be referred to as ‘ pulses’ and an event which 
actuates the recorder will be termed a ‘count’. ‘The efficiency E of the counting |} 
system is defined as the ratio of the number of particles counted in a given interval | 
of time to the number of particles arriving at the counter within thistime. In this 
connection it should be remarked that it is possible for some of the incoming | | 
particles not to actuate the counter for reasons other than that it is recovering from | 
the effect of a previcus particle. . In a Geiger-Miiller tube, irradiated by X rays, | 
for example, only a certain fraction of the incoming photons will eject an electron) 
from the wall of the counter and so be capable of producing a discharge within it. }} 
This cause of inefficiency is not dealt with in these calculations which are concerned 
solely with the effect of finite resolving times within the system. The incominp| | 
particles are assumed to be distributed at random and arriving at an average rate bb | 


§3 
Consider first a system in which no scaling circuit is employed. The recorde 
is assumed to have a constant resolving time 7 and the counter and amplifier to have: 


* Results essentially in agreement with those of the present paper have recently been derived b H 


Jost (1947) by a very elegant method. Jost’s paper, however, emphasizes the mathematical rather} 
than the physical aspect of the problem. 
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a resolving time negligibly small compared with this. Suppose that N counts 
occur in a certain experiment of duration 7. After each count the recorder is 
inactive for a period 7, Let p be the average number of particles which occur 
within a period 7 after a given particle. ‘Then, on average, this number of particles 
is lost at each count. Hence the total number of particles arriving during the 
experiment is N+pN and the efficiency of the system is 


E=N/((N+pN)=1/(1+p). we... (4) 


It is now necessary to calculate p. It is well known that for particles distri- 
buted at random, the probability W(m, t) of m particles arriving in an interval ¢ is 


W(t. 0) ee iar) 7 ee en 8 eek (5) 


Consider a time interval t commencing immediately after the arrival of a particle. 
The above formula is still valid since there is no correlation between the time of 
arrival of particles. Hence the probability of m further particles following the 
first particle within the resolving time 7 is W(m,7). In such a case the number of 
particles lost is m. Hence p the average number lost in time 7 is 


p= X& mW(m,7) 


m=1 


— eH, [Lt SS (7) /(m = 1) ! SSH, (6) 
1 


The efficiency is therefore given by 
En) sige a eae ee (7) 


This formula is one of three quoted in the introduction and has been derived 
in various ways in the literature, for example by Ruark and Brammer (1937). 
‘The formula (2) has been derived by Schiff (1936) who employs the following 
argument for a recorder with a constant resolving time. ‘The probability of no 
particle arriving within a time 7 after a given particle is e“” (see equation (5) for _ 
m=(). This is then put equal to the probability of counting the next particle 
which in turn is put equal to the efficiency. ‘This expression gives the probability 
that no loss will occur, but it is not immediately related to the magnitude of the 
probable loss which does occur. At very small counting rates, the number of 
particles which arrive in an interval 7 will be either 0 or 1 to a very good approxi- 
mation, and in this limit the efficiency (1 —yr)is correct. Inall other cases it bears 
no relation to the true efficiency. It will be seen later, however, that this formula is 
correct for another type of recorder in which the resolving time is not constant. 
Formula (3) has been given by Locher (1933). He finds 


Pali eaece yf yt {PL Biss. (8) 
He takes a series of intervals of duration 7 at random, and assumes that if there 
is only one particle within such an interval it is counted. If there are two, one is 
lost and so on. With this assumption 
na 
p= &X (m—1)W(m,7), 

m=2 

_ which can readily be seen to lead to (8) if r is identified with the resolving time of the 


recorder. He then considers an “‘ experiment” of duration 7 in which the average 
number of particles arriving is 7 and hence deduces the average fraction lost and 
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the efficiency, obtaining equation (3). This result differs from the correct result 
given by (7), particularly for very low counting rates in which region the loss given 
by (3) is half that given by the correct formula. The error in Locher’s argument 
is that he has ignored the fact that the intervals 7, having been defined as times. 
during which the recorder is inactive, start only after a particle has been recorded. 
Consequently if, as in Locher’s calculation, such periods are chosen at random, the 
first particle, which actuates the recorder, will not always occur at the beginning of | 
the period and hence the effective resolving time assumed in the calculation is in 
general less than 7. Evidently, if the counting rate is very low, the probability of | 
the first particle occurring anywhere within a period 7 will be constant and the 
effective resolving time assumed will be 7/2. If the counting rate is very high, the 
probability of the first particle occurring very near the beginning of the period will 
increase and Locher’s formula should then tend to the correct one, which in fact it 
does. | 
§ 4 
Consider now a counting system in which a scaler of scale-factor m is used. | 
The scaler and all components of the system previous to the recorder are again 
assumed to have resolving times negligibly small compared with that of the recorder 
which has the constant valuer. As before, an experiment is considered in which V. 
counts are obtained in an interval T. Let p again be the average number of | 
particles lost during an interval 7 after acount. The number of particles recorded | 
is nN and the total number arriving ismN+pN. The efficiency is therefore given 
by | 


E=nN|(nN+pN)=n[(n+p)=1f(l+pjn). —..... (9) 


In this procedure the number of pulses(<x) left in the scaler at the end of the 
experiment has beenignored. ‘This is a sufficiently good approximation when Nis. 
very large compared with. In this connection it should be noted that the formula 
derived for the efficiency will in any case be valid only if the number of counts is 
large. ‘This follows from the fact that the number of particles lost after a count is. | 
obtained by taking a statistical average over the possible cases which can occur ina | 
large number of counts. If only a small number of counts are taken, the loss will | 
show fluctuations about its mean value p which may be quite large. In experi- ) 
mental work, it is usual to count a number of particles sufficiently large for the | 
statistical fluctuation in this number to be small. If, however, a scaler of large i 
scale-factor is employed, the number of counts may not be sufficiently large for the | 
fluctuation in p, the average number of pulses lost per count to be-small. In| 
applying the expressions for the efficiency of the system, it is therefore necessary | 
to consider both the value of p and its fluctuations. | 

It is now necessary to calculate p in this case. Suppose that m particles arrive | 
at the system within a time 7 after a count. If m is less than n, these are stored in | 
the scaler and none are lost. If mis equal to or greater than n but less than 2m, then |} 
n are lost. In general, if gn>m=>(q—1)n, then (q—1)n are lost. Hence the| 
expression for p is 


p=n{W(n, 7) + W(nt+ 7) tose + W(2n —1,7)} 
+ 2n{W(2n, 7) + W(2n+1,7)+ .... +W(3n-1,7)} 


I 
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A function L(n,7) can be defined by 
L(n, 7) ={W(n,7)+W(n+1,7)+....to } 
ACES A (Un yC Haran! ise acon (eae OMe ee ar en ee eee (10) 
In terms of this function 
p=n{L(n, 7) —L(2n, r)} + 2n{L(2n, 7) —L(3n,7)} 4+ .... 
Saft Lehtin gr pleomery oat Qh a (11) 


The function L(v,r) can be expressed in terms of tabulated functions (Molina 
HAZ) 

It is evident that if r/n is small compared with unity, the series L(n,7) will 
converge very rapidly. In the case where piz/m is very small, an approximate 
expression can be obtained from (11), 


pene" {(yr)"n}} 
~and E=1—e'{(pr)"/n!} TG TUN soe ad (Re ei Moceicle (12) 


‘Vhe expression (11) does not, as far as we know, appear to have been given in 
the literature. Lewis (1937) has given an expression for the loss in the above type 
of system in the region of small losses. His formula gives the fractional loss in the 
form of an expression which is equivalent to the function L(n, 7) instead of p/(p +n). 
This is a valid approximation. A somewhat similar expression has been given by 
Alaoglu and Smith (1938) who introduced the function L(v,7). They use it in 
conjunction with a different formation of the efficiency of the system. ‘They 
calculate the probability that the time ¢ in which 4 group of n particles follow a 
given particle (that is, after a count) islessthanz7. This probability they find to be 
the function L(n,r). The probability G(n,7)=1—L({(n,7r) of the time ¢ being 
greater than 7, they call the efficiency of the system. ‘This is not a satisfactory 
definition of efficiency for a recorder of constant resolving time, but the formal 
result applies to a certain type of recorder of variable resolving time. ‘This type 
of recorder is discussed in § 6. 

§5 

The system considered can now be further generalized to include the case in 
which the counter and amplifier have a constant resolving time o. ‘The recorder 
has a constant resolving time 7. It is assumed that the resolving time of the first 
stage of the scaler, if not negligible, can be incorporated in that of the counter and 
amplifier and that those of the succeeding stages of the scaler can be neglected. In 
such a system, not all the particles arriving at the counter will produce pulses which 
are fed into the scaler and thence to the recorder. Owing to the effect of the re- 
covery time of the counter, the pulses which reach the recorder are no longer 
randomly distributed in time. Alaoglu and Smith (1938) have introduced a 
probability function L*(m, t, 0) giving the probability that the time of arrival of a 
group of m pulses, after a given pulse, is less than a specified time t. ‘This function 
can be calculated in the following way. [tis first necessary to obtain the probability 
that a pulse will occur at a time between f and ¢ + d¢ after a given pulse, there being 
no intervening pulses. The probability that no particle will arrive within an 
interval ¢ is, for a random distribution, equal toe“. Now the probability that 
no pulse will occur in a period ¢ following a pulse is the same as the probability 
that no particle will arrive in the interval between o and ¢, since no pulse can occur 
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between 0 and, and for any time after o an incoming particle will produce a pulse. 
That is, the probability that no pulse will occur in the interval ¢ is et eo But 
the probability of a particle arriving in the subsequent interval dé is pdt. Hence 
the probability Q*(¢, c) dt of a pulse occurring in the interval between ¢ and t+ dt 
after a given pulse, with no intervening pulse is 


Q*(t,0)dt=e" "ndt, to 
=0 EG; : 


The function L*(1, t,c) which expresses the probability that the time interval 
between a pulse and its successor is less than ¢ is therefore given by 


LA(1,t, a= O*is)d | = eee (14) 


In order to calculate L*(m+1,t,c) from L*(m,t,c), one can proceed as 
follows (Alaoglu and Smith 1938). Consider a pulse occurring at a time between 
x and x +dwx after a given pulse («<?), there being no intervening pulses. The 
probability of obtaining m pulses in a time less than (t—x) is by definition 
L*(m,t—x,c). Hence the probability of a pulse occurring at x followed by m 
pulses within the succeeding interval (¢ —x) is 


Q*(x, 0) dx L*(m, t —x, 0). 
rt— mo 
Therefore L*(m+1,t,0)= | Q*(«, 0) L*(m, t —x, 0) dx. 
0 


The upper limit in the integral is determined by the fact that it must be possible to 
fit at least m pulses of durationo in the interval (¢ —x); alternatively this can be 
stated in the form 


L*(m, t,o) =0 if t<mo. 


Tt has been shown by Alaoglu and Smith, using a proof by induction, that the 
function L* is given by 


L¥(m, t,o) =({u(t —mo), m}/(m —1)! 


=e enn—mor{ = ma)|™ [ult—mo) |e i 
m | (m+1)! dead it 


Here, as before, I'(«, m) is the incomplete gamma function. 

Now L*(m, t, a) can also be considered as the probability of m or more than m 
pulses occurring within the interval ¢ following a pulse. There is, however, an 
upper limit M to the number of pulses which can possibly occur within a given 
interval owing to the finite recovery time o following each pulse. This limit is 
given by t=Mo+f, f<o. When this limit is exceeded, L*(m, t,o) is zero. 

Let W*(m, t, a) be the probability that exactly m pulses occur in an interval f 
following a given pulse. Then L*(m,t,c) can be written as 


L*(m, t, 0) = W*(m, t,o) + W*(m +1, t, o)+....+W*(M, to). 


‘That is, the probability of m or more pulses occurring within an interval t is the 
sum of the probabilities of exactly m, m+1 etc. pulses arriving within the interval ¢. 


Henc 
: W*(m, t,o) =L*(m, t,0) —L*(m+1, 1,0). 
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If mis greater than M, L*(m, t,o), and thus W*(m, t, c) is zero. 
Let p* be the average number of pulses lost by the recorder in an interval 7 
following a count. Then 


p* =n{W*(n,7,0)+W*(n+1,7,0)+.... +W*(2n—-1,7,0)} 

+ 2n{W*(2n, 7,0) + W*(2n+ 1,7, p)+..... +W*(3n—1,7,0)} 

aP oooc 
This can be written 

p* =n{L*(n, 7, 0) —L*(2n, 7, 0)} + 2n{L*(2n, 7, c) —L*(3n,7,6)}4+ . 
that is p* =n{L*(n,7, 0) + L*(2n,7,0)+.... +L*(qn,7,0)}, ...... (17) 
where gun<M <(q+1)n. 
If, as before, N counts are recorded, the number of pulses leaving the amplifier 

isnN+p*N. From (6) the average number of particles arriving during a pulse 


is (1+yo). Hence the number of particles arriving during the experiment is 
N(n+p*)(1+ 0). The efficiency of the counting system is therefore given by 


E=nN[N(n+p*)\(1+yo)=1/(1+p*/n)(ltyo). se... (18) 


This formula is seen to reduce to those given for the simpler cases if first o is put 
equal to zero (cf. (10) and (15)) and then z is put equal to unity. Also if 7 is less 
than no, p* is zero and the only loss is that due to the counter and amplifier. This 
should be so, since in this case the recorder has already recovered before another 
group of pulses can have arrived from the amplifier. If p(t —7o)/n is much 
smaller than one, the expression for the efficiency can be reduced to the approxi- 


mate form 
h=(1 — eM) (7 =o) jn \/\(Te yo)... ae (19) 


In figure (1), the fraction of to ———— 
pulses counted 1/(1+p*/n) has 
been plotted as a function of 
u(t —no)/n for various scale- 09 
factors. In computing values of 
p*, use has been made of tables 
by Molina (1942). If o is zero, 
these curves give the efficiency 
of the system directly. In the 
general case, the value obtained 
from the curves must be divided 
by (1+ 0) to give the efficiency. 

Feather (1943, 1945) has also, 441 
by a different method, derived 
expressions for the statistical 


° 
co 
T 


Ss 
~n 
= 


Fraction of Pulses Recorded 


distribution in time of pulses 05 ? an - ' 
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from a counter, the results being Wain 

generalized to include the case Figure 1. Fraction of pulses recorded as a function 

when the incident stream of of the product of ((r—70)/n), the effective 


articles has a non-random dis- resolving time, and the average rate of entry 
BP of particles into the counting system, for 


tribution. He has applied his various values of n, the scale factor of the 
formulae to the calculation of the system. 
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losses in a system consisting of a counter and recorder, but has not, however, 
extended them to take into account the effect of a scaler. 


§6 

It remains to consider the situation which arises when the resolving time of 
the recorder is not constant. The efficiency must then depend upon the way in 
which the resolving time varies, that is, upon the mechanism of the recorder. One 
case which arises in practice is that in which the arrival of a pulse within the resolv- 
ing time re-excites the recorder without, however, recording that pulse. ‘That is, 
the recovery time begins again after any such pulse. This type of recorder is 
termed by Ruark and Brammer (1937) a type I recorder and that with a constant 
resolving time a type II] recorder. This isa convenient form of nomenclature. A 
common example of a type I recorder is the so-called Cenco counter. 

Consider now a simple counting system consisting of a counter and amplifier 
feeding directly a type I recorder of natural resolving time 7. ‘The resolving time 
of the counter and amplifier are assumed negligibly small. 

Consider a stream of pulses entering the recorder. Ifthe interval between two 
successive pulses is less than 7, the recorder is re-excited and the pulse lost. If 
the interval is greater than 7, the pulse is recorded whether the previous one has 
been recorded or not. Hence the probability of recording a pulse is equal to 
the probability, e-’’, that the interval between two successive pulses is greater 
than 7. This will be equal to the fraction of pulses recorded and consequently 
the efficiency is 

Bat, . 3 wel ee eee (20) 
This argument js due to Ruark and Brammer (1937) who attribute it to Volz (1935) 
and Schiff (1936). . 

The argument can be immediately extended to the general counting system 
in which the counter and amplifier have a constant resolving time o and feed a 
scaler of scale-factor n actuating a type I recorder. Here the recorder registers 
only groups of 7 pulses. A group of will be recorded only if the last pulse of the 
group occurs in a time greater than 7 measured from the last pulse of the previous 
group. ‘This is the case whether the previous group has been recorded or not. 
That is, 1f G*(n, 7, a) is the probability that the group of ” pulses after a given pulse 
occupies an interval greater than 7, the fraction of pulses recorded is G*(m,7, ¢). 


As before, every pulse is, on the average, associated with (1 + ) particles. 
Hence the efficiency is 


E=nG*(n,7,0)/n(1 + po) =G*(,7,6)/1 +o) nee (21) 
The function G* is given by 
G*(n, 7, ¢)= 1 —L*(n,:7, 0), 
where L*(n,—7, c) is the function considered above. That is 


G*(n, 7, 6) =e 4" 1 + le no)... [u(r —no)]"4/(n —1) 1}. 


Equation (21) applies equally well in the case where o is zero. 

The above formula for E was given by Ruark and Brammer (1937) fore =a 
and by Alaoglu and Smith (1938) in the general case. Ruark and Brammer 
explicitly derive their formula for a recorder of constant resolving time, that is, a 
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type Il recorder, for which it certainly does not apply. Alaoglu and Smith refer 
to a “recorder which will not respond unless the interval between the count and 
its predecessor is greater thant”. This statement applies equally well to type I 
and type IT recorders. Their calculated efficiency, however, is valid only for a 
type I recorder. 


§7 

An experimental verification of the expression given by Ruark and Brammer 

has been published by Lifschutz and Duffendack (1938). This expression is 

Bee no): 

What was actuallv done experimentally, however, was to verify the formula in a 
case for which it was not derived. In view of the confusion which seems to exist, 
it is as well to quote the remarks of Ruark and Brammer (1937, p. 324) : “If7, is 
greater than 7, and the recorder of a kind which cannot be re-excited whiie it is in 
action (type II), the efficiency of the apparatus is f’/f=(exp [ —f(r,.—7,)])/(1+f;)”. 
Here f’/f is our E, 7,=7, 7;=0 and f=p. 

In the experimental paper, the formula was interpreted as being applicable to 
recorders which are re-excited during recovery (type I). This is shown by the 
following quotation (p. 717) : “‘ In order to test the theory, the circuits used must 
satisfy the assumptions made in deriving the formulae. In referring to equation 
1 (a) it is clear that with r>o and o> 0 equation 1 (a) becomes n, = ny exp ( —M7). 
In other words the recorder circuit used to test equation 1 (a) must be a type I 
recorder’. Equation 1(a) referred to is Ruark and Brammer’s formula given 
above. A slightly different notation is used : m,/ny=E and ny=p. 

The experimental work was carried out employing a type I recorder (a Cenco 
counter) and very satisfactory agreement was obtained with the theoretical formula. 
As shown above, the formula in fact should apply only to a type I recorder. 

§8 

In practice, experimental conditions are usually adjusted so that the counting 
rate is sufficiently low for losses to be small. In these circumstances the efficien- 
cies of both type I and type II recorders are the same. In figure 2 the region of 
low counting losses is plotted in some detail. ‘The ordinates in this diagram are 
values of 100p*/(p* +7), that is, the percentage loss of pulses in the recorder. 
The abscissae are rates of entry of particles into the system. Logarithmic scales 
are used on both axes. Curves are drawn for various combinations of selected 
values of the scale-factor n and the effective resolving time (r —70). 

It can be seen that for a given value of », the curves for different values of 
(+ —no) are similar in shape, and can be obtained from one another by a simple 
displacement in abscissae. This is a consequence of the employment of 
logarithmic scales and the fact that the losses are functions of the product 
u(t —no). 

The losses due to the counter and amplifier can be obtained from the formula 
(18) for x=0. The total loss is the sum of those due to counter and recorder 
separately in the region of small losses. Inspection of figure 2 shows that the 
curves are all very nearly straight lines whose slopes increase with the scale- 
factor. The logarithmic plot emphasizes this effect rather strongly. ‘The steep- 
ness of the curves for moderately large scale-factors does, however, indicate that 
where there is uncertainty in the value of (t —”o), any attempt to make a correction 
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for appreciable counting losses ” 
would fail. Further, any fluctuation 
in the rate of entry of particles will 
tend to introduce large errors into 
any such correction. For large 
scale-factors therefore, it would 
appear to be desirable to arrange 
working conditions so as to make 
recorder losses negligible (e.g. less 

than 0-1%). Thatis, the maximum 4, 


8/a)| 10) 86) 32; 50, 


Average Loss (per cent) 


i00 i000 


0 
counting rate employed should be Particle Entry Rate (number per second) 
less than that corresponding to the Figure 2. Logarithmic plot of the percentage loss of 
point on figure 2 at which the pulses as a function of the average rate of 
relevant curve cuts the axis. hese entry of particles into the counting system. 


The curves are drawn for various scale factors 


PPA AMIS GUC. AY essential for indicated in the diagram, with a value of 0:1 


scale-factors of 1 or 2, and here it sec. for (r—n«), the effective récording time of 
may be possible to make small cor- the recorder. The two curves 8 (a), 8 (0), with 
rections if the parameters of the values of (r—no) equal to 0:1 and 0:05 sec. 


: respectively, and scale factor 8, show that a 
Sy Steps ate sufficiently accurately change in (t—7¢) results in a parallel displace- 


known. ment of the curve. 

In addition to fixing a safe 
counting rate, it is also necessary to fix a minimum total counting time. This 
must be sufficiently large to ensure both a small fluctuation in the total number 
of particles entering the system and a sufficient-number of counts. 


STALE PEN Dace 


The two tables given below summarize, in a form suitable for practical appli- 
cation, the data concerning the performance of various kinds of counting-devices. - | 
The authors are greatly indebted to Dr. C. E. Wynn-Williams for suggesting the 
form in which these tables should be compiled. 


Table 1 


The table gives critical mean rates of entry of particles per second for scalers of various 
scale factors associated with recorders of various resolving times. For entry rates less than 
these, the losses due to finite resolving time of the recorder are less than 0-194 and can 
usually be ignored in comparison with the losses due to the counter or other device at the 
input of the counting system and the statistical fluctuations. As the particle entry-rate is. 
increased above these values, the fraction lost increases so rapidly, particularly for high 
scale factors, that it is unwise to attempt to apply a loss-correction. The working rate 
should therefore be kept well below the critical rate. The critical entry-rate for other 
values of the scale factor can be determined sufficiently accurately by interpolation on a 


graph of the values given in the table. Note that the critical rates are inversely proportional 
to the resolving times. 


Recorder Scale factor 
resolving time 

(seconds) 2 4 8 10 16 B2 64 100 
1/5 0-25 2 10 15 32 86 210 360 
1/10 0:50 4 20 30 64 172 420 710 
1/25 1-2 10 50 75 160 430 1000 1800 
1/50 2:3 21 100 150 320 860 2100 3600 
1/75 3:5 32 150 225 480 1300 3200 5300 
1/100 4:5 42 200 300 640 1700 4200 7100 
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Table 2 


The table gives maximum permissible mean rates of entry of particles per second into: 
a scale-of-one device for different values of the mean fraction of particles lost and for various. 
resolving times of the device. In brackets are also given the minimum time of counting 
(in seconds) at these maximum rates for the statistical fluctuation in the count not to exceed 
the instrumental losses. The table shows that it is difficult to work to an accuracy much 
better than about 0-5°% unless long counting times are employed. 


Type of Resolving Percentage of count lost by scale-of-one device 


device time (sec.) 01 0-2 05 1 2 | 5 | 10 
Input stage iQ) 1000 2000 5000 | 10,000 At higher entry-rates an 
of hard (1000) (125) (8) (1) adequate count is se- 
valve scalers 7 a a 1000 2500 5000 cured in less than 1 sec. 
limit about (250) (16) (2) but special timing de- 
O-? sea. SOR? 400 1000 2000 vices are required. 
| (625) (40) (5) 
Input stage Ome 500 1000 2000 
of thyratron (80) (10) (2) 
scalers BENO== 250 500 1000 
limit about (160) (20) (3) 
10-5 sec. 
' Geiger- Sal 100 200 400 
Miller (400) (50) (7) 
counters IO 50 100. 200 
limit about (800) (100) (13) 
DGlOmesec: 2am At lower entry-rates count- 50 100 260 
ing time must exceed (200) (26) (2 
S10 1000 sec. (about 16 min.) 20 40 100 
otherwise statistical fluc- (500) (64) (4) 
tuations may exceed 
instrumental losses. 
Ome 10 20 50 
(1000) (130) (8) 
ENO 10 25 DD 
(260) (17) (2) 
Modified Be alms a 10 22 
ratchet (640) (40) = (3) 
uniselectors IO) 5 11 
limit about (80) (10) 
Seal Omerseea ia <ul Om? Mek 5 
(160) (20) 
Modified 5x MOS 1 2 
telephone call (400) (50) 
counters 104 0-5 1 
limit about (800) | (100) 
A Os2'sec. 
108 Dre 102 4210" 104 2-596 108 4107 | 02 


Minimum desirable count. 
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DISCUSS 1 ON 


Dr. Denis Taytor. The computation of counting losses is important both from the 
point of view of arranging the conditions of the experiment so that the losses may be 
reasonably small and also so that the true counting rate may be calculated from the observed 
counting rate. The analysis given by the authors allows this computation to be made for 
most of the practical cases which arise in nuclear physics. One case which is of some 
importance, and which seems to have escaped mention, is the case when the resolving time 
of the scaler itself must be allowed for. In the simple case when only the resolving time of 
the input stage of the scaler is important, it is presumed that the analysis given by the 
authors applies, but that o, the resolving time of the detector (ion-chamber or Geiger- 
Miller counter), must now include the resolving time of the input stage of the scaler. Other 
cases can arise which are more complicated. A typical example is a pulse analyser which has 
been designed recently at T.R.E. This consists of a pre-sorting stage, a “ scale of two ae 
a ‘ring of five’? and an electro-mechanical register. The resolving time of the input 
circuit was 100 psec., that of the second stage (i.e. after the “‘ scale of two’’) was 2 msec., 
and that of the electro-mechanical register was 0-1 sec. The second stage had a high 
resolving time because cold cathode valves were used for the “‘ ring of five” to economize 
in power supplies (30 such chains are used in the complete equipment). In a case like this, 
the finite resolving time of the second stage of the scaler must be allowed for. Have the 
authors considered extending their analysis to cover this case ? We have made some vain 
-efforts to do this at T.R.E., but eventually obtained the information we required by connect- 
ing a Geiger-Miiller counter to two counting chains, one of the type described above and a 
second one in which the resolving time of the second stage of the scaler was sufficiently small 
to have a negligible effect on the counting rate. The counting losses for the second chain 
could be computed and, therefore, a comparison of the observed counting rates for the two 
chains allowed the counting losses for the former to be determined. 

Another case of considerable importance is the computation of counting losses in the 
counting of particles from pulsed sources. The use of particle accelerators which give short 
pulses of particles, usually of a few seconds’ duration, repeated from a few to a few hundred 
times a second, as, for example, with a synchrotron, frequency modulated cyclotron, etc., is 
likely to increase very much during the next few years, and the importance of being able to 
compute the counting losses in such cases cannot be overestimated. Have the authors given 
any attention to this problem ? 

Finally, I should like to call the authors’ attention to a further paper on the computation 
of counting losses which may have escaped their attention. This is Hole, Arkiv Mat. 
Astr. Fys., 1946, 33 A, No. 11. 


Dr. G. BuRNISTON Brown. I should like to ask whether any experimental confirmation 
of the formula arrived at by the authors has been attempted? Physicists must not be 
prevented, by the fact that amplifiers etc. were mentioned, from noticing that this is a 
mathematical paper whose conclusion is nothing more than the consequence of assumptions 
made at the start and during its progress. 'The result would have been just as correct for 
elephants as for amplifiers, provided we assume that if m stimuli are given to an elephant, a 
fraction p of them produce a “ kick ”’, etc. Anyone making use of the theory of probability 
should bear in mind that it is the most obscure of all theories, and that Leibniz, the late 
Lord Keynes, and Dr. Harold Jeffreys all agree in admitting that they cannot define the 
concept of probability. 

Why should modern physicists be excused from the venerable custom of calibrating 


their instruments ? Arguments from the theory of probability can never take the place of 
experiment. 


AvuTHors’ reply. In reply to Dr. Taylor, we have not so far considered the extension of 
our calculations to the case where the second (and further) stages of the scaler have appreci- 
able resolving times, but we think it probable that the methods used are capable of this 
extension. The problem of counting losses from pulsed sources is of a different type. It is, 
for instance, obvious that when the resolving time of the counting system is longer than the 
duration of a pulse only one particle can be counted per pulse. It is therefore necessary to 
modify even the most elementary formulae for counting losses. We are considering this 
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calculation in view of its general interest and are indebted to Dr. Taylor for calling our 
attention to this case. 

With regard to the point raised by Dr. Brown, we have drawn attention in our paper to 
an experimental verification of one of the formulae discussed, namely the work of Lifschutz 
and Duffendack. It would be possible in principle to calibrate a counting system empiri- 
cally if standard sources of suitable types were available in every laboratory. In some cases, 
however, this calibration would involve considerable difficulties, e.g. the efficiency of a 
Geiger counter for y rays varies with the energy of the y ray, and a whole range of standard 
y-ray sources would be necessary. 


Certain Aspects of the Mechanism of Spark Discharge 


(BYol a Bas IORB 
University of California, Berkeley, California, U.S.A. 


MS. received 12 November 1947 ; read 23 April 1947 | 


PARKS area class of transient occurrence in which a given existing conduction 
current in a gas suddenly and irreversibly changes to a current of higher 
magnitude, operating more efficiently by different mechanisms under 

conditions which rendered the previous lower current unstable (Loeb 1939, p. 408). 

Certain general features of the threshold for spark breakdown, common to all 
mechanisms but differing in detail for the different processes, will first be presented. 
The more complete equation for sparks in one regime will then be used to describe 
the general properties which may be expected. 

Spark breakdown in general requires fields of sufficient intensity to ensure 
that two processes, a primary and a secondary, are activated. ‘The primary 
process is one which causes electron multiplication in the field by electron impact 
to reach such values that after traversing the gap of length 6 the exp(«), or 


0) 
exp (J ads), electrons created by one initiating electron are adequate. ‘This 
/ 


is the primary, and predominating, process, with «, Townsend’s first coefficient, 
representing the number of new electrons created per centimetre of path in the 
field direction by an electron (Loeb 1939, Chapter VIII and p. 369). Alone this 
process leads to field intensified currents which are proportional to the initial 
ionization present. ‘The progeny of exp (ad) electrons produced by one electron 
in crossing a gap is called an electron avalanche. ‘The average distance traversed 
in the field direction to produce an ionizing act is then 1/«. ‘The function of the 
secondary process is to furnish secondary electrons to replace the initiating 
electrons and to maintain or build upa discharge. It can take one of two possible 
forms (Loeb 1939, Chapter IX, pp. 377 and 403): 
(a) A secondary process at the cathode leading to Townsend’s type of break- 
down, (i) by impact on the cathode of positive ions that liberate electrons ; 
(ii) by photoeleciric effect on the cathode of photons excited in avalanche 
formation; (iii) by action of metastable atoms on the cathode (limited to 
certain gases). 
(b) A secondary process im the gas that leads to sparks by the streamer 
mechanism, i.e. photo-ionization. 
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To explain the behaviour observed, the liberation of electrons at the cathode 
by positive ion impact on the cathode is indicated as the most likely process. 

As Townsend has shown, this leads to a current 

i=ile@/1 —ye)|,— = 3 eee (1) 
when expressed in a simplified form (Thomson 1911, Townsend 1915). 

When ye=1 this expression becomes indeterminate mathematically. 
However, as Holst and Oosterhuis (1923) long ago indicated, ye” =1 has another 
interpretation. ‘This becomes clear if one reasons as follows :— 

(a) For ye <1 the discharge follows equation (1) and the current is not self- 
sustaining; i.e. it depends on 7 and ceases with z. It yields merely a 
field intensification of Zp. 

(b) For ye* =1, on the average, each avalanche of e”’ electrons by the primary 
process is multiplied by a value of y large enough,to give one new secondary 
electron when the e” positive ions return to the cathode. ‘Thus one 
initiating electron at the cathode is able to maintain its succession inde- 
finitely. This is then the threshold for a self-sustaining discharge inde- 
pendent of 7. It marks the sparking threshold. 

(c) For ye*®>1 the ionization of successive avalanches is cumulative and more 
ions are created than start. As electrons have a mobility 100 times that of 
positive ions, positive space charges accumulate in the gap near the 
cathode. These produce a very efficient secondary mechanism that 
allows of the discharge continuing as glow or arc at a higher current. 
The space charge will grow and the spark will materialize the faster the 
greater the amount by which ye*? exceeds unity. 

The mechanism of such a.spark then proceeds as follows. An avalanche 
crosses a gap of length 6 at a speed of the order of 10’ cm/sec. in sparking fields at 
atmospheric pressure in a time 7. Arrived at the anode the e electrons are 
absorbed and the positive ions, of which 50% are created within the last 1/xcm., 
start back to the cathode. Arrived there in a time T.,, roughly 100 times as great 
as T., these ions yield ye electrons which start a new chain of avalanches, again 
taking a time T, and T,. At this time ye?” electrons begin a third sequence. 
When, as a result of 7 trips, the positive ion space-charge of y"1e” positive ions 
produces a self-sustaining field for a glow, the spark is complete. Its time of 
formation, or formative time lag is 

P= (Lit Te) eee (2) 

If the secondary mechanism had involved photoelectric action at the cathode, 
the time 7; would not have involved T,, for the photons reach the cathode without 
delay. ‘Thus 7; would be merely 77,. What evidence there is, however, 
indicates that photoelectric emission at the cathode is largely ineffective because of 
the wide geometrical dispersion of photons, for it will be rare indeed that a photon 
is liberated near the origin of the initial avalanche. Thus, in general, the y 
mechanism is responsible for sparks and the value of 7; is given by equation (2). 

The question has been studied by Schade (1937) using equation (2) with 
‘Townsend’s mechanism active at lower pressures and has been verified as being 
in accord with observation. For gaps in Ne and H,, 7; has a value of 10-1 to 
5 x 10~°, and has been shown to depend on 7, and 7. 

Near sparking in air at atmospheric pressure the quantity «5 has the relatively 
small value of about 17. Furthermore, ionization is a chance phenomenon and 


Certain aspects of mechanism of spark discharge 563 


it is not certain that every electron will ionize each time that it has traversed 1/« cm. 
There is also nothing to prevent electrons by accident from ionizing twice in a 
distance 1/«, even though on the average it does so only once. Thus the individual 
avalanches will yield an electron multiplication which fluctuates above and below 
e* in a purely random fashion, i.e. if «5 varies by +1 the avalanche may vary bya 
factor of e or 1/e. Again electron liberation by positive ion impact is purely a 
chance phenomenon with an average value y under any set of conditions. 
Thus the quantity (ye*’),, as yielded by individual avalanches will fluctuate about 
the average value ye” in considerable measure. 

Accordingly even though the potentials are chosen so that the average ye*® = 1, 
individual avalanches occurring will yield other values of the product. Thus if 
V, is the potential corresponding to ye” = 1, an occasional fortunate avalanche will 
build up a spark at a potential V less than V,, and at a potential V greater than V, 
- an occasional unlucky avalanche will not give the necessary perpetuating electrons 
and there will be no spark. It is clear then that in individual avalanches sparks 
occur at values of V below and fail at values above the value of V fixed by ye” = 1, 
through the relations of «/p to X/p. 

If now there be applied to a gap a potential V with initial value less than V, 
and the fraction of P, initiating electrons which give avalanches resulting in the 
formation of a spark be plotted against V, the curve A of figure 1 would be observed. 
This curve is the integral 
from V=0 to V=V of the 
fluctuations of the value of 
(ye), for individual elec- 
trons about the average 
value for which ye” =1 and 
V=V,. Thus differentia- 
tion of.the eurve A_of 
figure 1 gives the chance of 
a variation of (ye); in a 
sequence of avalanches 
about the mean value as 
expressed in terms of V. 
The differential curve is 
shown in figure 1, curve B. 


Fs 


: Figure 1. 
(A) Probability of sparking Ps plotted against applied potential 


If the derivative gives a cer Ons, 

symmetrical distribution, (B) Derivative of sparking probability against voltage indicating 

V. will be the maximum the distribution of sparking events about the threshold 
8 


value, and in the integrated DUN Seta SOROS 


curve it will represent the point of inflection at P,=0-5. With values of «5~20 
this is nearly true. 

Curves of the type of curve A in figure 1 will be obtained if the fraction of 
individual avalanches leading to a spark are observed and plotted against V. 
By measuring time lags using varying but known ultraviolet illumination of the 
cathode and different potentials V, Wilson (1936) obtained curves analogous to 
figure 1 curve A and used these to evaluate V,. ‘The width of curve B of figure 1 
depends on the magnitude of the fluctuations in y and e” and on the number of 
sequences 7 to build up the space charge. While any one of the steps in the 
sequence of 7 successive events could interrupt the spark, it is probably only the 
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early ones where not many positive ions are present that are most critical. ‘Thus. 
at V, where 7 is large the curve may be somewhat broader than at higher values of V. | 
The width of curves calculated for Schade’s data by Hertz (1937) using fluctuations 
of y only lie around 1%. Actual widths may be much larger as fluctuations. 
in are larger than in y. No experiments using individual electrons have been 
made. 

Whether the distribution be wide or narrow, it is important to know that this 
distribution is always present and that sparks will sometimes occur below a 
potential V, based on average values and will sometimes fail above. ‘The value of 
V., unless the curves are very distorted, will be that at which P,=0:5. Thus, 
to be certain of the sparking threshold experimentally, a (P,—V) curve must be 
obtained from individual avalanches and V, chosen as above. It should further be 
noted that 7; varies with V about V, and that properly to interpret formative time 
lag studies these lags must be established relative to V,. 

When a potential V is applied, as in sparking potential measurements, it is 
observed that even at V,, or above, the spark does not occur at once. ‘The time 
between applications of V and the appearance of the spark is called the t7me lag of 
sparking. We shall call it the observational time lag, Ty. In 1925 Laue and 
Zuber showed it to be composed of a statistical time element T,, depending on the 
chance of arrival of an initiating electron and the chance that this electron causeda 
spark, and of a formative time element, 7; which is the time taken for the spark to | 
grow, i.e. 7y>=7,+T7;. Inturn, T, depends on T, = 1/n, the average time between 
liberation of initiating electrons, which is the reciprocal of the number of electrons. | 
per second, and the chance P, that these events will give a spark. That is, 
T,=T,/P;=1]/moP,. Thus Ty>=1/n)P,+T:. But by equation (2) T; on Town= 
send’s theory is 7;=(T,+T7,)y. With gap length 6, field strength X, and ion |} 
mobility K,, 7, is given by T, =8/K _X, and T, is given by T,=5/K_.X, where K, | | 
is the electron mobility. Thus , 

T= (8/ KX +07 KX jp (3) 

It must be also recognized that 7 decreases rapidly with increase of V above V, as 
a/p changes rapidly with X/p. ‘Thus for overvoltages (_V —V,)/V, which are not | 
very large 7 approaches unity and 7; approaches 7;=8/K,X+8/K,X. Hence | 
for Townsend sparks we can write | 

Dyo=1/mpPs OK 2X +0/K Xa ee (4) 
If now T,, = 1/n) be made small by large m), and if P, and 7 be made unity by using | 
a sufficiently large overvoltage (V -V,)/V,, then TJ, becomes sensibly 


Ty=3{KyX, ..- eee (5) 


since K,>K,. For rare sparks below V, observed with large no, the value of T; | 


will be slightly larger than that of equation (3) due to the decrease in XK,. The |f} 


conditions given by equation (5) are those under which most experimental studies | 
of both V, and 7, have been made—except those whose purpose it was to study T,. 
This means that true values of V, have rarely been observed, but with large ‘i | 
and the narrow statistical range of V about V, the values are not seriously wrote | 
It also means that by using large values of (V —V,)/V, the quantity T, observed 

has been the characteristic value of T>=8/K.,X for the Townsend mechanism. | 
Observations of 7) near V, would give exceptionally long and non-informative | 
values of T;. It may ke noted here, however, that values of 7) must not be too high 
since various workers (see Loeb 1939, pp. 386, 448-9 Loeb and Meek 1941, 
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pp. 32, 172) have shown that a heavy ) with much field intensification can 
cause a space charge distortion that lowers the sparking potential or of itself leads 
to a spark. 

It may also be noted that mp) can be determined, even at breakdown fields, by 
measuring the value of the photoelectric current at two potentials below cumu- 
lative ionization and evaluating the constants in J. J. Thomson’s equation for 
photo-ionization currents in gases (Loeb 1939, Chapter VII, Sect. 3, p. 310 and 
equation, p. 313; Johnson 1948). This enables the value of the current from 
the cathode, 1, to be calculated for purposes of study. 

It must finally be concluded that to evaluate V, one must observe 7, under 
the following conditions : 7; must be kept much less than 7, and 7’, is increased by 
a controlled my which must be sufficiently small. Then sensibly 7) =1/n)P, and 
by observing 7) for different values of V about V,, P, can be evaluated as a function 
of V and by differentiation the peak value will locate V,. 

Practically all studies except Schade’s have so far been made on sparks passing 
by the streamer mechanism and not by the Townsend mechanism. It is thus 
essential to see how the considerations above apply to the streamer. 

The streamer starts with an avalanche of e** electrons where x is nearly equal to 
the gap length 6 at V,._ If e*” is great enough, electrons of the avalanche passing 
to the anode leave behind an adequate positive space charge to initiate a streamer. 
Meek (1940) and Raether (1941) independently assessed this as occurring when 
the space charge tip field X’ equalled the imposed field XY. ‘The density is deter- 
mined by e*” and the radius of the avalanche head p produced by diffusion of 
charges during advance. ‘The diffusion is restricted in scope so that e*” and the 
pressure primarily determine the density of charge. ‘To enable a streamer to 
advance, the space charge tip field plus the imposed sparking field must be great 
enough over an adequate small volume element Ax’ deep and solid angular 
aperture 77/2, so that a photoelectron produced in this volume in advancing to the 


0 ae 
tip can extend the space charge by making exp (| a’ a) =exp (| o dv) s abiere 
Av’ 0 


w is the value in the undistorted field X where the electron avalanche travels a 
distance x to start a streamer and «’ is the much higher value in the vector field 
X +X’ between Ax’ and p._ If several photoelectrons can be produced in Ax’ the 
volume can be smaller or X’ less. ‘The production of photoelectrons in a distance 
Ax’ depends on the density of photon production in the avalanche tip and on the 
absorption coefficient of the gas for wavelengths which produce photo-ionization, 
as well as on Ax’. Photon production varies directly with the number of electrons 
e*” and can be set as fe%”, where f is a nearly constant multiplier which may be 
greater or less than unity, but probably less. It should be noted in passing that 
the density of the molecules in these mechanisms will play a very critical role since 
p and pare pressure dependent, so that pressure will be an important parameter, 
apart from the role of X/p. fed 

It is now clear that while the avalanche formation 1s common to both Townsend 
and streamer mechanisms so that only one initiating electron is needed in both, the 
probability y of the Townsend mechanism is replaced by the chance of photo- 
ionization in restricted volumes Ax’ in the streamer mechanism which might be 
lumped in a factor e. Thus streamer progress likewise depends on fe*” for photon 
production, multiplied by a probability factor ¢ akin to y, which is now complex and 
involves yz, Ax’ and the criterion of Meek and Raether. ‘Thus one may set the 
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sparking threshold by streamer mechanism as <fe*” = 1 for streamer propagation. 
. It is at once clear that with streamers the statistical time lag T, will, by its depend- | 
ence upon m, and P,, play the same role as with Townsend’s mechanism except 
that now P, will depend on other factors lumped in « such as p, Ax’, and pressure. 
T, will show differences with streamer formation. As, however, the condition 
cfe*” =1 is radically different from ye”? =1, the analogy is largely formal. 

It was the early discovery, 1923-27, that the value of 7; for higher pressure 
sparks is much less than is possible on Townsend’s picture (Loeb 1939, pp. 32-33, 
Loeb and Meek 1941, p. 449) that, together with the help of corona studies — 
(Loeb 1939, pp. 433, 440-450) and cloud track pictures (Raether 1935, Flegler | 
and Raether, 1936, 1938, 1940), led to the streamer theory. In this case, the time 
scales indicated that 7; depended primarily on the time of electron crossing, 7... 
For ence the electron travels a distance x, to near the anode, the streamer forms and 
progresses back to the cathode with a velocity several times that of the electron. 
This high speed of the streamer advance results from two conditions. ‘The 
enhanced tip fields of positive and negative space charges in the streamer increase 
electron velocities, and photo-ionization in advance of both space charges leads to 
an increased velocity, because photons move with the velocity of light. These 
increases are of the order of the ratio of the square roots of the enhanced tip fields X’ 
to the applied field X, from 4/2 to 1/10 according to Raether’s computations _ 
(Raether 1940), and of the order of the ratio (Ax’ +p)/Ax’, where p is the radius 
of the avalanche tip and Aw’ the width of the sensitive part. These factors can 
increase the velocity of the streamer by a possible factor of 2 to 10 fold above that 
of the avalanche advance. ‘They depend critically on pressure and on applied 
field strength. Accordingly for streamer advance T;=x/XK,+.x/bXK,, where 
is the multiplication factor y/(X’/X) . (Ax’ + p)/Ax’. The final stage of spark 
breakdown after the streamer has crossed the gap is very short, for it is determined 
by the velocity of re-ionization of an ionized streamer channel by the steep potential 
wave proceeding up the streamer channel from the instant of junction of streamer 
head and cathode. ‘The velocities, as shown by Snoddy, Dietrich and Beams 
(1937), by Allibone and Meek (1938a,b), and by Schonland et al. (1934, 1935, jf 


1937) are of the order of 10*°cm/sec. The time lag of the return stroke is 7). 
‘Thus for a streamer, 


Tp=x/XK,+x/bXK,+ T,. 


Since b>1 and 7, is very small, one has sensibly T;=x/XK,. Here again with | 
increased potential, i.e. overvoltages, 6 is greatly increased. Furthermore, with | 
overvoltages streamers form in mid-gap at a distance x from the cathode much less | 
than 6 as « increases rapidly with V. Such streamers have been observed with |] 
Kerr cell shutters by Dunnington (1931) and White (1934) in overcharged gaps. | 
Hence, 7 will decrease with (_V —V,)/V, very rapidly and continuously in the | 
proportion «/6. 
Thus with the streamer mechanism statistical and formative time lags appear. | 
There is, however, one difference apart from the fact that the value of T; for | 
streamers is much less than that for Townsend’s discharge (8/X,K,<8/X,K,). ] 
This is that 7; has a fairly sharply defined upper limit, 8/XK,+5/bXK,+ T,,.forg| 
streamers, and decreases rapidly with V above V,; in a Townsend discharge on |} 
the other hand 7; is indefinitely great at V, and, for slight increases in V above Vs 
decreases rapidly towards a minimum value 5/XK _,, the characteristic value. 
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In a further development of the streamer theory of spark discharge, one of the 
critical data yet to be determined is the condition for the change over from a 
Townsend to a streamer mechanism (Loeb and Meek 1941, pp. 49, 71, 76). From 
a knowledge of this much information on streamer propagation could be derived; 
it is hoped that some of this lacking information may be provided by the researches 
now being carried out by Meek and Craggs. It is indispensable for sparking 
theory, however, that the transition be observed as a function of pressure and 
gap length. ‘The only way in which this can be done today is by means of a 
careful study of time lag, thus establishing the value of P, and the position of V, 
from a knowledge of Ty and mp. Then by increasing ny so that 1/7)P,< 7; one can 
observe Ty when Ty = T;, and so determine 7; as a function of (_V —V,)/V,. Then 
from the value of 7; observed in relation to 6/XK, and 8/XK, the operative 
mechanism is clear. ‘his must be done with 8 constant for different pressures p 
and with p constant with different gap lengths 56. 

That such transitions will not occur at a single value of the product pd for a given 
gas is indicated by the fact that unlike the Townsend mechanism, neither the 
streamer mechanism nor a pure space charge breakdown follow Paschen’s law 
(Loeb 1939, p. 410, Loeb and Meek 1941, pp. 8, 46, 72, 81, Varney et al. 1935) 
which states that V,=/(pd). Physically interpreted, the failure of Paschen’s law 
and the principle of similitude with streamers follows from the fact that a Paschen’s 
law breakdown depends on the number of ions created in the gap, while the streamer 
mechanism depends on zon concentrations. ‘The failure of Paschen’s law has not 
been observed at ordinary pressures, as the present day accuracy of measurements 
of V, with the peculiar mathematical character of the sparking equations, either 
‘Townsend or streamer, are not adequate. ‘Thus characteristic data in the tran- 
sition region are urgently called for. It is essential, however, at this point to study. 
the character of the departure from Paschen’s law in streamer formation as shown 
by Meek’s theory alone in order to proceed further.* | Meek’s equation may be 
written as 


a) Bi Gln ein) ain Maa (7) 


Vs 
(2) vom) -00 
where V,=X,5 for a uniform gap. Here it is seen that V, is a function of pd as 
Paschen’s law predicts up to the last term Ind. ‘The Meek equation again 
takes the form, 


ad+Ina/p=14-46+InX/p+tind/p — ...... (8) 


with ; = (5) = (5) =i (55) cht: ae ry (9) 


The form of the function f(X/p) for lower values of X/p in the streamer region is 


either 
a/p=Aexp(BX/p), or a/p=A((X/p)+B)?, ...... (10) 


the latter holding for the higher range of X/p. What is done in solution of Meek’s 
equation for V, by trial and error is as follows : For a given p and 8 a value of X'/p 
is chosen and the corresponding value of «/p obtained from the curve for 
a/p=f(X/p). The values are then inserted into equation (8) and the relative 

* It must be noted, however, that apart from the concentration effects in Meek’s equation the 
value of p is also independently pressure dependent so that the pressure dependence may be much 
greater than calculated above. 
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magnitudes of the two sides of the equation observed. A value of X/p is next 
chosen to make the other side of the equation larger. Thus by successive approxi- 
mations, or by the use of Simpson’s rule, a value of X/p can be obtained to satisfy 
the equation. This value at X/p is X,/p and yields V, from X,=V,. Itis at once 
seen from equation (8) that if one calculates V, for 5 constant and p varying or for 
p constant and 8 varying, different values of V, will be obtained for the same: 
product p3. This is, of course, because if 6 is increased and p decreased, In (3/p) 
increases, while for the same pd if p is increased and 6 decreased, In (6/p) decreases.. 
The changes however, are small because of the logarithmic function. ‘Therefore, 
for any one value of pé there will be two values of X,/p and hence of V, depending 
on whether the value pé is reached by decreasing 6 and increasing p or vice versa. 
The difference between the two values will be the greater the greater the range of 
variation in 5 and p, but it will be relatively small in its effect on V,, for V, depends. 
on « (equations (9) and (10)) which varies rapidly with X,, while compensating 
-changes in 5/p in equation (8) would have to be very great since they alter the’ 
left-hand side by only In8/p. Experimentally, therefore, it would be only for 
larger ranges of variation that deviations from Paschen’s law would be detected 
experimentally. 
Failure to recognize this fact has in the past led to misunderstanding. When 
Meek first studied his equation he put it to test by calculating V, as a function of 
po, using the value of 6 at 760 mm. Hg initially and decreasing p. In consequence, 
he obtained a V, curve in which In (6/p) increased as p decreased. This required a 
larger value on the right-hand side of the equation with the result that he plotted 
a (V,, pd) curve which was high. He noted a marked departure of the (V,, pd), 
curve from the experimentally observed values of Whitehead, at p~200 and 
assumed that this was the point at which Townsend’s mechanism changed to. 
streamer mechanism (Meek 1940, Loeb and Meek 1940, 1941, p. 49). Roughly 
Dunnington’s spark photographs with Kerr cell shutter also indicated a similar 
value as the point of appear- 
ance of mid-gap streamers 
in an overcharged gap 
(Dunnington 1931). Actu- 
ally the appearance of mid- 
gap streamers depends not 
on the product pd alone but 
also on p and/or 6. Later 
L. H. Fisher (1946) com- 
puted V, as a functien of 
po from Meek’s equation 
using constant p and vary- m—eaeee 
ing 6. He thus found a 
curve which started at the Gos Gon 
same 6 as Meek’s curve for eet ted aa 


atmospheric pressure Figure 2. Ppehae the sparking potential plotted against 
ie 760 butwhich beeen fa ee nes: ‘Townsend 's equation and Meek’s. 
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data taken from Schumann’s book (1923). These data are, however, variable 
in their reliability, and it is necessary to take into account the method by which 
they were obtained. Presumably, they were obtained with decreasing 8 to conform 
to Fisher’s calculation. Thus Fisher noted no departure from observed curves, 
even at values where streamers must have given way to Townsend’s mechanism. 
This precipitated a rather complicated situation whose interpretation may be 
clarified by reference to the curve shown in exaggerated form in figure 2. The 
chosen point of departure for calculation is pS). __D, is the point at which the 
Townsend mechanism and Paschen’s law take control for p decreasing. Meek’s 
theory no longer holds. D, is the point at which Townsend’s discharge and 
Paschen’s law take control for 6 decreasing. Meek’s law no longer holds. The 
dotted curve is the Townsend curve with a y following Paschen’s law and repre- 
sents the observed curve up to streamer formation at the two different transition 
points. 

The minima (p35), for Meek’s two curves fall at the same value of pd in both 
cases as they reflect merely the point where «/p=/(X/p) passes through its point 
of inflection. The value of (pd),, will be lower for ‘Townsend’s theory as the 
quantity y is increasing markedly with X/p in this region of values of p6 and will 
certainly be more effective than the photo-ionization in the gas. 

With this understanding of what may be expected, one can say that the differ- 
ences 1n the calculations of Meek and Fisher indicate the character of the departure 
from Paschen’s law. None of the data which they use are of any value in this 
region of po to indicate this and thus accurately to confirm theory. » 

Fortunately, recent results of Skilling and Brenner (1941), and particularly 
of Howell (1939) and of Trump, Stafford and Cloud (1941) with a large Van de 
Graaff generator have unequivocally shown the failure of Paschen’s law. Mr. C. 
G. Miller has kindly recalculated and plotted some of their data to show this. 
Curves of V, plotted against p for constant pd should be horizontal if Paschen’s law 
holds. The curves of Trump, Stafford and Cloud (figure 3) show that with 
increasing pd and in- 
creasing V, there is a 
marked and increasing 
slope to the curves. 
This linear departure 800 
is difficult to corre- y 
late with the Meek w 606 
equations (7) and (8) KY 
above. If, however, 
we calculate 6 from 
the data of these 
curves and plot V, 200 
against 6 for constant 
ps, we obtain the 
curves shown in figure (o) 30 60 90 120 150 180 210 
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with po constant and Figure 3. Curves of Trump, Stafford and Cloud showing sparking 
: : Sl ee potential in ky. plotted against pressure in pounds per square 
MCLE ABLE z ee inch for constant values of p6. If Paschen’s law held strictly, 
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(7) indicates.* Here, then, we have our first indication of the deviation from 
Paschen’s law predicted by the streamer theory, and it is observed in regions 
where p and 6 are varied over large ranges of values. 

It should be added that in Howell’s data with very high p and high V,, above 
400 ky. from the gassy surface used, the appearance of field emission of electrons 
from the cathode could lead to more efficient breakdown mechanisms and lowered 
V.. All Howell’s data at above p5=35 have such fields and since these are too- 
large they distorted his gap and possibly lowered V,. The character of the curves 
reported and the 
consistency of the 
data with those of 
Trump (below field 
emission) indicate, 
however, no de- 
parture ascribable wy 
to field emission in VY 600 
this region. ‘This 
is doubtless due to 400 
Howell's experi- 
mental care and 
choice of results. 
The findings thus 
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tionto the streamer Figure 4. Curves of Trump, Stafford and Cloud as calculated by 
theory and the use C. G, Miller with sparking potentials in kilovolts against gap length 


: in inches for various value of pd for comparison with Meek’ - 
of Meek’s equation. ‘ cote Ss oe 
; tion. Meek’s equation will give analogous curves, but the variation 
There is, how- will be distinctly less than observed as indicated by later considera- 
ever, a definite tion of Meek’s sparking condition. 


limitation to be ex- 

pected in the application of Meek’s equation beyond that occurring at low pd and p 
where ‘Townsend’s mechanism is more efficient. ‘This applies particularly to large 
values of 5, where p remains relatively low, i.e. to long sparks at atmospheric 
pressure. ‘The limitation will also appear when somewhat lower pressures and. 
longer gaps are considered, i.e. for V, curves where the initial value of 6 is large and 
pis decreased. Meek’s present equation takes no account of the effect of lowering 
the density of photo-ionization through either the reduction of « at low X/p or in- 
creased diffusion (Loeb and Meek 1941, pp. 49,71, 76). Further, both density of 
ionization and absorption in the sensitive zone Ax’ rapidly decrease as p decreases. 
‘Thus before Townsend’s mechanism at low p3 would normally be more effective 
than streamers according to Meek’s theory, streamers could cease to form because 
of factors not included in the theory and sparks will occur on the Townsend mechan- 
ism. Should this factor cause transition as p is decreased or at low p as 8 is 
increased then the (V’,, 5) curve would show a slight increase in V, with scattering 
of data and lowering of pd as indicated by A in figure 2. Experimental data are 
as yet too crude to show this. 


* Recent calculation (Miller and Loeb 1947) indicates that in fact Meek’s theory predicts the 
sort of departure observed. It is, however, smaller than the observed departure (1% as against 


14%). This is because the pressure sensitive effect present in ¢ is not included in Meek’s equation, 
to its full extent. 
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It is much more important, however, to consider what happens from this 
circumstance in very long gaps and thus also in lightning discharge (Loeb and Meek 
1941, p. 79). If either Townsend’s or Meek’s equation be studied it is clear that 
except for slowly varying logarithmic terms, the chief variation is given by 
f(V,/p8)p5 which is nearly constant. It yields V,/p5=f-1(1/p3). This indicates 
a hyperbolic type of decrease of V,/p3 with increase in pS. Hence, since V’,/8 = X,, 
the sparking field strength, one obtains X,/p=f'(1/p8). Thus, the important 
sparking parameter X,/p should decrease approximately hyperbolically with 
increasing pd. In fact, from what data we have (taken mostly at atmospheric 
pressure with increasing 5), up to p5=104mm. x cm., this relation is observed as 
seen in figure 5. Carrying the data up to 3 x 104mm.xcm. from the data of 
Howell and of Trump, Stafford and Cloud, the value X,/p decreases nearly 
linearly, within observational error, as p increases to 40 atmospheres. Corre- 
lation of the data with pd, as required by the argument above, is not very good but 
there is a similar trend with large scatter. In this region, however, Paschen’s law 
deviations are extreme and the scatter with pd is caused by the failure of this law 
in keeping with Meek’s equation. For lower pressure regions where Paschen’s 
law deviations are not serious, the value of X,/p is 33 with X,=25000v/cm. 
At 40 atmospheres, the value X,/p is 21 with X,=650000 v/cm. 

The cause for this behaviour is seen if one writes the Meek positive space charge 
field equation in the exponential form, 

X _ A(a/p) exp ((«/p)ps} ( 
? (ip) Soe feed weet we etcwse: 
This indicates that, neglecting the non-exponential factors, the equality, and thus 
the condition for sparking, is maintained if the value of «/p decreases in proportion 
to increase of pd. ‘Thus, by increasing 6 at p constant as in the conventional 
studies, «/p can decrease, and hence the 
value of X,/p needed can decrease and 
yet give a spark by Meek’s equation. i 
This, interpreted physically, means 
that with increasing gap length owing 
to cumulative ionization the number of 
ions formed per centimetre of length 
to give the space charge field can de- 
crease and yet yield a spark in a gap of 
length 5. It is true that «/p as a 


iF = > |x 
a 
° 


30 


: 20 

function of X/p increases rapidly with 
X/p, so that the decrease in « is not 10 
: : : : 2000 4000 5000 8000 i0000 
linear with 5 in the exponential term. i 
None the less, it is there. Hence, as Figure 5. Sparking field strength divided by 
a) or pd increases, 1/x, the mean distance pressure in ky/cm. plotted against the 
for ionization in the avalanche, in- product p64 in millimetres times centimetres 

Thi ih jonificance obtained from various sparking data below 
ee: Ne Arne tt the critical gap length 4 , indicating the 
on Mecek’s theory as formulated. decline in sparking field strength with in- 
However, as stated, this theory regards creasing gap length at constant pressure 


owing to increased ionization. 


only the space charge field which, from 
the slow increase of channel radius due 
- to diffusion, should always be adequate under equation (8). It does not regard 
the density of photo-ionization (Loeb and Meek 1941, pp. 49, Tite 10) AS. 
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however, 1/x increases, that important density decreases. ‘Thus the photo- 
ionization in advance of a streamer tip can fall below the critical value if 6 is in- 
creased at constant pressure. This will prohibit the spark occurring at the values 
of X predicted by Meek’s theory. On the other hand, if path length, and pressure 
or gas density, i.e. pS, increase with increasing p, Meek’s theory should apply. 
It is thus expected to be that in long sparks and in lightning at atmospheric 
pressure at some critical path length 59, Meek’s equation should break down 
(Loeb and Meek 1941, p. 79). What this critical length is we do not know, but 
Edwards and Smee (1938), with alternating 60c/s. potentials, observe that at 
gap lengths between 15 and 20 cm. at 760 mm. Hg the (V,, p3) curves pass through 
a point of inflection, increase and are different in trend for higher values of V. 
This is what might be expected on the anticipated breakdown of Meek’s theory. 
For an avalanche crossing a gap greater than some critical length 6) can no longer 
produce the photo-ionization necessary, and sparks will not pass until X, is in- 
creased to a point where 1/« has the values required to give the needed photo- 
ionization. ‘Thus, beyond 89, X,/p cannot decrease further and X, will remain 
at a constant value X5, corresponding to the limiting field at 6). Sparking will 
then occur only at or above a field X;,, and V, will be given by X55 above 4, 
irrespective of 5. The sparking potential curve given by Meek’s equation will 
follow the full curve of figure 6. At 5, Meek’s equation no longer holds and sparks 
will pass at a constant field 
X;, such that V,=X;,6, 
which increases linearly 
with 6 as shown by the 
straight dashed line which is 
extrapolated back through 
the origin; the actual 


TOWNSEND | 


transition will not be sharp. REGIME | 
This, then, indicates the A 
trend observed by Edwards 
and Smee. 


What actually happens Figure 6. Sparking potential plotting against gap length at 
in this region is that when constant pressure, the full curve indicating the actual 

= ~ G > sparking potentials to be expected through the whole 
V, nie X5,0 an avalanche from régime. Note the transition at 9, which extrapolates back 
the cathode proceeds a dis- through the origin. The dot-dash curve between the two 
tance do where its streamer indicates the probable shape of observed sparking 
starts; at the same time, potential in the transition region. ‘The Townsend régime 


; : holds f * i i - 
with cesuléantacld diate olds for the shorter gaps, Meek’s equation for the inter 
mediate, and at longer gap lengths the sparking potential 


tion, its electrons pass to the is a constant field strength times the gap length. 


anode, initiating at a dis- 
tance 5) or less a new streamer directed towards the cathode, and soon. At path 
lengths 5~20 cm. the streamer channel will be some 0-2 to 0:5 cm. in radius, which 
facilitates ionization according to Raether’s calculations. The statistical fluctu- 
ations in photo-ionization will smooth the range of critical potentials, fields and 
path lengths about 6, so that the transition is not sharp. The sparking potential 
curves should then appear as in figure 6. 

Considerations of the effect of pressure indicate that , will be a critical function 
of pressure and rough calculations by Loeb (Loeb and Meek 1941, p. 81) placed 
dy at 38cm. at half an atmosphere. L. H. Fisher (1947), in Loeb’s laboratory, 
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attempted to observe this in a carefully controlled chamber during the War. His 
investigation failed owing to field distortion and corona just below the critical gap 
length. It is urgent that this study be undertaken in some laboratory where high 
potentials and adequately large spark discharge chambers are available. In 
carrying out the study extraordinary care must be taken to avoid field distortion. 
The urgency of the problem will be understood when it is realized that it is the 
condition that determines flash-over for long static sparks and establishes lightning 
discharge breakdown fields. 
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ABSTRACT. Phosphorescence and thermoluminescence emission from photoconducting 
impurity activated phosphors have been satisfactorily explained by the storage of electrons, 
freed from luminescence centres or other atoms of the solid, in metastable energy levels 
known as electron traps. Electrons escaping from these traps give rise to emission when 
they recombine with luminescence centres but there is a probability that they may be re- 
trapped in empty electron traps before their final recombination with centres. The present 
theoretical and experimental studies attempt to determine the extent to which retrapping 
does occur and what effects it will have in modifying the phosphorescence and thermo- 
luminescence characteristics. Theoretical treatment shows that there are marked differ- 
ences in these characteristics for conditions when the retrapping process is present and for 
those when it is negligible. Experimental investigations of the characteristics of specimens 
of zinc sulphide, zinc silicate and strontium silicate phosphors indicate that, except under 
special conditions, retrapping of electrons is negligible. These results together with other 
work can be explained theoretically if it is assumed that electron traps operative in the 
luminescence process are spatially associated with the immediate neighbourhood of the 
luminescence centres formed by activating impurities. This new concept is also supported. 
by the relations found between the luminescence characteristics and the dielectric changes 
in phosphors of the zinc sulphide type. 


~ 


§1. INTRODUCTION 


T is now generally accepted that the phosphorescence and thermoluminescence 
characteristics of such phosphors as zinc sulphide and zinc silicate activated by 
suitable specific impurities are governed by the storage mechanism of electron 

traps which are thermally metastable. Previous workers have developed theoreti- 
cal models of a luminescent solid which interpret many of the luminescence 
properties in terms of the trapping of electrons (Johnson et al. 1939). Experi- 
mental and theoretical studies of a more quantitative nature have been carried out 
by Randall and Wilkins (1945) and by Garlick and Wilkins (1945). Recently 
their hypotheses have been extended to explain the origin of the relatively large 
dielectric changes occurring in zinc sulphide phosphors during the emission of 
luminescence (Garlick and Gibson 1947). 

The theoretical concept of a crystalline phosphor now generally accepted is. 
that based on the “collective electron” model of Bloch which has been developed 
by various workers (see Mott and Gurney 1940). The allowed energies for elec- 
trons in such a phosphor consist of bands of energy states separated by forbidden 
energy bands as shown in figure 1 for the outermost electrons responsible for 
luminescence processes. Impurities, lattice defects or other perturbations of the 
ideal crystal lattice can give rise to discrete energy levels like those of an isolated 
atom and these may lie in the forbidden regions as shown. While impurity atoms, 
ions or larger complexes can function as luminescence emission centres the discrete 
levels of more minor disturbances can provide electron trapping levels just below 
the conduction levels of the system. These traps have been assumed to be those 
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responsible for the phenomena of phos- te 
phorescence and thermoluminescence. Goi 
: ‘ = a a 
It is also generally assumed that in Se ee 
i =f aS is 


phosphors of the zinc sulphide type the 


electron traps are separate entities from 
the luminescence centres in contrast to E| ‘ os i) 
the case of potassium chloride activated 
by thallium where the trapping states are a 
within the thallium centres. Photo- Oy on Ge Se 
Pe : , ; ay F / 
conductivity present in zinc sulphide Uy Y oe Th 
phosphors supports this distinction. 
The mechanism of electron trapping - 
has been investigated in a relatively Figure 1. Schematic electron energy diagram 
cas eee Bey for an impurity activated crystalli 
quantitative way by Randall and Wilkins oNemetey ees ee eee 
; tee f phosphor. 
(1945), and from their experiments they Riese , 
ig is C. Conduction energy band. 
have postulated that, although the lumin- Re aed cneroe band: 
escence emission is due to recombination H. Positive hole. 
of excited electrons with empty lumin- L. Luminescence centre. 
¢ TT’. Electron traps. 


escence centres, yet the experimental evi- ere ee 
dence favours assumption of a monomole- hv’, Emitted quantum. 
cular process. As they have pointed out, 

this assumption is possible in explanations of phosphorescence and thermo- 
luminescence, which are determined then only by the rate of escape of trapped 
electrons and their return to empty centres, if there 1s no retrapping of escaping 
electrons before their recombination with centres. Using the schematic model of 
figure 1 it is evident that the degree of retrapping will depend at any instant on the 
relative numbers of empty electron traps and empty luminescence centres and on 
their relative capture cross sections. Randall and Wilkins have shown that these 
capture cross sections are approximately equal and for normal phosphorescence 
conditions the numbers of empty electron traps and empty luminescence centres 
are about equal. ‘Thus if electron traps and centres are independent entities the 
absence of retrapping cannot be explained by the model of figure 1. 

Some of the conditions under which retrapping is very likely have been given by 
Garlick and Wilkins (1945) and depend on a high proportion of empty traps 
compared with the number of empty centres. ‘These conditions occur at high 
temperatures where all traps are effectively empty, at the beginning of excitation of 
a previously unexcited phosphor and at very long decay times. ‘The present 
studies are concerned with the investigation of the occurrence of retrapping under 
normal conditions as defined above and attempt to provide from experimental data 
a more correct theoretical model than previous ones for phosphors of the sulphide 


and silicate classes. 
A note on the following studies has already been published (Garlick and 


Gibson 1946). 


§2. THEORETICAL STUDIES 


For the sake of completeness it is necessary to include and enlarge upon the 
theoretical developments of Randall and Wilkins (1945) following their assumption 
that the retrapping of electrons does not occur to any significant extent except for 
special conditions. Ifa phosphor contains only traps of one depth E and at any 
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time these contain n electrons, then the luminescence emission intensity / is 
determined by the rate at which electrons escape from traps, and the process is 
represented by the following equations: 

[= —dn/dt=nsexp(—E£/RT), 79 > ee (1) 
where sexp(—E/RT) is the probability per unit time that an electron escapes, 
s being a constant of the order of 108*!sec-! for most phosphors, E being the 
electron trap depth, k being Boltzmann’s constant and 7 the absolute temperature. 
If the equations (1) are solved for a uniform rate of rise of temperature after 
excitation of the phosphor at alow temperature the variation of thermoluminescence 
intensity with temperature is given by the expression 


P sexp(—E/kT 
[=ngsexp(—E/RT) exp | sexpert, ea (2) 
0 


where f is the rate of warming in deg/sec. Randall and Wilkins (1945) have dis- 
cussed the effect of the warming rate on the thermoluminescence variation. ‘The 
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Figure 2(a). Variation of thermoluminescence with temperature for phosphors with single trap 


depths after excitation at 90° K. Curves calculated theoretically assuming retrapping to be 
absent. Warming rate=2°5°/sec. 


(a) E=0:3 ev.:s=10% sect. (d) E=0-4 ev.:s=10" sec. ~ (f) E=0-8 ev. :s=10° sec-}. 
(6). E=0-4 ev. :s=10° sec—. (e) E=0-6 ev.:s=10*% sec.  (g) E=0-4 ev. :s=108 sec—. 
(c) H=0-4 ev. :s=108 sec. 


graphical form of equation (2) is shown in figure 2(a) for various values of s and 
FE, my and B remaining fixed (8 =2-5°/sec.). The chief characteristics of the 
curves of figure 2(a) for the variation of thermoluminescence with temperature 
when E is single valued are enumerated as follows :— 


(a) For given values of s and my) and £ the temperature at which maximum 
emission occurs is proportional to the trap depth E. 


(5) For given values of E and mp) the temperature of maximum emission varies 
with s and with 8. For a change in s value from 10%sec-! to 107sec when 
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E=0-4ey. the position of the emission peak moves to higher temperatures by 
about 50°k. A hundred fold increase in the warming rate will cause the same shift. 


(c) In all cases the area under the emission—temperature curve is proportional 
to the number of electrons initially trapped before warming begins, that is to mp. 
However, the actual shape of the emission curve is independent of 1, which also 
means that the height of the curve at any point is proportional to mp. 


(d) The initial rise of the thermoluminescence curve is exponential with 

temperature following the simple relation 

De igS OX LHe LZ) Poke a PA SS he lat (3) 
This affords a means of determining F from experimental curves of the same form 
when s is not known. 

By means of the above characteristics the experimental results from studies of 
phosphor specimens may be tested to see whether the phosphors have only one 
trap depth or a distribution of trap depths according to the above hypothesis. 

The above theoretical developments rest on the assumption that retrapping 
is a negligible process in the phosphor. It is essential however to consider the 
modification of this theory when an electron escaping from a trap has the same 
probability of being retrapped as it has of recombining with an empty luminescence 
centre. ‘There is evidence from previous studies that the two probabilities are 
about equal (Randall and Wilkins 1945, p. 405). Consider a phosphor containing 
a total number of electron traps N of which 7 are filled by electrons at any instant. 
Then there will be (N—n) empty traps and m empty luminescence centres 
previously vacated by the trapped electrons (the number of electrons in the con- 
duction band at any instant is assumed to be small compared with m during the 
phosphorescence and thermoluminescence). The probability that an escaping 
electron will recombine with an empty luminescence centre and not be retrapped is 
given by n/{(N—n)+n}=n/N and the luminescence process will be given by 


I= —dnjdt=(n2/N)sexp(—-EJRT) ses (4) 


with the notation used above. ‘These equations should be compared with the 
previous equations (2). 

The decay of phosphorescence at a fixed temperature is given then by solution 
of equations (4) and is represented by the equation 


[=ns exp (—E/kT) / N{1-+(n/N)stexp(—E/RT)2.  ...... (5) 


The decay is hyperbolic in form showing that retrapping if present can produce a 
marked modification of the exponential decay given by equations (1) obtained by 
neglecting its effect. It will be seen from equation (5) that the form of the decay of 
phosphorescence will be dependent on m, and thus on the conditions of filling of the 
traps by excitation at the commencement of decay. 

If equations (4) are solved to give the thermoluminescence variation with 
temperature when the phosphor is warmed at a uniform rate after excitation at a 
low temperature, then the resulting variation is given by the equation 


Geren p(s ~HikT)[n{i+ Bf SPEER gpl Peace (6) 


This equation corresponds to equation (2) and its graphical forms for various values 
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of E and ny when s and f remain fixed are given in figure 2 (4). The main char- 
acteristics of thermoluminescence-temperature variation when retrapping occurs 
are enumerated as follows :— 

(a) For fixed values of s, ) and f the temperature at which the emission reaches 
its maximum is proportional to the electron trap depth Z. For saturation of the 
traps ()=.N) and the same values of s the equations (2) and (6) give maxima at the 
same temperature. 

(b) For given values of E and ny the temperature position of the emission 
maximum varies with s and f in much the same way as for the case when retrapping 


is absent. 
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Figure 2(b). Variation of thermoluminescence with temperature for phosphors with single trap 
depths after excitation at 90° k. Curves calculated theoretically assuming retrapping to be 
present. Warming rate=2:5°/sec. s=10%sec-!; N=total number of traps of depth E ; 


Ny=number filled. 


(a) E=0-4 ev. : m=N. (ec) B=0-6 ev. 275=0°5N- (@)) E==0°8 ev; 2 725—0-25.Ne 
(6) E=0-6 ev. : m=N. (CPE =0:3 ev-7)— 0-5 Ne (j) #=0-4 ev. : m=0-1N. 
(c) -£=0°8 ev. 2, —N- (g) E=0-4 ev. : 295=0:25N. (k) H=0-6 ev. : n7=0-1N. 
(d) E=0-4 ev. : m9=0°'5N. (h) E=0-6 ev. 2 m)=0-25N. (2) £=0°8 ev. : m9=0-1N. 


(c) In all cases the area under the emission—-temperature curve is proportional 
to the number of electrons (mp) initially trapped. However, the actual shape of 
the curve is dependent on mas shown by equation (6). ‘This is essentially different 
from the case when retrapping is negligible. 


(d) 'The initial rise of luminescence before peak emission is reached follows 

the simple relation 

l=(n?/N)s exp (.-E/RT). © Ae ee (7) 
‘Comparing this equation with equation (3) it is seen that the power of the 7) term 
is different in the two equations, the exponential term being identical. This 
‘provides a further means of experimental correlation with theory. 

The above characteristics obtained by the two different assumptions of re- 
trapping and no retrapping show the marked effect which retrapping may have on 
the phosphor characteristics. For most phosphors which have complex electron 
trap distributions the above treatment must be extended to traps of many depths. 
‘This is relatively simple if retrapping is neglected as shown by Randall and Wilkins 
for particular trap distributions. For example if the number of traps decreases 
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exponentially with depth (i.e. N oce~*”) then the phosphorescence decay 
for initial saturation of the traps involved is of the form J oct~@"2+» where ¢ is 
the decay time. he trap distributions may be obtained by means of thermo- 
luminescence experiments and in the above case the theory verified by correlation 
of the results with phosphorescence measurements at different fixed temperatures, 
since the index of ¢ should vary linearly with the temperature. This particular 
case shows that the hyperbolic decays characteristic of zinc sulphide phosphors 
can be explained by the sum of exponential terms whose constants are known from 
experiments. ‘The derivation of phosphorescence characteristics for phosphors 
with complex distributions when retrapping is present becomes difficult as 
electrons escaping from traps of one depth may be retrapped in traps of different 
depths. Solutions have not yet been obtained for these conditions. 

It is necessary to consider the way in which the trapping processes are affected 
by the conditions of excitation of the phosphor both with and without the presence 
of retrapping of escaping electrons. Asa simple case that of steady excitation of a 
phosphor at fixed temperatures will be considered. 


(1) The filling of traps by constant excitation with no retrapping 
The considerations below make a fundamental change in the theoretical model 
of figure 1. ‘They assume that each luminescence centre has an electron which on 
excitation is raised into electron traps which are in the immediate neighbourhood 
of the centre or are metastable trapping levels within the centre. Then if ultra- 
violet light of intensity J excites aJ electrons per second and in equilibrium the 
number of electrons trapped is 7) in the N available traps of depth F then the 
equilibrium is represented by. 
dn,/di=0=abJ(N —ny)—nmysexp(—E/RT) -...... (9) 
where 6 is the probability of electron capture by a trap. ‘This gives the relation 
between m, and N for different intensities of ultra-violet radiation, temperatures 
and trap depths as 
y= NAPE (siabd)expCWiRE)}- 9) © Ranati. (10) 
This equation may be tested experimentally by both thermoluminescence and 
phosphorescence experiments. It may be extended to apply to complex trap 
distributions as has been done in the theoretical treatment of dielectric changes 
due to electron traps in phosphors (Garlick and Gibson 1947). For the purposes 
of the present experimental studies equation (10) is rearranged to show the simple 
relation between nm, and J for fixed trap depth and temperature 


CN ing) —l=(s/abd yexp(— ERT). 2 5 aca. (11) 


(ii) The filling of electron traps when retrapping 1s present 
With the above notation and assuming retrapping to be present as in the deri- 

vation of equations (4), (5) and (6), the filling of the N available traps by constant 
excitation J at a fixed temperature may be determined. If aJ electrons are raised 
into the conduction band per second by the excitation and in equilibrium mp 
electrons are in the band, 1, being in traps, then (7) + my) centres are empty and the 
equilibrium equations are 

dm, [dt =0 = aJ — Bim (mtg + 19), 

dny|dt =0 = bm (N —ny) —nos exp ( —E/RT), 
where b is the probability of capture of an electron by an empty trap and f the 
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probability of electron recombination with an empty luminescence centre. 
Solution of these equations shows that below saturation level for the traps (i.e. 
N>n) both m, and ny are root functions of the excitation intensity J. Thus when 
retrapping is present the simple relation between m, and J as in equations (10) and 
(11) is not found and the examination of trap filling thus provides a means of 
testing the above relations for validity. 

The above theoretical developments have been compared with the results of 
experimental studies on typical phosphors as described below and these compari- 
sons provide evidence for the relative importance of retrapping. 


§3. EXPERIMENTAL STUDIES AND RESULTS 


(i) Experimental methods 

The experiments consisted essentially of measurements of thermoluminescence 
and phosphorescence characteristics of selected specimens under many different 
but known conditions of excitation and temperature. For all these experiments. 
the phosphor was mounted 7” vacuo in a thin layer on a rhodium-plated copper 
surface forming part of a metal Dewar vessel described in a previous paper (Garlick 
and Wilkins 1945). The specimens could be excited by monochromatic radiation 
of wavelengths 2537 a. or 3653 a. by means of amercury arcsource. The methods 
of measuring thermoluminescence curves and phosphorescence curves. have also. 
been previously described (Randall and Wilkins 1945). The rate of warming in 
all thermoluminescence experiments was kept fixed at 2:5°/sec. 


(ii) Experimental results and theoretical correlations 

For the experimental investigations three phosphor specimens have been 
chosen with narrow trap distributions giving a peak in their saturated thermo- 
luminescence curves which are shown in figures 3 (a), 3 (6) and 3 (c). The specimen 
corresponding to figure 3 (a) is a zinc sulphide phosphor activated by copper and 
excited by 3653. radiation. It has two peaks in.its thermoluminescence curve 
but in experimental studies use was made of the peak due to deeper traps only. 
Figure 3 (b) is the thermoluminescence curve for a zinc silicate phosphor activated 
by manganese and excited by 2537 a. radiation having a single peak in its emission 
curve. In both figures 3(a) and 3(6) the theoretical curves obtained by substi- 
tution of suitable values of s and £ in equation (6) are shown by broken curves. 
The experimental curves show remarkably good agreement with the calculated 
curves for which retrapping has been assumed. However, on the basis of “‘ no 
trapping”? these phosphors would have a distribution of trap depths. Both 
specimens are typical of their class and are photoconductors. The third specimen 
studied is a strontium silicate phosphor with europium impurity as activator. 
Its thermoluminescence curve has one main peak as shown in figure 3 (c) and this 
represents a single trap depth as shown by the agreement with the theoretical curve, 
the broken curve of figure 3 (c), derived from equation (2) assuming no retrapping. 
This phosphor is excited by ultra-violet light of both wavelengths used but 
thermoluminescence and the corresponding phosphorescence of long duration 
only occurs after excitation by 2537. radiation. The unique characteristics of 
the luminescence of this phosphor are described and discussed later. The follow- 
ing experiments attempt to correlate the characteristics of each phosphor with one 
or other of the theoretical hypotheses. By this means the extent to which 
retrapping occurs may be defined. 
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Figure 3(a). ZnS—Cu phosphor. Thermoluminescence characteristics after excitation 
at 90° x. 


(a) Theoretical curve assuming retrapping to occur. (6) Experimental curve. 
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Figure 3 (6). ZngSiO,—Mn phosphor. Thermoluminescence characteristics atter excitation 
at 90° kK. 


(a) Theoretical curve assuming retrapping to occur. (6) Experimental curve. 
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SrSiO,—Eu phosphor. Thermoluminescence characteristics after excitation 
at 90° K, 


(a) Theoretical curve assuming retrapping to be absent, (b) Experimental curve. 


Figure 3 (c). 
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THE DECAY OF PHOSPHORESCENCE 
AT FIXED TEMPERATURE 


§ 4. 


(i) Copper-activated zinc sulphide 


The decay of phosphorescence of the zinc sulphide phosphor has been mea- 
sured at several different fixed temperatures. ‘The decay as measured at room 
temperature is shown in figure 4, curve a, the coordinates of the graph being 
logarithmic to show the hyperbolic form of the decay, at longer times. At room 
temperature for these longer times the intensity of phosphorescence is inversely 
proportional to the 2:2 power of the decay time. This power is found to vary 
linearly with temperature as shown by the measured values at different temper- 
atures given in table 1. 


Table 1. Power law indices for the phosphorescence decay of ZnS—Cu phosphor 
Temperature (° K.) 265 292 318 355: 
Power law index («kT-+1) 2:07 2D Dy, 25 
x (ev-) 47 47-7 48 49 
% from thermoluminescence 49 aN 49 ao) 


measurements (ev?) 


These experimental results indicate that the phosphor does not exhibit the 
phosphorescence characteristics predicted by equation (5) although its thermo- 
luminescence curve may be interpreted by equation (6) since the final asymptotic 
form of the decay varies with temperature. SS eee 
However, if the thermoluminescence 
curve and phosphorescence decay curves 
are interpreted according “to the 
hypothesis which neglects retrapping, 
then the linear variation of the power 
law index is expected from equation (8), 
and from the decay curves « of the equa- 
tion may be obtained. ‘The constant « 
may also be obtained from the exponential 
form of the thermoluminescence curve 
at its higher temperature end given in 
figure3(a). Values of« derived fromthe 
decay curves are given in the third line 
of table 1 and the value from the curve of 
figure 3 (a) in the fourth line of the table. emer 
In the latter case the electron trap con- 
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stant s had the value 10° sec~! found by 
other experiments as previously described 
(Randall and Wilkins 1945, Garlick and 
Gibson 1947). Good correlation is ob- 


Figure 4. Decay of phosphorescence with time 
for zinc sulphide and zinc silicate phosphors 
at room temperature (291° k.). 


(a) Copper activated zinc sulphide specimen. 


(6) Manganese activated zinc silicate specimen. 
tained between the differently derived 


values of a. ‘hus experimental evidence from decay measurements shows 
that the deeper traps in this phosphor, while giving rise to a relatively narrow 
thermoluminescence-temperature curve, have a distributed range of depths and 


that agreement of experimental results with theory necessitates the assumption 
that retrapping is negligible. 
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(11) Manganese-activated zinc silicate 


‘The decay of phosphorescence of this specimen at room temperature is shown 
by curve (5) of figure 4. It is of the similar hyperbolic form to curve (a), at longer 
decay times its power law index being 1:7. By use of equation (8) this gives the 
trap distribution constant « as 28-Qev-!. Derivation from the thermolumin- 
escence curve of figure 3(b) gives « as 25ev-!. Again agreement with theory is 
good when retrapping is neglected. The value of the power law index is found to 
increase linearly with temperature as in the case of the zinc sulphide phosphor. 
These facts provide the data supporting the theoretical hypothesis that retrapping 
is not important. 


(111) The phosphorescence of europium-activated strontium silicate 


This phosphor may be excited to give fluorescence and phosphorescence by 
radiation of either 3653 a. or 25374. However, there is no long duration phos- 
phorescence or any thermoluminescence when the longer wavelength radiation is 
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‘igure 5 (a). Phosphorescence decay of Figure 5 (b). Phosphorescence decay of SrSiO;—Eu 
SrSiO;—Eu phosphor at various phosphor at various temperatures after exci- 
temperatures after excitation by tation by 2537. radiation. 

3650 A. radiation. ° (a) 387°K. (a) 340° x. 
(a) 476°K. (d) 192°. (b) 373°. (e) 320° x. 
(6) 398° k. (e) 90° k. (c) 358° K. 

(c) 294° kK. 


used. The phosphorescence decay with time has been measured after excitation 
by each of these wavelengths of ultra-violet radiation. Figure 5 shows the results 
obtained. Figure 5 (a) shows the form of the phosphorescence decay after excita- 
tion by radiation of 3635. wavelength at different temperatures from 90°K. to 
476°«. ‘The initial decay is rapid, occupying a small fraction of asecond. ‘The 
subsequent decay has an exponential form, the decay rate being independent of 
temperature. ‘The half-value period of the decay at different temperatures from 
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90° K. to 476° K. is included in the figure, its mean experimental value being 2:25 sec. 
Since this decay rate is independent of temperature it cannot be due to thermally 
metastable electron energy levels. This is also supported by the absence of thermo- 
luminescence and photoconductivity when the 3650. is used. ‘Thus the phos- 
phorescence must be due to a forbidden optical transition within the europium 
luminescence centre, the life time of the state being abnormally large. Similar 
decay characteristics have been observed in other studies of the phosphorescence 
of adsorbed dyestuffs made according to the method of T'ravnicek (1937). 

When radiation of wavelength 2537 a. is used to excite the specimen the latter 
exhibits both photoconductivity and thermoluminescence. Its phosphorescence 
decay at various temperatures is given in figure 5(b). The initial decay is due to 
the same two components which are present when 3653 A. irradiation is used, but 
the decay includes also phosphorescence of long duration decaying exponentially 
with time over the range measured but having a decay rate markedly dependent 
ontemperature. This decay must be due to electron traps with a single depth or at 
the most a very narrow spread of trap depths as indicated by the thermolumin- 
escence curve of figure 3(c). The absence of a hyperbolic decay due to these 
electron traps provides evidence that the hypothesis of Randall and Wilkins 
applies in this case. The phosphorescence decay is thus of the form :— 


awa) Se ee ee PP et (14) 
From the change of decay rate the value of E has been calculated and is found to be 


about 0:32ev. This value is not in agreement with those obtained from thermo- 
luminescence measurements and is discussed in the section below. 


I=cexp(-—ste 


§5. STUDIES OF THERMOLUMINESCENGE CHARACTERISTICS 


The measurement of the variation of the thermoluminescence characteristics 
of a phosphor with different intensities of initial excitation provides a means of 
determining how the electron traps are filled for the different excitation conditions. 
Each phosphor specimen is therefore excited at liquid air temperature with various 
intensities of ultra-violet light of suitable wavelength for a sufficiently long time so. 
that equilibrium conditions are reached. Then after excitation is removed the 
thermoluminescence-temperature curves are measured. ‘These curves are given |] 
in figures 6 (a), (b) and (c). ‘The area under each curve represents the number of |} 
electrons initially trapped, denoted by mp, under the equilibrium excitation con-_ 
ditions at 90°K. It will be seen from these curves that the temperature at which 
maximum emission occurs does not alter appreciably as the excitation intensity is 
varied. 

In order to compare these experimental results with the theoretical predictions] 
represented by equations (11) and (14) the area under each curve is taken as 7, | 
the number of electrons initially trapped, and the area under the curve for very high 
intensities as the saturation value N. Then for each specimen the experimentally |} 
derived values of ((N/n )—1) are plotted against the inverse of the excitation | 
intensity J. ‘These plots are given in figures 7 (a), (b) and (c) for the three speci-_ ] 
mens respectively. From inspection it is found that the experimental results. |) 
follow the relations given by equation (11), i.e. retrapping assumed to be absent. | 

A further correlation of experiment with theory can be made by considering | 
the portion of the thermoluminescence curve at some temperature lower than that} 
at which peak emission occurs. As shown by the equations (3) and (7) the shape | | 
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(a) Zinc sulphide—Copper-activated phosphor. 
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Figure 6. Variation of thermoluminescence of sulphide and silicate phosphors with intensity of 
excitation (J) at 90° k. 
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(a) Zinc sulphide—Copper- (6b) Zine — silicate—Man- (c) Strontium sulphide— 
activated phosphor. ganese-activated phos- Europium-activated 
phor. phosphor. 


Figure 7. Variation in the number of trapped electrons (m9) in the number of available traps (NY 
with excitation intensity (J) at 90° K. 


of this portion of any curve will depend upon the degree of retrapping present. 
From this shape the value of trap depth may also be obtained for each specimen. 
Figures 8 (a), (6) and (c) show how the height of the initial portion of the curve at 
some selected fixed temperature varies with the number of electrons initially 
trapped, as given by the total area under the curve. Results for each of the three 
specimens show that equation (3) is adhered to. Calculations of the trap depth 
from the curve shape have been made for the strontium silicate specimen where the 
depth is almost single-valued. ‘The depth thus obtained is 0-68ev. By using this 
value and the temperature of peak emission the value of the constants may be found. 
This is calculated to be 10° sec~! which is of the usual order for silicate and sulphide 
phosphors (Randall and Wilkins 1945, Garlick and Gibson 1947). Again by 
using this value of the constant the value of E may be calculated from the 
phorescence decay curves of figure 5(b) since the decay rate is given by 
sexp(—E/RT). ‘The values thus obtained are collected in table 2. 
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Figure 8. Variation in height of initial portion of thermoluminescence curve with the number of 
electrons (mp) initially trapped by excitation at 90° K. 
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Table 2. Values of electron trap depth for specimen of europium-activated 


strontium silicate obtained from phosphorescence characteristics with 
constant s=10°sec"!, 


Temperature (° k.) 320 340 358 373 387 
Decay rate s exp (—E/kT) 0:00935 0-013 0-021 0-031 0-05 
Electron trap depth E (ev.) 0-74 0:74 0:75 0:77 0:77 


As shown by the phsophorescence curves the values of / given in table 1, 
while agreeing approximately with that obtained from the thermoluminescence 
curve, do not correlate with the value obtained from the change of decay rate with 
temperature. This may be due to a small spread of E values for this specimen 
similar to that found by Randall and Wilkins (1945) for potassium chloride 
activated by thallium. Otherwise if F is considered single valued the value of the 
constant s must be of the order of 10*sec™ according to decay rate changes, which 
does not agree with the value obtained from the thermoluminescence measure- 
ments. 


(ii) Changes in thermoluminescence characteristics during phosphorescence decay 
and after pulse excitation 


The following experiments provide further data to distinguish between re- 
trapping and no-retrapping processes in the phosphor specimens and follow the 
method used in earlier studies (Garlick and Wilkins 1948). Each specimen is 
excited at liquid oxygen temperature by a pulse of ultra-violet light whose duration 
canbe varied. After such excitation the variation of the thermoluminescence with 
temperature is measured. The curves obtained are shown in figures 9(a), (c) 
and (e) for the three phosphor specimens respectively. ‘The specimens are also 
excited at room temperature (291°k.) and allowed to phosphoresce after removal 
of excitation for various selected times after which periods the specimen ts suddenly 
cooled and then warmed at the unifrom rate of 2-5°/sec., its thermoluminescence 
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Figure 9(a). ZnS—Cu _ phosphor. Figure 9 (0). ZnS—Cu_ Phosphor. 
Variation of thermoluminescence Variation of thermoluminescence 
with duration of pulse excitation after different times of phos- 
at 90° kK. phorescence decay at 291° kK. 
(a) Steady excitation. (a) Excitation at 90° Kk. 
(b) Medium length pulse. (6) Decay for 1 sec. at 291° k. 
(c) Short pulse. (c) Decay for 2 sec. at 291° kK. 
(d) Very short pulse. (d) Decay for 3 sec. at 291° k. 


(e) Decay for 10 sec. at 291° k. 
(f) Decay for 20 sec. at 291° k. 
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Figure 9 (c). Zn ,SiO,—Mn phosphor. Figure 9 (d). Zn ,SiO,—Mn phosphor. 
Variation of thermoluminescence Variation of thermoluminescence 
with duration of pulse excitation after different times of phos- 
at 90°K. phorescence decay at 291° kK. 
(a) Steady excitation. (a) Excitation at 90° kK. 
(b) Long pulse. (b) Decay for 1 sec. at 291° k. 
(c) Medium pulse. (c) Decay for 4 sec. at 291° kK. 
(d) Short pulse. (d) Decay for 10 sec. at 291° kK. 
(e) Short pulse. (e) Decay for 20 sec. at 291° K. 
(f) Very short pulse. (f) Decay for 60 sec. at 291° x. 
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Figure 9(e). SrSiO;—Eu phosphor. Figure 9(f). SrSiO3—Eu phosphor. 

Variation of thermoluminescence Variation of thermoluminescence 
with duration of pulse excitation after different times of phos- 
at 90° K. phorescence decay at 291° kK. 

(a) Steady excitation. (a) Excitation at 90° K. 

(6) Medium pulse. (b) Decay for 30 min. at 291° x. 

(c) Short pulse. (c) Decay for 1-5 hr. at 291° x. 


(d) Decay for 16 hr. at 291° x. 


temperature variation being measured during warming. The curves so obtained 
are given in figures 9(b), (d) and (f) and represent the distribution of electrons in 
traps at different stages of phosphorescence decay. Cross comparison between 
curves 9 (a), and (4), (9c) and (d) and (9 e) and (f) respectively shows that change in 
duration of pulse excitation affects the number of trapped electrons but the shape 
of the thermoluminescence curves and their maxima positions are not markedly 
altered. However, phosphorescence decay does cause a shift of the maximum 
position of these curves at higher temperatures as the decay proceeds for the 
specimens of zinc sulphide and zinc silicate. The strontium silicate phosphor 
characteristics in curves 9(e) and 9(f) show no change in the maximum position 
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with the number of trapped electrons either during phosphorescence or due to 
pulse excitation. 

These results confirm the hypothesis that retrapping is a negligible process in 
these phosphors. As stated in the theoretical section above the thermo- 
luminescence curve shape for a single trap depth when retrapping is present 
depends on the number of trapped electrons (figure 2 (b)) but is independent of the 
process of filling or emptying, whether by pulse excitation or phosphorescence 
decay. ‘Thus for a given area under the curve the shape should be fixed. As 
shown by comparison of curves (a) and (b), (c) and (d), (e) and (f) respectively 
of figure 9 the experimental results do not agree with this prediction but are in 
agreement with the supposition that a distribution of trap depths is present and 
that retrapping is negligible. This means that pulse excitation of the previously 
unexcited phosphor will fill traps of all different depths with equal probabilities, 
the number filled thus depending only on the number available. Thus the shape 
of the thermoluminescence curve is little affected. However, when phosphor- 
decay proceeds the shallower traps are emptied leaving only the deeper ones filled 
with the consequence observed in the curves of figure 9(b), (d) and (f), that is, 
a temperature shift of the emission maxima. 


§6. CONCLUSION 


The experimental results described above have been obtained from studies of 
three typical phosphor specimens of three different classes. In each case the 
specimen has been selected because it affords suitable thermoluminescence and 
phosphorescence properties to determine how far retrapping of electrons escaping 
from electron traps was a major part of the mechanism. ‘The evidence presented 
has shown that the characteristics measured are in agreement with the simple 
theoretical treatment originally due to Randall and Wilkins which assumes 
retrapping to be negligible. This could be understood in terms of the model of 
figure 1 if, as in the case of the thallium activated alkali halide phosphors, the 
electron traps are within the luminescence centres to which the phosphorescence 
processes would then be confined. However, the specimens all exhibit photo- 
conductivity which is associated with the thermoluminescence and phosphor- 
escence processes and therefore electrons must leave luminescence centres during 
the process. ‘This means that the luminescence process involves a recombination 
of electrons with luminescence centres after escaping from electron traps which 
are not in the luminescence centres themselves. Previous experiments due to 
Randall and Wilkins (1945) have given estimates of the capture cross sections of 
empty luminescence centres and electron traps which appear to be of the same order. 
Thus retrapping of escaping electrons from traps will be an integral part of the 
mechanism if traps are quite independent of luminescence centres. It is therefore 
more probable that the impurity centres while not containing the trapping energy 
states cause perturbations in their environment (i.e. in the surrounding local 
lattice) which give rise to the trapping states. If this is so then the absence of 
retrapping of electrons may be more readily understood since electrons escaping 
from traps adjacent to empty luminescence centres will have a higher probability of 
recombining with the centres than of moving to other empty traps and being recap- 
tured. It is not possible to provide a detailed treatment of this assumption as no 
studies of the complex potential configurations around impurity centres in a matrix 
lattice have been made. However, the postulate of association between traps and 
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luminescence centres is supported by other and independent studies of the effect 
of one or more different impurities and also of the preparation conditions (firing 
temperature, fluxes etc.) on the electron trap distributions of phosphors (Garlick 
1947). These experimental investigations show that introduction of specific 
impurities gives rise to particular electron trap configurations. Other studies of 
the dielectric changes in zine sulphide and other phosphors show good agreement 
between theory and experiment when retrapping is neglected (Garlick and Gibson 
1947), ; 

To summarize, the experimental and theoretical developments described 
above show that the simple electron trap theory due to Randall and Wilkins (1945) 
which neglects retrapping is in better agreement with experiment than that in 
which retrapping is assumed to be present. This evidence, combined with that 
from other investigations leads to a new concept of the origin and spatial location 
of electron traps in phosphors. They appear to be formed by the perturbanons 
due to the impurity activator in the matrix crystal lattice immediately surrounding 
the centre. 
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The Moving Coil Galvanometer as a Circuit Element 


By N. F. ASTBURY 


G.K.N. Research Laboratories, Wolverhampton 
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ABSTRACT. An account is given of a simplified calculus applicable to a moving-coil 
galvanometer in which use is made of the concept of motional impedance. The fundamental 
galvanometer equations are presented in a novel and very simple form. Examples of the 


use of the formulae are given, notably with reference to the ultimate “noise level” of a. 
galvanometer. 


$1. INTRODUCTION 
HE moving-coil galvanometer, one of the oldest and most familiar of our 
laboratory instruments, has recently found considerable application 
in various systems for direct-current amplification, amongst which may 


be noted those described by D. C. Gall (1942, 1945). It is hoped to publish 
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shortly descriptions of other systems developed by the present writer, but in 
the meantime it may be noted that the problem of p.c. amplification is one of 
very great technical importance and one which is most admirably met by the 
use of a moving-coil galvanometer in a preamplifying stage. In the course uf 
development work on problems of this type it occurred to the writer that the 
mathematics of a moving-coil galvanometer could be considerably simplified 
for many purposes of circuit analysis, and it is the purpose of this paper to show 
how the behaviour of a galvanometer as a circuit element can be succinctly ex- 
pressed in terms of its motional impedance, which can in turn be expressed in 
terms of experimentally determined constants. 

While the concept of motional impedance applied to galvanometers is by 
no means new, it is believed that the present treatment is novel. Butterworth 
(1912) showed that a galvanometer can be represented by a parallel combination 
of a pure inductance, a capacitance and resistance placed in series with a resistance 
equal to that of the instrument coil, and that in this equivalent circuit the current 
in the inductance branch is homologous with the angular deflection of the coil. 
Kennelly showed that the impedance circle diagram was applicable to a galvano- 
meter. A note on these points is given in the Dictionary of Applied Physics 
(Vol Li 2-976). 

It is proposed in the present paper to put forward the motional impedance 
concept and to illustrate it by applications to various circuit problems. 


§ 2. PRECIS OF STANDARD RESULTS 
It is convenient to start with a précis of well-known results for the galvanometer. 
Let x, B and y be respectively the moment of inertia of the moving system, 
the mechanical damping factor and the torque control constant of the galvano- 
meter, and let k be the product of the total area of the coil and the control field, 
i.e. the ‘“‘area-turns-field”’ constant. Then, if the galvanometer form part of 
a closed circuit of total resistance R (and of zero reactance) in which an E.M.F. é is 
acting, the current z and the angular deflection 6 of the moving system are related. 


to e by the torque and E.M.F. equations 


(oD? BD y= hp SS (221) 
Ri+kD6 = nde) 
where D stands for d/dt. 
Eliminating z, we have 
{aD?+(B+R/R)D+yjO=he/R, sae (2.3) 
from which we find S, the critical resistance, as 
Sekula) ef | 6 anes (2.4) 
where wy is the free frequency of the instrument and is given by 
Wy =(y/a)t(1 —B2/4ay)P(y/a)P. nee 2) 
The logarithmic decrement, A, is, if T=27/w, 
NAB RADE IZa) oreo one (2.6) 
and 
cel hig key a i ie eer (227) 


where A, is the value of A for the open circuit condition. We may define the: 


damping factor, 6, by 
De Ng TE Cp aatasale (2.8) 
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Finally, we may define the direct current sensitivity, u, by 
pekRiy 0) le eee (2.9) 


and then express all our coefficients in terms of the experimentally observable 
quantities S,w, and pw, obtaining 


Pe (2.10) 
B=LSIGTOy i) en ee (274.15) 

= 2S0/Wo7[i75. 4 a ae (2.12) 
y= 2S /i705,. = el ee (25h) 
R=2S/ ag. > Oe eee Qi) 


These results can be used to express Butterworth’s (1912) formulae in simple 
forms. If L, C and P are respectively the inductance, capacitance and shunt 
resistance of his equivalent circuit, we have 


Li=2Slaepit 8 4 te ee (2.15) 
C=1/2Sag2 0 ee (2.16) 
P25) -= 2S ee (2.17) 


§3. THE MOTIONAL IMPEDANCE 

We can now proceed to a discussion of the motional impedance of the galvano- 
‘meter, but we may note before so doing that it is proposed to regard the circuit 
in the static condition as non-reactive. There is no loss of generality in this, 
because the dynamic reactance of the instrument is always so much greater 
than the static value that the inclusion of the latter in the discussion would be 
a mere algebraic tour de force. . 

Eliminating 6 from equations (2.1) and (2.2), we find 


={R-+k?D)/(ob2+ BD+y)\ | ae oes (3.1) 
It is evident that the second term in the bracket represents the contribution 
‘made to the effective impedance of the circuit by the moving galvanometer 
coil. ‘This term we call the motional impedance (more precisely it is the motional 
impedance operator) and we denote it by J. Then 
l=PD/(aD*? + BD+y) ~ 9 ee (3.2) 
= 2S0,5D/{(D? +.0)57)+5a9D}, esse (3:3) 
‘using equations (2.4), (2.5), (2.8) and (3.1). 
From (2.1) and (3.2) we find 


f fe es ae 
O= 5 5 i= z lide Me (3.4) 
‘while from (3.1) and (3.4) we have 


tte, ahh ot! 1 MN 
=E-TER De TRO alielte,e "ele (3.5) 


Equations (3.4) and (3.5) present the equations of the moving-coil galvanometer 
in a novel and simple form which is of considerable value in circuit analysis. 
Equation (3.5) is of particular interest, because it holds the key to the problem 
of the behaviour of a galvanometer as a fluxmeter and shows at once that for 
flux measurement the condition is that the motional impedance shall dominate 
‘the static resistance of the circuit. 
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The expression J/(J+R) may be regarded as the “response operator”’ of 
the instrument and its circuit, and may conveniently be rewritten as a function 
of D, wo and o, where o=S/R may be called the “degree of damping” of the 
instrument. We have 

te as 2w aD 
T+R”~ D?farfayD@o+s) "oo * 

While we shall discuss later the application to particular forms of time- 
variation of e, the form taken by J for sinusoidal variations, which is useful in 
connection with vibration galvanometers, is interesting. If w be the frequency 
of e, then we have at once from (3.2) 


2]mywS 
(wo" —w?*) +]Jww,d , 


l= 


where j is the operator rotating through a right angle. This may be written 


285  2%Sz 


i 22462 22482 


where 2=(w )/w) —(w/wo), showing the existence of a motional resistance term, 
Ry, and a motional reactance term, X,, given by 


Ry =283/(22+32) ; Xy=2Se/(a?-+8%), 

It is easy to show that the locus of the end of the vector Ry +jX}, is a circle of 
diameter 2.5/5, passing through the origin and with its centre on the R axis. 
Except very near to resonance, we can disregard 6 and then Xy, assumes the very 


simple form X},=2S/{(w9/w) —(w/@ )}, and this quantity can, in fact, be taken 
in many problems as representing the motional impedance of the instrument. 


§4. ILLUSTRATIVE EXAMPLES 


Two of the most interesting examples of the application of the galvanometer 
equations in the form derived above are to be found in the calculation of the 
inherent ‘‘noise level’? of a galvanometer (that is, the background deflection 
due to spontaneous voltage fluctuations in the circuit) and the calculation of the 
response of an amplifier using a galvanometer. 


(i) Noise level 
From equation (3.5) we may write 
—RI+R~ 


Here e will be the spontaneous E.M.F. of the circuit, and in any frequency range 
df we may write 


De 


aie 
R7\1+R 
where mean square quantities are indicated by.the bars. 
By Nyquist’s theorem (Moullin 1938, p. 13) we may put 
C2 ARR Uidfis 98 a Pate ey (4.2) 


where Rk is Boltzmann’s constant and 7' is the absolute temperature. We 


| D6|? = 
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therefore have, for the mean kinetic energy of the suspended system arising 
from spontaneous fluctuations in the range df, 


a| fer. 
1q| DO\? =4 R2 TER e- 
ae eit = ARRT = 
* R® (w? —wo")? +40?w 20 
_ 2wyckT w - (4.3) 


a  (w*—a,")? aS 
from (2.10) and (2.14). 
The total energy, E, is thus given by 


This may also be set equal to the mean potential energy, $y| 6|?, so that 


[Olt =Ri ly. = 7 eee ere (453) 
Using (2.13) we may write, for the mean deflection 0, 
G=p(nRT/St)t oe (4.6) 


where +=27/w, is the free period of the instrument. ‘This corresponds to an 
R.M.S. current z given by 
i=B)wa(akT/Sr) Oe eee (4.7) 

Ising (1926) has given a similar result (but restricted to the case of critical 
-damping) deduced from considerations of the effect of Brownian movement 
of gas molecules on the suspension. As Moullin (1938, p. 217) has pointed out, 
we must suppose that the total background noise arises from the combination 
of the two effects, and we would expect therefore a background corresponding 
to 14/2, that is, to 


inf 2A (2aRTIST)N = 1 | eee (4.8) 
Taking k= 1-4 x 10-1 ergs per degree and T=290, we have 
i\/2™5x10-7(Sr)-? absolute amperes. —S....... (4.9) 


(ii) Galvanometer amplifier 


If we consider an amplifying system such as that described by Gall ie 
1945), the essential circuit of 
which is shown in figure 1, the 
response can be calculated very 
readily by the present method. 

Referring to figure 1, the 
galvanometer is represented by 
I+r, where r is its D.C. resist- 
ance, Z and Zy, are elements of a 
feedback circuit, and 7; is the 
amplified current arising from 
the application of an E.M.F. e through a resistance P to the galvanometer. 


Figure 1. 
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We then have the following equations 
e= Pi, +(1+r)(2, —15), 
O=(L+7)(t, -—1)) + Zt, + Zy(7, —7g), 
3=g0 5 O=(/RD)(4 —2), 


where g is the gain factor for the amplifying device. Solving for i3, we get for 
the current amplification 


Zo va +Zy, 1 
= ae er ee er ee 4.10 
a 2x 1+é ( ) 
where 
G=RD Ir EZ PLZyligiZye. 9 sen ts (alt) 


The quantity € is a non-dimensional complex number. In any practical 
case it is easy to evaluate, and generally the two components are of the order 


of 0-01. 
The effective input impedance, J’, of the aibemnareiges is given by 
a4 +). L+iy 1 ; 
Pee ieSE NA eae oS 


pp) 


‘Thus in general I’<J+r. 

In Gall’s amplifier, Z and Zy, were pure resistances, Z/Z,, being of the order 
of 104. If Zy be of the form S+ MD, i.e. a combination of mutual inductance 
M and resistance S, and if Z be a large resistance P, we have, disregarding &, 


ts/t)=(P+S+MD)(S+MD), —...... (4.14) 
so that 7, can be made to record the zntegral of 2, over certain ranges of D. It is 


easy to make S/M about 0-001 (say with M=1n, S=0-001 Q) and thus to secure 
an integrating time constant of the order of 1000 seconds (Astbury 1948). 


Effective current gain 
Gain without feedback ° 


(iii) Galvanometer response problems 


The present calculus leads very readily to the operational method for the 
solution of certain response problems. We may take as simple examples the 
problem of determining the response to suddenly applied £.M.F. or flux signals. 
In the former case, let the E.M.F. be H(t).e, where H(t) is the unit operator. 
‘Then from (3.5), we have 

(ae f 1 


Cie oR aay) 


Writing J/(I+ R) IED) we have heals 1931) 


g=1 f(D). rs .H(t).e 


(2 (Cy Sa (4.15) 


1 e ets® e At 

|. es pas eee (4.16) 

For the second case, if the applied ale is H(t).®, we have 
La ae 
ag tiee oun ae 

1 ) C+ jo eft 
- — {ae eee 4.17 
~ kh’ Qny Ae r gaan i) 
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The complex integral in both cases is taken round a curve from ¢ —jao toc +jo, 
c being positive and finite. The evaluation of the integrals is straightforward, 
since f(A) has simple poles at the roots of A? + 20m A + 9" =0. 


(iv) Two galvanometers in series 


Finally, we may note the easy way in which we can write down the solution 
to the problem of two galvanometers in series, an arrangement suggested by 
Butterworth for filtering. If 7, and /, are the motional impedances of the instru- 
ments and R the static resistance of the circuit, the current z is given by 


t=el(Ra1,+1,) ee (4.18) 
and the deflections 6, 0, are given by 
veal 1 es 
O10= had Recieen edt 
and 
gos, es Sha Ree (4.19) 
R+T,+1,|  (R?+4(S1/% +S2/25)?5? 


where the zs are as defined in (3.8). 
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DISCUSSION 


Mr. J. F. W. BELL. I have recently been making use of the equivalent electrical circuit 
concept of the moving coil galvanometer as deduced by Butterworth. When a condenser is 
placed across the terminals there is an increase in the period, as would be expected if the 
galvanometer is looked upon as a parallel LCR network. The change in period gives the 
ratio k?/gamma. From a knowledge of the unloaded period and the current Sensitivity the 
constants a, 8 and k can be calculated. This method of determining the constants of a 
galvanometer has given good agreement with nominal values. 

There are possible applications of the galvanometer as an electromechanical resonator 
at power and low frequencies, similar to the use of the piezoelectric crystal resonator at 
radio frequencies. Its use in an electronically maintained sine wave generator is of 
particular interest. 


Mr. J. C. Wittmotr. I would like to ask Mr. Astbury if the formula @=p VtkT/sT 5 
for the mean deviation is applicable whatever the resistance of the circuit connected across 
the galvanometer. 


AuTHor’s reply. I had originally intended to include as an example in the paper the 
problem quoted by Mr. Bell, and perhaps I may be allowed to show how simply the result 
can be obtained. If r be the D.c. resistance of the galvanometer, then the impedance, Z, 
obtained by connecting a capacitance C in shunt with the instrument, is given by 
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1/Z=1/(r+-D+j@C or Z=(r+D)/[1+joC(7r+D]. The maximum deflection will occur 
when jwCI+1=0, from which w?= w,2/(41+ 2SCw 9). Here w» is the natural frequency of 
the instrument and w the frequency as modified by the presence of the capacitance. 

I feel that measurements of p, S and 9 Should be sufficient to provide estimates of all the 
other constants required, although I might mention that during the war I did devise a 
method for testing galvanometers based on a potentiometric method for measuring Ry and 
Xv at various frequencies. I agree with Mr. Bell that applications do exist for the galvano- 
meter as he suggests, and I would like to add that it is good to see an appreciation of the 
galvanometer as a scientific instrument and not merely as a hack for null-detection work. 

In reply to Mr. Willmott, the analysis given in the paper certainly leads to the result 
that the noise level is independent of R. We have made records of background signals 
which show that the level is substantially independent of R for values between critical 
resistance and one-fifth of this value, but the level observed was always higher by about 50% 
than that given by equation (4.9). 


FE BEER SSO TEE DLO 


The Structure of Carbon Monoxide 


In both the molecular orbital and pairing approximations it is easily shown that the 
four 7 electrons of carbon monoxide are completely equivalent, so that the distinction 
between bonding and non-bonding electrons recently introduced by Long and Walsh (1947) 
is invalid. The symmetry of the molecule is most simply exhibited by Mulliken’s formu 
lation (1932) : 

(K)(K)(20)"(yo)*(wor)*(xo)?, 14, 


where (wz) is doubly degenerate. 

The existence of the (w7)* shell, while indicating that four equivalent 7 electrons are 
concerned in 77 bonding, by no means implies that the molecule is best represented by a 
triple-bonded structure, as Long and Walsh seem to suppose. Indeed, analysis shows 
(Moffitt, to be published) that of the conventional structures 


C=OFIF C—Onulr C—Onrlit 


II is the most important, and I is less important than III. 

The rotational analysis (Schmid 1932) of the (CO)* produced by the first ionization of 
the molecule has shown it to have the symmetry 7%. Consequently the first ionization 
potential of CO, 12+ corresponds to the removal of a o electron. ‘The removal of a 7 
electron, as postulated by Long and Walsh, would lead to (CO)*, 7iJ, and Mulliken (1928) 
has cited evidence to show that this ionization occurs near 17 ev. and not near 14 ev. 

It is difficult to follow Long and Walsh in their reasoning that the removal of a 7 
electron should lead to a strengthening of the bonding in the resulting (CO)*. Introducing 
the wave-function symmetrization necessary to give significance to their structure, the ioni- 
zation of an electron which accounts for a quarter of the considerable 77 bond strength in 
CO would surely lead to a weakening; and this is actually observed for CO*, a*II (Sponer 
1935). The ‘‘ three-electron bonding function” involving one 7 and two o electrons, 
invoked in this connection, has, as far as I am aware, no precedent. 

The low dipole moment of CO is well known ; calculations (Moffitt, to be published) 
based on Long and Walsh’s structure indicate, however, that the molecule should be highly 
polar. Further, the electron donor properties of CO are hardly to be attributed to the 
“inert pair” (2s), whose ionization potential is expected to be near 20 ev., as assumed by 
these authors. a 

The real difficulty in describing the carbon monoxide molecule does not lie in the 
specification of the LCAO forms for the (zw7)* shell, but in characterizing the nature of the 
three different o-type orbitals. A plausible assignment of the Six o electrons to there 
localized molecular orbitals has been devised (Moffitt, to be published). Retaining Mulli- 
ken’s notation, (zo)? now corresponds essentially to (2s9)” and (ya)” refers to an electron pair 
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| 
bond between some sp-type hybrid carbon orbital and a predominantly 2p¢ oxygen orbital. | 
The third pair of electrons (xc)? are localized near the carbon atom and are nearly non- 
bonding. ‘The ionization potential of a (xc) electron is estimated at near 14 ev., and, there- 
fore, accounts for the first ionization potential and electron donor properties of carbon 
monoxide. The (yo)? electron pair bond is shown to be largely homopolar. 

The effect of the s-p hybridization gives rise to a compromise between strong oo-type 
bonding and the promotion required to produce the appropriate carbon valence state. 
This compromise has been shown to lead to an asymmetry in ¢o electron distribution, whose 
direction is opposite and whose magnitude is approximately equal to that of the 7 electrons ; 
the low dipole moment of carbon monoxide is no longer surprising. 

The ionization of a non-bonding (xe) electron leads to a strengthening of the 77 bonding 
in as much as (CO)+, x2% now becomes more nearly triple-bonded in the conventional sense. 
Further, the properties of the excited a‘ll state of CO (Sponer 1935), in particular its weak 
bond strength, are accounted for by the formulation 


(K)(K)(2’0)?(9'0)"(wo'7)*(x"o)(w'7), *H, 
where (zw’z) is the antibonding orbital corresponding to (zw’z). A full account of this work 
is being prepared for publication. 
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A Note on the Shape of the Polepieces of a 
Synchrotron Magnet 


The Birmingham synchrotron magnet has recently been described (Oliphant, Gooden 
and Hide 1947), and an interesting problem connected with its magnetic field has arisen. 
For certain reasons it is desirable that the magnetic field strength H along the central 
radius of the gap between the polepieces of a synchrotron ring magnet should, in any 
normal section of the ring, be related to the radius vector r drawn from the centre of the 
ring by an equation of the form H=ar~-*. What shape should the polepieces have to give 
such a field? If the direct field due to the current in the coils be neglected, an approximate 
solution of the problem can be obtained by regarding it as a two-dimensional one. This 
is justified by the large mean radius of the ring, which is 450 cm. in the Birmingham 
magnet. : 

Once the method of conformal transformations has been adopted, the fact that the 
integral of x—$ dr is equal to 374 suggests the trial of 


WaAst 2. | eee (1) 


aS a possible solution. A is a constant and W is the complex potential of the form 
W=U-+1V and z=x+iy=re®. Applying equation (1), we get 


W = Arrel3 = Art cos 6/3+1Ar sin 6/3. 
Take V as the potential and U as the stream function, then, since 
U=Ar* cos 6/3 ; V=Ar* sin 0/3, 
equipotentials are given by 
rt sin 0/3=B or. r=C cose? 6/3. ~~ =) Gea (2) 
and lines of force are given by 


r> cos 6/3=D or r~Esec3 6/3. BA oe (3) 
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Figure 1 represents sets of such curves. 


Figure 1. Section of field of “‘ ideal”? synchrotron magnet (right side only). Lines such as 
B,B.B;, CyC,C3, are equipotentials and possible polepiece profile curves. Lines such as 
B,A,C,, BzAsC, are lines of force. Beams of particles pass normal to the plane of the figure 
through points on A,A,A3. 


The radial field strength H,=—2aV/ér=—(Ar-? sin 0/3)/3, and the tangential 
component Hy= —@V/r00= —(Ar-* cos 6/3)/3. The resultant field strength is H= — Ar~3/3. 
Along the central radius 6=0; thus along that line H,=0, Hp =H=-— Ar-?/3. Hence 
the field satisfies the desired condition. At any point H,/Hj=tan 6/3, showing that for 
small values of 0, H, is very small compared with Hp. 

The polepieces are to be made with profile curves y=C cosec? 0/3, for they are them- 
selves magnetic equipotentials. 

The complete curves given by (2) and (3) consist of closed loops with “ wings ”’ extending 
to infinity. The two sets of curves are of the same shape but the lines of force are obtained 
by rotating the other set through 270° about an axis through the origin normal to the plane 
of the figure. The necessity of making the real polepieces finite in length instead of 
infinite increases the discrepancy between theory and experiment, but in the middle of 
the gap the difference should not be great. 

The ‘ basic”? shape of the profile curves of the Birmingham synchrotron magnet, 
obtained empirically by the aid of the Peierls tank, appears to be very close to that suggested 
by equation (2). 


University College, Nottingham. N. Davy. 


23 February 1948. 
O.ipHant, M. L., Goopsn, J. S., and Hipe, G. S., 1947, Proc. Phys. Soc., 59, 666. 


REVIEWS OF BOOKS 


The A Priori in Physical Theory, by ARTHUR Pap. Pp. x+102. (New York: 
King’s Crown Press, 1946.) 13s. 6d. net. 


This book deals with a problem which, mainly through the assertions of the late Sir 
Arthur Eddington and Professor E. A. Milne, has in recent years become of much interest to 
physicists, namely the question whether or not the laws of nature are a priori, i.e. derivable 
by reason apart from experience and so not liable to falsification by experience. ‘The author 
does not, however, attempt to solve the ultimate problem, but points out that in practice 
physical laws are often treated as “ functionally ’’ a priori ; thus, if an experiment or 
observation appears to violate a law we sometimes do not abandon the law but regard the 
experiment as not providing a fair example of it. Thus, when Uranus appeared to violate 
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Newton’s law of gravitation, the law was not changed but Uranus was held to be disturbed. 
Laws sufficiently well grounded can acquire a high degree of immunity in this way. It was 
on such grounds that Poincaré maintained that one could always regard space as euclidean by 
attributing any apparently non-euclidean qualities to fields of force or non-rectilinear 
propagation of light. The book is clearly and interestingly written, though one feels that 
the author tends to continue his argument after the point has been adequately explained. 
With the main thesis it is impossible to disagree, for it simply records what every 
physicist knows to be true. It inevitably raises another question, however, which the 
author does not consider—namely, what determines the choice between abandonment of the 
law and rejection of the observation as irrelevant ? This is a pertinent question because, for 
example, while the anomalous behaviour of Uranus did not upset Newton’s law of gravita- 
tion, the anomalous behaviour of Mercury did so. The difference was not merely that 
Neptune was observed and “ Vulcan ”’ was not, for the law was retained in spite of the 
irregularities of certain stars to which were assigned purely hypothetical companions. It 
would be interesting to see how the author would deal with this extension of his problem. 
Mr. Pap is evidently better acquainted with physics than some philosophers who write 
on physical theory, and his views command respect. There are a few slips, however, 
which might be worth mentioning. The principle of conservation of mass is not rendered 
invalid by the discovery that mass depends on velocity (p. 61) ; it is, in fact, in order to 
maintain that principle that we must make mass depend on velocity. Again, neither in 
Newton’s nor in Einstein’s language is it true to say that force means m du/dt (p. 65) ; it 
means d(mv)/dt. An injustice is done to Galileo on p. 85 where it is said that he regarded 
the string of his pendulum as having no mass. He considered the effect of its mass very 
interestingly in the Dialogues Concerning the Two Great Systems of the World. Finally (p. 22), 
the author, apparently following Lenzen, seems to think that there is a “‘ vicious circle ”? in 
making the length of a standard measuring rod depend on its temperature and defining a 
temperature scale in terms of length—a circle which needs a series of “‘ successive approxi- 
mations ”’ to break it. This is not so, however, for a thermoscope is sufficient to make the 
prescription of the standard measuring rod definite, leaving us free to choose whatever 
type of thermometer we wish. Nor, in fact, do we ever measure temperature in terms of 
_measures of length ; the marks on eur glass tube containing mercury are arbitrary so far as 
the unit of length measurement is concerned, which is why it is possible (though not 
necessary) to regard length and temperature as independent fundamental measurements 
in dimensional theory. HERBERT DINGLE. 


Five-figure Tables of Natural Trigonometrical Functions, prepared by H.M. 
Nautical Almanac Office. Pp. 124+iv. (London: H.M. Stationery 
Office, 1947.) 15s. net. 


Tables of the natural values of the four trigonometrical functions: sine, tangent, 


cotangent, cosine, for every 10 seconds of arc, and an auxiliary table for the cotangent for 
every second of arc to 7° 30’. ACerS 


Logarithmic and Trigonometric Tables, by J. B. Date. Pp. 42+vi. (London: 


Edward Arnold and Co., 1947.) 2nd Edition. 2s. éd. 


A new edition, completely reset in very clear form. ANCaSs 


Note. Details of binding cases will be given in the July issue of the Proceedings. 
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